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In materials with strong local Coulomb interactions, simple defects such as atomic substitutions 
strongly affect both macroscopic and local properties of the system. A nonmagnetic impurity, for 
instance, is seen to induce magnetism nearby. Even without disorder, models of such correlated 
systems are generally not soluble in 2 or 3 dimensions, and so few exact results are known for 
the properties of such impurities. Nevertheless, some simple physical ideas have emerged from 
experiments and approximate theories. Here, we first review what we can learn about this prob- 
lem from ID antiferromagnetically correlated systems. We then discuss experiments on the high 
Tc cuprate normal state which probe the effect of impurities on local charge and spin degrees 
of freedom, and compare with theories of single impurities in correlated hosts, as well as phe- 
nomenological effective Kondo descriptions. Subsequently, we review theories of impurities in 
d-wave superconductors including residual quasiparticle interactions, and compare with experi- 
ments in the superconducting state. We argue that existing data exhibit a remarkable similarity 
to impurity-induced magnetism in the ID case, implying the importance of electronic correlations 
for the understanding of these phenomena, and suggesting that impurities may provide excellent 
probes of the still poorly understood ground state of the cuprates. 
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I. INTRODUCTION 

A. Aim and scope of this article 

The study of "strongly correlated electron systems" . 
materials where the typical electronic interaction ener- 
gies are of order or larger than the kinetic energies, un- 
derwent a renaissance in the 1980's and 1990's after the 
discovery of superconductivity in organic conductors, in 
heavy fermions and then in copper oxides. In particu- 
lar the cuprate family has had an enduring fascination 
for the condensed matter community, and not simply be- 
cause of the allure of room temperature superconductiv- 
ity. Most of the important unsolved problems of quan- 
tum condensed matter physics - metal insulator transi- 
tion in low dimensions, breakdown of Fermi liquid thcjory, 
origin and behavior of unconventional superconductivity, 
quantum critical points, electronic inhomogeneities and 
localization in interacting systems, and quantum antifcr- 
romagnetism in low dimensions-are represented in the 
set of questions currently under active consideration in 
this lively field. Two factors are currently believed to 
underly the anomalous properties of these materials: low 
dimensions and strong correlations, both of which appear 
to be necessary to stabilize the new types of electronic 
states present. 

The cuprates contain as common elements Cu02 
planes (Fig.l), which are sometimes considered to con- 
tain all the important physics. Their structure is a stack- 
ing of such planes separated by other oxide layers, which 
maintain charge neutrality and the cohesion of the struc- 
ture essentially through ionic interactions. High tem- 
perature superconductivity (HTSC) has been of course 
the physical property which has driven the interest for 
these systems initially both for its fundamental aspects 
and its potential applications. Another aspect respon- 
sible for their appeal has certainly been the ease with 
which the carrier concentration can be changed by chem- 
ical hcterovalent substitutions or oxygen insertion in the 
layers separating the Cu02 planes, which play the im- 
portant role of charge reservoirs. Electron or hole dop- 
ing can then be fairly continuously controlled from zero 
to about 0.3 charges per unit cell, which allows one to 
study the evolution of the physical properties of these 
materials with doping. Despite the fact that the phase 



diagrams of electron- or hole- doped materials look quali- 
tatively similar, their electronic properties in the normal 
and in the superconducting states show significant dif- 
ferences. To date, the hole doped family has been more 
extensively studied, particularly as far as the issue of im- 
purities is concerned and in this review we chose to focus 
on this type of systems. In most other classes of corre- 
lated electron materials such as quasi -ID or -2D organic 
compounds, ID chains, or fuUcrides, the number of car- 
riers is difficult to modify, except close to some specific 
compositions, where pressure must be applied to modify 
electronic properties. 

Cu O 

-Ot*-tH>-CH>0-oG^>-0^>G- SrjCuOj 




FIG. 1 A selection of actively studied inorganic oxides. The 
figure only shows the chain (Sr2Cu03, Y2BaNi05, CuGeOs), 
ladder (SrCu203, Sr2Cu305) or planar (cuprates) unit in the 
molecular structure. 

Impurities were used initially to test their pair- 
breaking influence on superconductivity. The observa- 
tion that nonmagnetic impurities such as Zn reduced 
drastically triggered a further interest in the electronic 
properties of impurity states in the normal metallic 
phases as well, and led to the discovery that even non- 
magnetic impurities induce local magnetism, the central 
theme of this review. In the following introductory para- 
graphs, we highlight general features which character- 
ize the physics of correlated electron systems, illustrating 
them using the phase diagram of the cuprates. We then 
recall why point defect impurities can be used as spe- 
cific probes of the physical properties of these complex 
materials, where understanding the role of disorder has 
proven so crucial. 



1. Correlated electrons systems at a glance 

Generically at low T an insulating, magnetic phase 
is found in some region of the phase diagram of corre- 
lated systems, while a metallic or superconducting state 
is present in a neighboring (or in some cases overlap- 
ping) region. This situation can be illustrated with the 
cuprate phase diagram of Fig. 2, which displays the dif- 
ferent phases versus the hole content in the Cu02 plane. 
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Notice that here zero doping corresponds to half fiUing, 
i.e. a situation where one hole is localized on each Cu 
atom, yielding an AF insulating ground state. Such a 
state is usually nonmctallic when the Coulomb interac- 
tion is large enough with respect to the bandwith. In 
some easily compressible materials, for instance organic 
systems such as the k-(ET)2CuX salts, applying pres- 
sure increases the bandwidth and yields a Mott insulator 
to metal transition. In the case of cuprates and some 
heavy fermion compounds, such as CeRhi_^Ira;In5, the 
AF order is suppressed by doping and is replaced by a 
metallic or a superconducting state. These types of metal 
insulator transitions (MIT) or superconductor insulator 
transitions (SI) are archetypes of the problems encoun- 
tered in correlated electron physics, which are not yet 
fully understood. 

The phase diagram of correlated electron systems then 
presents many novel states of matter, which in the case 
of cuprates includes a spin glass, a pseudogap regime, a 
"strange metal" , a d-wave superconductor and a putative 
Fermi liquid, with only the last thought to be reasonably 
well understood (see Fig. 2) theoretically. The nature 
of these phases has been reviewed elsewhere and will be 
briefly discussed in Section IV below. At T = the 
boundaries between some of these phases may correspond 
to Quantum Critical Points (QCP). In the cuprate phase 
diagram of Fig. 2 QCP might be present at the end points 
of the superconducting dome, or at the extensions to T = 
of the dashed lines shown in Fig. 2, which have however 
not yet been established as true phase transitions. 
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FIG. 2 Schematic phase diagram of cuprate HTSC materi- 
als: temperature T versus hole doping Uh- AF refers to a 
state with long-range antiferromagnetic order, which disap- 
pears at the Neel temperature Tjv. Superconductivity with 
d-wave symmetry sets in below a critical temperature Tc. 

Much can be learned from studies of systems with di- 
mensions less than two, since many theoretical results are 
available for ID and quasi-lD systems. Furthermore, it 
is natural to imagine approaching the 2D system, where 
few exact results are known, by adjoining chains to form 
ladders and to successively increase the number of legs. 
Many physical situations encountered in one dimension, 
such as the separation of spin and charge degrees of free- 



dom yielding the concepts of spinous and holons, which 
prevail in ID can be conceptually followed if one progres- 
sively increases the dimensionality of the system from 1 
to 2. Although this was initially considered as a purely 
theoretical approach, the skills of researchers in materials 
science have allowed the synthesis of systems implement- 
ing the theoretical strategy. 

Indeed, inorganic realizations of quasi-lD spin chains 
are represented, e.g. by Sr2Cu03 (spin-1/2), Y2BaNi05 
(spin-1), and CuGeOs (spin-1/2, spin-Peierls). The 
Sr2Cu03 material contains CuO chains which are sim- 
ilar to those found in the 2D cuprate YBCO. The ladder 
compounds consisting of coupled spin-1/2 CuO chains, 
e.g. SrCu203 (2-leg) and Sr2Cu305 (3-leg) are natural 
low-D analogs of the cuprates (see Fig. 1). We discuss 
these systems in Sec. III. 



2. Impurities as probes of correlated electron systems 

To gain insight into the complex problems which arise 
in correlated electron systems, all possible tools and tech- 
niques have been extensively used and improved. We 
shall show below that a specific approach, given a privi- 
leged place in this review, has been to study the modifica- 
tions of the physical properties induced by the controlled 
introduction of point defects in the materials. We shall 
report studies of impurity substitutions, of defects cre- 
ated by electron irradiation, or of defects which occur 
naturally even in the purest available materials synthe- 
sized so far. This approach combines material science 
investigations with the use of the most advanced local 
probe techniques. 

Impurities are known to modify superconductivity in 
crucial ways. In conventional s-wave BCS materials, dis- 
order causes changes in the penetration depth A and in 
the coherence length ^. For superconductors with non s- 
wave symmetry, very different effects are expected since 
impurities break Cooper pairs. It is therefore quite nat- 
ural to probe the influence of defects on most matters 
linked with superconductivity. 

In conventional metals, impurity effects constitute a 
relatively well understood subject. The essence of the ob- 
served phenomena is that an impurity is a local screened 
Coulomb potential, which ideally is a uniform perturba- 
tion in q space, inducing a response which is inhomo- 
geneous in real space but which reflects the response to 
all q values. So if some singularity occurs at a specific 
q value for a system, this singular response will domi- 
nate the modifications introduced by the impurity. In a 
classical metallic system, since the response is homoge- 
neous up to |q| = kp, the main detected feature comes 
from this truncation of the response at fc^?, which yields 
the well known Friedel oscillations in the local density of 
states for a charge defect, and the RKKY oscillation for 
a spin defect. This idea can be naturally generalized to 
any type of material. For instance, in magnetic materials 
for which AF correlations can be important at a wave 
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vector q = qAF, a staggered magnetic response at this 
wave vector is expected (Fig. 3). This is exemphfied in 
the case of an AF spin chain, in which a staggered para- 
magnetic response appears near chain ends. So quite 
generally an impurity potential is a fundamental tool to 
probe the specific response of a system. 




FIG. 3 Cartoon of the staggered paramagnetic polarization 
induced by a nonmagnetic defect (gray dot at center) in a 
correlated 2D material. Such oscillations of the local suscep- 
tibilities are directly detectable by spin probe experiments 
such as NMR, since the NMR shifts of the various shells of 
neighbors of the defects can be resolved. 

The potential importance of this property was cap- 
tured by experimentalists who developed local probe 
techniques, such as NMR and STM, to map out the 
spatial changes occurring around extrinsic defects. For 
instance, NMR experiments performed on the cuprates 
over the past fifteen years or so, using a variety of nuclei, 
have developed a fairly clear picture of how a nominally 
"nonmagnetic" impurity creates a cloud of local stag- 
gered polarization on nearby Cu sites, and have charac- 
terized the response as a function of doping and temper- 
ature, for different materials. 

These striking results triggered a growing interest in 
the theory community as well. The enormous theoretical 
effort devoted to solving problems of homogeneous 2D 
correlated Fermi systems by approximate methods has 
been occasionally applied as well to impurity problems 
or even to finite disorder. While at this writing there is 
no generally accepted theoretical description of an im- 
purity in a strongly correlated 2D system (as there is 
no generally accepted theory for the homogeneous prob- 
lem) , a qualitative picture captured by many different ap- 
proaches has emerged and will be summarized in this re- 
view. By contrast, the impurity problem has been quan- 
titatively studied in the case of ID systems, since exact 
analytical techniques are available, and the defect's influ- 
ence on transport, magnetism and thermal properties is 
known. The best studied and understood example of this 
phenomenon is the nonmagnetic impurity in spin chain 
systems, which is equivalent to a chain end in a non in- 



finite spin chain. Here it can be shown in a controlled 
fashion that a single (even weak) potential defect induces 
a paramagnetic net moment and staggered magnetism in 
its vicinity, and that the decay of this magnetism di- 
rectly measures the spin correlation functions of the pure 
state. Observations done on 2D cuprates have been suc- 
cessfully extended to quasi- ID systems where impurities 
have been used to probe magnetic correlations in a wide 
variety of materials. The behavior of these states depend 
on whether the pure host system manifests a spin gap, 
which in turn depends on the number of coupled chains 
and whether the spins are integer or half-integer. 

Studies of impurities in correlated electron systems are 
also important to help us understand the influence of 
(intrinsic) disorder that is present in nominally pure sys- 
tems. As pointed out before, in cuprate materials and in 
heavy fermion systems, one can easily change the carrier 
concentration and thus explore a wide range of the phase 
diagram using the same system. However the very dop- 
ing process introduces disorder. The importance of this 
phenomenon has only been recently appreciated by much 
of the community, as local techniques such as STM have 
imaged impurity states directly and brought to the fore 
the study of disorder as a key issue in correlated elec- 
tron systems. The efforts which have been made so far 
to understand the influence of controlled disorder and 
impurities are very good reference points to try to un- 
derstand those due to the intrinsic disorder. So although 
this is not the main aim of the present review, we shall 
try hereafter to recall as well some of the established 
facts which should help to separate out in the materials 
problems those which might be attributed to the intrinsic 
disorder. 



B. Structure of the manuscript 

We will first give some preliminary background in Sec- 
tion II, discussing the basic ideas of strongly correlated 
electron physics, introducing the simplest and most pop- 
ular models like the Hubbard and t — J Hamiltonians, 
along with some brief words about current methods em- 
ployed to solve them. We then begin our main theme 
in the strange world of one dimension: in Section HI we 
summarize progress which has been made in recent years 
on the conceptually simplest problem in the field, that 
of an impurity in a ID spin chain. We first attempt to 
explain qualitatively the physical phenomena expected, 
discuss the theoretical background for this expectation, 
and justify the use of NMR to study such effects. We de- 
scribe the most studied materials and models first in their 
idealized homogeneous form in terms of their physical 
properties (Haldane gap, quasi-long range order). In the 
disorder- free limit, one finds correlated quantum ground 
states, called spin liquids, and we define the concepts of 
Resonating Valence Bond (RVB) phase, and of spinous 
or holons excitations. These notions will be quite useful 
when we turn to the 2D case. We then introduce mod- 
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els of defects and the concept of end-chain excitations, 
and show why NMR is particularly useful for studying 
such effects. The discussion is then broadened to include 
ladder compounds and models, which are in certain well- 
defined senses steps on the road to two dimensions. For 
the most part we restrict ourselves to pure spin models, 
as there are few relevant experiments on doped ladder or 
chain compounds. 

In passing to two dimensions in Section IV, we fo- 
cus on the cuprate family, and first discuss the various 
types of point defects which occur or can be induced in 
these materials. We then attempt to make contact with 
what we have learned from spin chains, by focussing on 
the undoped antiferromagnetic cuprate insulators. The 
analogy is not completely straightforward, however, since 
although the cuprate planes are sometimes regarded as 
idealized 2D Heisenberg spin systems, the interplane cou- 
plings and small spin anisotropics must be accounted for 
as well. In addition, few experiments have been per- 
formed on local properties of defects in the insulating 
systems, and so we pass quickly to the underdoped pseu- 
dogap state, in the metallic regime. Here the suppression 
of the single-particle DOS and magnetic susceptibility 
suggest that impurity-induced magnetism should have a 
large relative effect, and indeed the first observation of 
moment formation around impurities was made in un- 
derdoped YBCO. We review these experiments, the mea- 
surement and interpretation of the temperature depen- 
dence of the local susceptibilities. We compare the sim- 
ilarities and differences between impurities, in the cases 
of Zn, nominally nonmagnetic, and Ni, nominally spin 
1, and Li, nominally nonmagnetic but charged. We re- 
view the changes in local susceptibility found to occur 
when the system is doped further to optimal doping, and 
discuss the behavior of the impurity Curie- Weiss tem- 
perature, and a simple interpretation in terms of Kondo 
screening. Finally, we discuss other measurements which 
do not give local information but which are relevant to 
understanding the physics of disorder in these materials, 
particularly transport measurements. 

In Section V we review theoretical efforts to under- 
stand the remarkable phenomenology just described in 
higher dimensions. Some of the ideas and relevant cal- 
culations about moment formation in correlated systems 
are not new, and we therefore begin by reminding the 
reader of early work on nearly magnetic alloys and phos- 
phorous doped silicon. More recent work related to the 
properties of Zn and Ni in cuprates which describes these 
ions as localized potentials embedded in a 2D Hubbard 
treated within the random phase approximation are the 
intellectual heirs of these early works. We describe the 
basic ideas, which involve the buildup of antiferromag- 
netic correlations on the lattice due to the Hubbard re- 
pulsion, and the mode-mode coupling induced by the im- 
purity between this staggered response and the homoge- 
neous response actually measured in NMR in the normal 
state. These theories provide a remarkably good descrip- 
tion at optimal doping, but it is not clear how to extend 



them to describe the effects in underdoped samples. We 
then discuss several theoretical approaches which begin 
at strong coupling with a t — J model, and then discuss 
moment formation in the spin gap state. For the most 
part, such theories have described spontaneous moment 
formation, enhanced by the presence of a pseudogap, and 
do not discuss the paramagnetic state in detail. 

Section VI tries to understand the impact of supercon- 
ductivity on the effects described above. To place the 
discussion in context, we begin by reviewing briefly the 
basic ideas concerning the properties of a single impurity 
and many-impurity interference in a d-wave supercon- 
ductor. We then describe how the predictions of theories 
described in Section V are modified in the presence of a 
gap. Finally, we discuss experiments on impurities in the 
superconducting state which directly probe local density 
of states and magnetic correlations and we attempt to 
compare the results to these ideas and predictions. 

II. MODELS OF CORRELATED FERMI SYSTEMS 

To set the stage for our discussion of defects, we first 
give a lightning account of the physics of homogeneous 
correlated systems and the models used to describe them. 
We refer the reader to more extensive reviews on the 
subject (Fazekas, 1999; Fulde et al., 2006). Prom a theo- 
retical perspective, transport, thermal and magnetic ob- 
servables can be obtained in strongly correlated systems 
once one has constructed a Hamiltonian containing all 
the essential ingredients characterizing these systems, yet 
simple enough to be amenable to analytical or numerical 
computations. By and large, the model of choice is the 
Hubbard model, which was introduced in the early sixties 
to describe magnetism and metal-insulator transitions in 
narrow band transition metals ((Gutzwiller, 1963a; Hub- 
bard, 1963; Kanamori, 1963)). It consists of a kinetic 
term describing the hopping dynamics of electrons (or 
holes) with site creation operators a'^i,^ and a potential 
term accounting for on-site Coulomb interactions. 

'HHubb = X] ^^'i ^J.<^ + ^•^) + X] '^i.'^^'Ml) 

In the simplest version, the hopping amplitude ti^j is lim- 
ited to nearest-neighbors and a is one of the spin-i eigen- 
states, while 'W stands for the other one. Ui = ni^„ 
fermions occupy site i. The case [/ > corresponds to 
repulsive interactions, while [/ < mimics effective at- 
tractive interactions (mediated for instance by electron- 
phonon or electron-paramagnon processes). For the spe- 
cial case J7 = the state is either metallic or (band) 
insulating, depending on the filling factor. The Hub- 
bard model was introduced originally to study the Mott 
transition which occurs when the bandwidth of a system 
narrows sufficiently in the presence of repulsive U. The 
basic physics can be understood by considering a lattice 
of atomic orbitals which can overlap to allow conduc- 
tion. If there is one orbital per site, and one electron 
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per orbital, the system is a metal if U = according 
to conventional electronic structure theory. Thus elec- 
trons may hop from orbital to orbital, provided all states 
satisfy the Pauli principle allowing only one t and one 
J, spin per orbital, i.e. the Hilbert space for each site 
consists of |0), | t), I i), and | ti)- As C/ > is in- 
creased, however, the Coulomb cost of doubly occupied 
sites I ti) increases, such that in the limit U +oo such 
configurations are no longer allowed, or are "projected 
out" of the Hilbert space. If no doubly occupied sites are 
present, electrons are effectively prevented from hopping 
in the half-filled case, and if a small number of holes are 
introduced one can easily see that, although conduction 
can now take place, it is drastically reduced due to the 
Coulomb repulsion. This is the primary effect of elec- 
tronic "correlations" which will be considered here. 

More generally, if U is non zero, depending on its sign, 
on the size of the ratio U /W {W is the bandwidth) and 
on the filling factor one expects to promote various insta- 
bilities of magnetic, metal-insulator or superconducting 
types, or even phase separation. These transitions may 
occur at finite temperature or at zero temperature de- 
pending on the spatial dimensionality. Only a few exact 
results are available for this model, mostly in the one- 
dimensional case (Eggert, 1993). In that limit, fluctu- 
ations prevent the occurrence of ordered phases at fi- 
nite temperature. At T = and for half filling, the 
ground state is insulating and has quasi-long range anti- 
ferromagnetic (AF) order when U > 0, or on-site singlet 
ordering when U < (Auerbach, 1994; Giamarchi, 2004). 
By quasi-long range order we mean that the expectation 
value of the staggered magnetization is zero in the ther- 
modynamic limit, but that the corresponding coherence 
length is infinite (the alterning part of the spin-spin corre- 
lation function decays as a power law at large distances) . 
For any other band filling, the chain is a non Fermi liquid 
paramagnetic metal, when Coulomb interactions are re- 
pulsive, or a quasi ordered s-wave superconductor when 
interactions are attractive. For higher dimensions, only 
a handful of methods allow analytic treatments of the 
Hubbard model. The simplest is the Hartree-Fock ap- 
proximation, including RPA extensions and beyond to in- 
clude fiuctuation effects (Aryanpour et aZ., 2002; Bickers 
et al., 1989; Hertz, 1976; Hirsch, 1980; Kanamori, 1963; 
Pao and Bickers, 1995). Its validity corresponds essen- 
tially to the weakly correlated regime {U/W much less 
than one). When U is large and positive, one may use 
the Gutzwiller variational wave function to determine the 
T = properties of the model (Gutzwiller, 1963b), the 
Hubbard X-operators (Ovchinnikov and Val'kov, 2004), 
or alternatively auxiliary field methods developed for in- 
stance by Kotliar and Ruckenstein (Kotliar and Rucken- 
stein, 1986) and by Lee and Nagaosa (Lee et al., 2006), 
based on slave bosons (Barnes, 1976, 1977; Bickers, 1987; 
Coleman, 1984; Read, 1985; Read and Newns, 1983). 
Note however that the Gutzwiller wave function accounts 
for on-site correlations but neglect inter-site correlations, 
such as the Heisenberg AF interactions, which play an 



important role in the large U limit (Kaplan et al, 1982). 
In the d = oo limit, many results pertaining to the Hub- 
bard model can be readily obtained by a diagrammatic 
treatment, thanks to the fact that many terms vanish 
(Metzner and Vollhardt, 1989), and the predictions are 
in agreement with those obtained by Gutzwiller and by 
Kotliar-Ruckenstein. Closely related is dynamical mean 
field theory, a combination of analytic and numerical 
tools allowing quantitative and accurate determinations 
of metal insulator transitions (Georges et al., 1996; Kake- 
hashi, 2004; Kotliar et al, 2006; Kyung et al., 2006). 

In the U = +00 limit, double occupancy of the sites 
becomes energetically unfavorable, and carriers perform 
only virtual hops onto neighboring sites. A canonical 
transformation then maps Hamiltonian Eq.(l) onto the 
so-called t—J model (Chao et al., 1977; Harris and Lange, 
1967; Hirsch, 1985) 

Tit- J = ^ U,j(a^i,a aj,a + h.cj 

<i,j>,cj 

+jY,(Si-^j- \ni n,) (2) 

Here, creation and annihilation operators are not the 
canonical ones, but are projected onto the subspace of 
non doubly occupied states. The exchange J = At^ /U 
is the AF coupling between the spins. In the half filled 
case, the kinetic term has to be zero, and one is left with 
a Heisenberg Hamiltonian for localized spins. As already 
mentioned, in ID one gets quasi-long range AF order at 
T = 0. After much debate, it is now generally agreed that 
the ground state is antiferromagnetic in 2D at T = 
(fluctuations destroy the order at finite temperatures). 
Away from half filling, one may treat the kinetic term as 
a perturbation (Brinkman and Rice, 1970; Kampf, 1994). 
The bandwidth is strongly narrowed due to the interac- 
tion between the carrier motion and spin dynamics, and 
the energy scale is not governed hy t, as one might ex- 
pect, but rather by J (Xu et al., 1991). This remarkable 
property is also seen in exact diagonalization of finite 
size clusters (Dagotto, 1994; Leung and Gooding, 1995) 
and experiments are consistent with the theoretical pre- 
diction (Damascein et al, 2003; Wells et al, 1995). We 
note that the Hamiltonian (2), if derived from Eq.(l), 
should also include an additional kinetic contribution of 
order t^ /U involving three sites < ijk > 

-ah,cT ak,aa^j,a aj^a ) (3) 

For non zero doping and when carriers are mobile, one 
can rewrite Hamiltonian Eq.(3) in terms of effective 
Heisenberg interactions involving second and third neigh- 
bor spin couplings (Inui et al, 1988). A wide variety of 
phases can be obtained when the model - in its Hubbard 
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or t — J form - is generalized to include more compli- 
cated kinetic or correlation terms, or else to include sev- 
eral orbitals or hybridizations between different bands. 
Restricting Coulomb interactions to on-site terms is a 
rather drastic approximation, and several authors have 
examined the influence of an additional coupling between 
nearest neighbors in Hubbard chains and ladders (Fab- 
rizio, 1993; Orignac and Giamarchi, 1997); in addition 
to Hamiltonian (1), where (in the case of ladder com- 
pounds) the hopping amplitude tij can take on different 
values along the chains and on the rungs, one includes 
an extra Coulomb contribution: 

'Hnnc = ^ X] '^i'^i+P i^) 

Here p is the neighbor of i on the same chain. For chains, 
this term stabilizes insulating commensurate phases for 
less than half filling. For ladders, it promotes an inter- 
esting orbital antiferromagnetic order. In organic ma- 
terials, the molecular bonding is responsible for a pro- 
nounced quasi-lD character (for instance in the Bech- 
gaard or Fabre salt series (Bourbonnais and Jerome, 
1994)), or a quasi -2D band structure character (for in- 
stance in the BEDT scries (Limelettc et ai, 2003)). In 
the former case, the hopping amplitude Uj in Eq.(l) 
has a value ta — lOOOK along the conducting chains, 
tb — ta/10 in one of the directions perpendicular to the 
chains and tc ~ ta/ Wi) in the third direction. In the lat- 
ter case ta — tb- Since doping has not been successfully 
achieved for these materials, the main tuning parameter 
is pressure: it is either applied externally or internally 
when (halogen) atoms arc chemically substituted into the 
structure and this changes the ratio t/U. 

Heavy fermions compounds are generally discussed in 
the framework of a periodic Anderson Hamiltonian which 
reads (Stewart, 2001; Sun and Kotliar, 2005) 

'Hhf = '^E^ N n^i,N + ^ i,N i,N' 

i,N i,N,N',N'ytN 

+ C^fe.JV %N + ^(^f\,N Ci,N + /l.c) (5) 

k,N 

The first two terms represent the limiting form of Ti.Hubb 
Eq.(l) for f-electrons states with creation operators 
f\ jy, since the band is so narrow that it can be approxi- 
mated by a simple atomic level. The last term models the 
hybridization of the f state on site i with a wide band of 
electrons with creation operators c']; pf- Their Coulomb 
interactions can be neglected and the kinetic part corre- 
sponds to the third term on the right hand side of Eq.(5). 
Note that the summation over N runs over an a priori 
degenerate spin-orbit multiplet (although crystal field ef- 
fects and Hund's rule might of course reduce the number 
of terms). 

An important extension of the previous model concerns 

cupratc materials. We will dwell upon this case in a 
little more detail, since these compounds will serve as 



our preferred example for our discussion of impurities in 
strongly correlated material. It is generally accepted that 
modelling the electronic states in the copper oxide planes 
allows one to capture the most remarkable and relevant 
properties of these high temperature superconductors. In 
this case, a narrow Cu d band hybridizes with a wide 
O p band. In the undoped compound, the (electron) 
dj;2_y2 orbital is half filled while the p orbitals are fully 
occupied and we denote their energies by Ed and Ep, 
respectively. Since the difference E^ — Ep is small, the 
dx2_y2, Pa; and Py orbitals hybridize strongly (Fig. 4) and 
form the basis of the 3-band model. The p^ and d^z^^,.^ 
wavefunctions couple to the buffer (or reservoir) out of 
plane oxide layers and are usually not included in a model 
of the Cu02 planes. However, they have been shown 
to significantly renormalize the effective parameters in 
a 1-band model (Pavarini et ai, 2001). Conventionally, 
one introduces hole creation operators d'^i^a- and P^ j^a for 
the Cu and O orbitals respectively. The most relevant 
site energies, hopping amplitudes and Coulomb U arc 
shown in Fig. 4; their values have been obtained through 
computer studies (Eskes and Sawatzky, 1991; Hybertsen 
et ai, 1989; Kampf, 1994; McMahan et ai, 1990). Using 

ed = -Ed -Ud- 2Upd (6) 
Cp = —Ep — Up — 2Upd 

the model Hamiltonian reads: 

Hsband = ^d^rii^cr'' + ip^rij^a^ (7) 
i,a j,(7 

<hj>,c 

<j,i'>,o- i 

+UpY,nj,fnj,lP + Upd E ''^^^ 

j <i,j> 




FIG. 4 Cu and O orbitals in the copper oxide planes; relevant 
site, hopping and correlation energies are included. Integrat- 
ing out the O variables yields a t,t'-J model for the effective 
Cu lattice (right part of the picture). The notation is defined 
in the text. 
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The on-site and kinetic parts of 'H^band yield three 
bands, as shown in Fig. 4. In the hole basis, the lowest 
energy band (called antibonding in reference to electron 
states) is mainly of d character with a small admixture 
of p. It is half full, since one has one hole per Cu02 unit. 
The highest (bonding) band, mainly of p character, is 
empty and the charge transfer gap between the two is 
A = Cp — Cd- An empty (nonbonding) p band is centered 
aroimd Ep. Since Ud is very large, the splitting of the 
antibonding states induced by the Coulomb energy 
pushes the 3d^ — > 3d^ band (again in reference to elec- 
tron states) far above the bonding band. Since Up and 
Upd are significantly smaller than Ud, they do not affect 
the bonding and nonbonding bands too strongly. Thus, 
the "upper" (Sd^ ^ 3di°) Hubbard band is filled with 
holes and constitutes the ground state; it is separated 
from the nonbonding excited states by a charge transfer 
gap of order A (Fig. 5). (Zhang and Rice, 1988) pro- 
posed a further reduction of the 3-band model to an ef- 
fective single band model: a hole injected into the plane 
goes into an oxygen orbital since Ud is more than two 
times larger than A. In fact, it lives on the four oxygens 
surrounding a copper site, but its spin forms a so-called 
"Zhang-Rice singlet" with the copper spin. On that par- 
ticular plaquette, one has a charge but no spin. Only 
one singlet can occupy a given plaquette at a time. The 
Zhang-Rice singlet then hops from one plaquette to a 
neighboring one, carrying charge and restoring a spin-^ 
on the first plaquette. The low energy physics is then 
minimally described in terms oi & t — J model Eq.(2) 
on a lattice (each site corresponds to the original Cu02 
unit), where t is of order 0.4 eV and J of order 0.1 eV, 
as shown in Fig. 4 (Dagotto, 1994; Kampf, 1994). The 
Zhang-Rice states are shown in Fig. 5. Such a reduction 




FIG. 5 A schematic rendition of the main energy hole states 
for the 3-band and t-J models. The UHB (upper Hub- 
bard band), NB(nonbonding band), and LHB (lower Hub- 
bard band) pertain to the former while the S (Zhang-Rice 
singlet band) and T (Zhang-Rice triplet band) refer to the 
latter. The effective ciiargc transfer gap between the UHB 
and S bands is about 2eV. For the undoped material pictured 
here, the UHB is completely filled with holes. 



is not universally accepted, and there are many proposals 

which require three or more bands (Emery, 1987; Varma, 
2006). On the other hand, early NMR experiments sug- 



gested that magnetic fluctuations on Cu and on O sites 
track each other, upon varying the temperature or the 
doping, lending support to the Zhang-Rice scenario (Al- 
loul et ai, 1989; Takigawa et al., 1991). Still, ARPES 
experiments which are used to assign values to the tij 
reveal that, in general, one needs to incdude at least sec- 
ond neighbor Cu-Cu hoppings (Damascelli et al., 2003; 
Schabel et al, 1998). Neutron and Raman scattering al- 
low one to determine J and show that nearest-neighbor 
couplings dominate. So a popular model used to discuss 
the strongly correlated limit is the t — t' — J Hamiltonian 

<i,j>,<T 

-t' ^ {c^i,cr Ci'^a+h.cj + J ^ ('S'i ' Sj - ^Tlj Uj^ 
<i,i'>,a <i,j> 

where i' is a second neighbor of i and where the c^i^^ 
are hard core operators indicating that any site can be 
occupied either by a charge or by a spin, but not both. 
The sign of t' controls the topology of the Fermi sur- 
face (Damascelli et ai, 2003; Martins et ai, 2001) and 
when t' < one obtains a shape qualitatively similar 
to that seen in ARPES for the hole doped cuprates. A 
negative value also lowers the probability for phase sepa- 
ration (between hole rich and hole poor regions) to occur. 
While exact solutions of the homogeneous strongly cor- 
related models introduced above are not available in 2D, 
most practitioners in the correlated systems field agree 
on the qualitative aspects of their intermediate to high- 
energy spectroscopic features. In the crucial low-energy 
regions, attention is currently focused on questions such 
as whether or not they support superconductivity with- 
out the addition of couplings to the third dimension, or 
other ingredients. Enough understanding has been de- 
veloped in recent years to embolden some authors to 
study inhomogeneous states within such models, using 
the approximate techniques cited above. The great need 
for such studies is driven by fascinating phenomena dis- 
played by defects in ID systems, as well as by a com- 
pelling series of experiments on 2D cuprates, which we 
discuss below. 



III. IMPURITIES IN SPIN CHAINS AND LADDERS 

A. Overview 

One of the biggest challenges that condensed matter 
physicists are still facing today is the lack of a clear un- 
derstanding of the low energy properties of correlated two 
dimensional electronic systems. A route to the study of 
correlations uses the spatial dimension as a parameter. 
In ID, generic interacting electronic systems belong to 
the strong coupling limit (Luttinger or Luther-Emery) 
(Luther, 1976; Luther and Emery, 1974; Luttinger, 1963). 
In 3D, weak coupling physics appears to describe elec- 
tronic systems well. Two dimensional systems stand at a 
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dimensional crossroads between the strongly and weakly 
correlated limits. Since the scope of this review is fo- 
cused on strongly correlated systems, we will work our 
way from ID to 2D: we recall key features of ID physics 
and, adding chains in parallel, forming ladders, we grad- 
ually construct planar systems. This approach will help 
us gain insight into the physics of the 2D case, providing 
properties evolve smoothly as a function of the interchain 
couplings. 



B. The ID world 

One of the great advantages of considering the ID case 
is that many exact or quasi-exact results can be obtained 
in that limit, using different techniques. Bethe-Ansatz 
or transfer matrix allow us to solve such seminal models 
as the XXZ chain or the Hubbard model, revealing the 
nature of the ground state and of the elementary excita- 
tions. They allow us to compute thermodynamic proper- 
tics, and in some instances correlations functions (Cardy, 
1990; Lieb and Wu, 1968; den Nijs, 1988). Bosoniza- 
tion or mapping to non-linear sigma models are meth- 
ods which give access to the low energy properties of 
ID quantum problems, when the energy spectrum is lin- 
earized in the vicinity of the Fermi points (Cabra and 
Pujol, 2004; Giamarchi, 2004); combined with conformal 
theory, these two methods allow one to identify univer- 
sality classes, to compute critical indices and correlations 
functions accurately, and to give expressions for the cor- 
rections to scaling due to boundary conditions or to ther- 
mal effects. Functional renormalization group gives ac- 
cess to multiple energy scales and is thus able to capture 
important features when impurities are present (Ander- 
gassen et al., 2004; Meden et al, 2002). In conjunction 
with the above analytical methods, a number of pow- 
erful numerical techniques have been developed: Den- 
sity Matrix Renormalization Group (DMRG) or exact 
diagonalization provide information on the ground state 
properties and on correlations of ordered or quasi-ordered 
phases. Quantum Monte-Carlo (QMC) yields the tem- 
perature dependence of various observables and correla- 
tion functions, but, for fcrmionic variables, the infamous 
"sign problem" prevents one from simulating tempera- 
tures much less than the bare kinetic energy (David- 
son, 1993; Grotendorst et al, 2002; Schollwock, 2005). 
/,From these studies, one may identify two generic types 
of systems: those which fall into the Luttinger univer- 
sality class and those which belong to the Luther-Emery 
class. The former case corresponds to a critical regime 
where correlation functions decay as power laws of the 
distance; the system is compressible and excitations are 
gapless. The latter case pertains to a gapped regime; 
some or all of the excitations are gapped (incompressible 
state) and the corresponding correlation functions decay 
exponentially with distance. Many of the concepts that 
apply to ID systems and that one may hope to extend to 
higher dimensions can be expressed fairly simply in the 



language of bosonization. We restrict our presentation to 
the material that is necessary to that end and we further 
refer the reader to some excellent reviews on one dimen- 
sional problems (Affleck and Marston, 1988; Giamarchi, 
2004; Gogolin et al, 1999; Lecheminant, 2004; Mattis, 
1993; Schollwock et al, 2004; Tsvelik, 1995; Voit, 1995). 
In the case of a single band fermionic problem, low en- 
ergy modes correspond to excitations close to the Fermi 
surface which consists of two points. One then defines 
(real space) anihilation operators with spin a {a =t,i), 
tpr,a{x), for right moving fermions whose momentum is 
close to+fci?(r = -|-l) and for left moving fermions whose 
momentum is close to —kp {r = —1). 
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The bosonic phase operators (^aip^^ and n(a:) = ^VQ^{x) 
are canonically conjugated. The Klein factors (Majorana 
fermions) ri^.a ensure proper anticommutation of the V^'s 
with different a or r. In Eq.(9) only the lowest harmon- 
ics {-\-kF and —kp) are included. To study commensu- 
rability effects or the influence of a non linear energy- 
momentum dispersion one would include higher harmon- 
ics. At this stage let us highlight some similarities and 
differences between ID and 2D: 

• Interaction terms in the Hamiltonian mix low energy 
fluctuations associated with topologically distinct por- 
tions of the Fermi surface {^kp and —kp)- These back- 
ward scattering contributions are often relevant, i.e affect 
significantly the thermodynamic and transport proper- 
ties of the system. 

• Introducing the particle density operator 'Pa[x) = 
{kp — W 4>a{x)) / 17 , one may rewrite Eq.(9) as 



V2 



i[r J dy Pa(y)-Oa{x)] 



(10) 



This form is the ID analog of the expression for fermion 
operators used in the context of the 2D fractional quan- 
tum Hall effect (Fradkin, 1991; Pham et al., 2000). The 
"Jordan- Wigner" term 



i[r J dy p„(y)\ 



(11) 



is equivalent to the topological Chern-Simons term. This 
analogy is made more precise in (Pham et al., 2000) which 
shows that the eigenfunctions of the Luttinger liquid are 
ID Laughlin wavefunctions (Seidel and Lee, 2006). 

• If one defines charge (c) and spin (s) phase operators 
for the (j) fields 

'Ac(a;) = ^['^T(a=) + <l>l{x)l Mx) = -^i'f'fi^) - ^i(^)] 

(12) 

and similarly for the 6 fields, one finds that the bosonized 
form of the low-energy Hamiltonian of the interacting 
fcrmionic system is simply the sum of two terms, one 
describing charge fluctuations 
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+7K^ I dx COS {s/S(l)c{x) - Akp x) (13) 
and one describing spin fluctuations 

+7T^ / cosx/8<^«(a;). (14) 
(27ra)^ 7o 

Since the charge and spin phase fields commute, we see 
that wc have spin-charge separation. Furthermore, one 
may rewrite Eq.(9) as 

i>rAx) = -^Vr,ae""''' (15) 

V27ra 

The fermion operator Eq.(15) is the product of a charge 
field times a spin field. In the Luttinger regime, if the spin 
part has SU(2) symmetry, one may identify the spin field 
with the spinon operator, such that applying it to the 
ground state wavefunction creates the elementary quan- 
tum of spin excitation. By contrast, the charge part is 
not a simple holon operator (applying the holon opera- 
tor to the ground state wavefunction would create the 
elementary quantum of charge excitation) . 

This feature is reminiscent of the slave-boson (co- (x) = 
&^(x)/(j(x)) or slave-fermion {ca{x) = f''' {x)ba{x)) de- 
composition of the electron in 2D, but we have to keep in 
mind the fact that the pseudo boson (b) and the pseudo 
fermion (f) operators are not independent (they satisfy 
local constraints) so that these representations do not 
necessarily imply spin-charge separation. Let us em- 
phasize that spin-charge separation in ID is mostly an 
asymptotic property. For instance, the Bethe ansatz so- 
lution of the Hubbard model shows that it only holds 
in the low energy limit, except when U goes to infin- 
ity (the t-J model) when it applies at all energy scales 
(Ogata and Shiba, 1990). In the low energy sector, we 
can discuss magnetic and (charge) transport properties 
separately and in the folowing we concentrate on the for- 
mer aspect. 

C. ID physics in the spin sector 

1. Spin-j systems. In the presence of a commensurate 
potential that pins one charge per lattice site, the spin 
Hamiltonian Eq.(14) can be mapped onto the low energy 
representation of the XXZ spin- 5 chain, 

Hxxz = Jj2^i^i+i + + ^^i^Ui (16) 



This is done using a gauge transformation that changes 
5^ and S"^ to —S^ and —S^ every other site. Since 2kF = 
TT, the gi term (see Eq. (14)) corresponds to an Umklapp 
term = —JA as in Eq. (13) (we set the lattice spacing 
to unity). Keeping the 2kF harmonics, the bosonized 
expression of the S^- (x) operator at point x reads (see for 
instance (Giamarchi, 2004)) 

S^{x) = --V(j>+^^^^cos2(j, (17) 

The first piece is the uniform {q ~ 0) part while 

the second piece is the staggered {q ^ tt) part. The 
isotropic limit corresponds to A = 1 where one recovers 
the Heisenberg model. For translationally and rotation- 
ally invariant spin- i systems with sufficiently short range 
antiferromagnetie interactions, the Marshall and Lieb- 
Schultz-Mattis (LSM) theorems state that the ground 
state has total = and is either non degenerate 
with gapless excitations of odd parity, or has degener- 
ate ground states and breaks parity (reflection about a 
site), with a gap to the flrst excited triplet state (Lieb 
et al., 1961). Clearly, the Heisenberg model falls into 
the former category and represents the paradigm of spin 
liquids for spin-i spins (Haldane and Affleck generalized 
the previous theorems to arbitrary spins, see below). El- 
ementary excitations of this spin liquid consist of pairs 
of (asymptotically free) spin-i spinous (Bethe, 1931; des 
Cloizeaux and Pearson, 1962). Spin rotational symme- 
try implies Kg = I; the gi term then corresponds to a 
marginal operator which gives logarithmic corrections to 
the staggered part of the correlation functions (Koster- 
litz, 1974), 

{S'{x,0)S'{0,0))=Ci\+C2{-ir(-)\ogx^/\{18) 
X'^ \x/ 

These correlations can also be calculated at finite tem- 
perature, and they allow one to compute the uniform sus- 
ceptibility per atomic site, x (Bonner and Fisher, 1964; 
Eggert et al., 1994). At low T, setting fcs = 1 and in 
units of {gjiBY /J one finds that 

x(r)~-l(i + — i^). (19) 

TT^ 2 log ^-hjF- 

The correlation length is given by (Greven et al., 1996; 
Nomura and Yamada, 1991) 

^-?( 2-|,o„U.)l- )- 

Very few experimental realizations of a truly ID 
Heisenberg S = 1/2 chain have been achieved up to 
now. The best known example is the Cu chain Sr2Cu03 
(Ami et al.. 1995; Keren et al., 1993; Motoyama et al., 
1996). Magnetic susceptibility measurements by magne- 
tization techniques and by NMR show good agreement 
with Eq.(19) at low temperatures and with a Bonner- 
Fisher behavior over the whole temperature range (Mo- 
toyama et al, 1996; Thurber et al, 2001). The T depen- 
dence allows one to extract an antiferromagnetie coupling 
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between adjacent Cu sites of J = 2200 ± 200K. Long 
range order is observed only below Tpf = 5K, confirming 
the strongly one dimensional character of the chain. As 
will be shown in the following section, impurity effects 
enable one to determine correlation lengths compatible 
with Eq.(20). 

The rotational symmetry of the XXZ spin-i chain 
(Eq.(16)) is broken when A 7^ 1 in the expression for 
Hxxz- For A > 1 a long range antiferromagnetic state 
is obtained in the z direction. For, |A| < 1, the XXZ 
chain is in the Luttinger liquid regime. Bethe ansatz 
methods allow one to obtain exact expressions for the 
spin velocity 



(21) 



and the compressibility = 7r/arccos (— A) for all val- 
ues of A (Haldane, 1980; Luther and Peschel, 1975). 
Equal time spin-spin correlation functions decay as power 
laws at large distances x, implying a critical behavior and 
infinite correlation length. 

2. Spin-1 systems. The situation here is qualitatively 
different from that of the previous section. Based on the 
study of the low energy physics of the spin 1 Heisenberg 
chain, Haldane (Haldane, 1983) showed that the groimd 
state of rotationally and translationally invariant integer 
spin systems is a non degenerate singlet and that exci- 
tations to the first odd-parity state (triplet) are gapped. 
He established these properties using a non linear sigma 
model which showed that the difference with the half- 
integer case stems from the presence of a Berry phase in 
the latter case (Affleck, 1989; Haldane, 1983). This Berry 
phase is the generalization of the ID Chern-Simons term 
(see Eq.(ll) that we wrote down in the introduction to 
Section(III.B). Correlation functions are given by 



< s{x,o) ■ s{o,o) >~ {-ly 

Here ^ = Vs/AhuI where Anai 
0, 



,exp-|a;|/C 



(22) 



0.4105 J and, as T 
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(Affleck, 1989; Golinelh et al, 1994). 
Accordingly, the susceptibility varies as 



xiT) 



1 Affal 

exp = — 



Vt 



(23) 



(24) 



(Damle and Sachdev, 1998; Kim et al., 1998a; Tsvelik, 
1987). 

Experimental realizations of the 5 = 1 an- 
tiferromagnetic chain could be obtained in the 
Ni chains [Ni(C2H8N2)2(N02)]CI04 (NENP) and 
YBa2Ni05 (Buttrey et al, 1990). In both cases, a gap 
could be observed, confirming Haldane's conjecture (Bat- 
logg et al, 1994; Darriet and Regnault, 1993; Renard 



et al, 1988, 1987; Sakaguchi et al, 1996; Xu et al, 1996). 
In YBa2Ni05, the large intrachain coupling J = 280-?^, 
the small single ion anisotropy and the small ratio be- 
tween interchain and intrachain coupling {J' / J < 5 
10^^) gives this compound a very good ID Heisenberg 
chain character, with a Haldane gap A = lOOK (Bat- 
logg et al., 1994; Darriet and Regnault, 1993; Das et al., 
2004). In both componds, the temperature dependence 
of the susceptibility follows that given by Eq. (24) . 

Haldane's conjecture received convincing theoretical 
confirmation via the mapping of the spin Hamiltonian 



TT T \ ^ nx nx 



qV qV 



+ ASfS, 



D 



(5-) (25) 



onto a classical 2D restricted solid on solid model (den 

Nijs and Rommelsc, 1989). These authors showed that 
the massive phase can be understood in terms of the or- 
dering of a Jordan- Wigner like order parameter. The 
nature of the ground state and of the excitations were re- 
vealed thanks to the bosonization analysis of this model 
(Schulz, 1986). Each spin is written as the sum of two 
spin-1/2, Si = (Ti + Ti (Fig. 6). Insofar as singlet states 
on any given site are not relevant in the low energy sec- 
tor, Hgi is now the sum of two coupled spin-| systems. 
The most relevant coupling is of the form 



g J dxcos {26o{x)) 



(26) 



where 9q is an "optical mode"- like phase field, related 
to an antisymmetric combination of 9 operators (see 
Eq.(9)). Since the renormalized g flows to minus infinity, 
^0 has to be either or tt in order to minimize the energy 
in the ground state. In the context of bosonization, the 6 
field characterizes a superfluid phase. When it acquires 
a finite expectation value, a Josephson superfluid current 
flows across the system. Thus we dub the term Eq.(26) 
a Josephson coupling. In spin language this term comes 
from kinetic contributions of the form afT~^i. The long 
range ordering of this quantity is responsible for a spin 
gap proportional to g. It signals the formation of an RVB 
phase, i.e a non degenerate, rotationally and translation- 
ally invariant valence bond gapped phase. 



S=l 



FIG. 6 Artist's conception of a valence bond state. Each 

spin-1 S is written as the sum of two spin-i states, a + f; 
valence bonds are formed between a a and a t belonging to 
neighboring sites. 



This term is the counterpart of the spinon hopping 
amplitude order parameter introduced in the framework 
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of the 2D slave-boson representation of the t — J model 
(Kothar and Liu, 1988; Lee and Nagaosa, 1992). Exci- 
tations break the valence bond singlets promoting those 
to the triplet state. This picture emerges from the ex- 
act solution of the AKLT model (Affleck et at, 1987). It 
is the spin-1 Heisenberg model with an additional term 

^JJ2i(^^i ' ^i+ij ■ Valence bonds are formed between 
nearest neighbor sites and spin-spin correlation functions 
decay exponentially with distance, with a correlation 
length £,AKLT = l/log3 (Arovas and Auerbach, 1988; 
Arovas and Girvin, 1991). The inclusion of the addi- 
tional term does not affect the low energy physics of the 
Heisenberg model. These particular valence bond states 
are sometimes considered to be characteristic of a bosonic 
RVB phase (Arovas and Auerbach, 1988; Ng, 1992) and 
to be relevant to the physics of cuprate or other strongly 
correlated materials. It is noteworthy that in the a-r 
language, the Heisenberg Hamiltonian can be viewed as 
a two-leg ladder model, with intrachain and interchain 
couphngs (Fig. 7). 




FIG. 7 Left: a two-leg spin-1/2 ladder. Dots represent the 
spins with couplings J along the legs and Jrung along the 
rungs. Right: a spin-1 chain can be viewed as a two-leg spin- 
1/2 ladder with diagonal couplings (dotted lines) 



3. Higher spin states. Affleck presented a physicist's 
proof extending the ground state and excitations proper- 
ties of the spin-i case to general half-integer spins (Af- 
fleck, 1989). Haldane also predicted that general integer 
spins fall into the same class as the spin-1 case. These 
characteristics are difficult to prove rigorously, for want 
of exact solutions for general spin S, but one can use the 
decomposition of 5' into 25* half-integer spins for integer 
S and 2(5 — 1/2) -I- 1 half-integer spins for half integer S 
(Schulz, 1986). Forming valence bonds between different 
"cr" and "t" species on neighboring sites, one is left with 
a "single" spin-i chain for the half-integer case. This 
gives an overall gapless regime, with massless spinons. 
For the integer spin case, on can generalize the AKLT 
model (Oshikawa, 1992) for any S, and this gives an ex- 
act RVB (valence bond) gapped state, with a correlation 
length £,AKLT = 1/ log (1 -I- 2/5). Again we can view 
this latter case as a model for an S'-leg ladder system in 



the low energy limit {S chains with interchain couplings 
along rungs). 

4-. Lord of the rungs. On the road to 2D, we now turn 
to spin-i ladders. Two and three leg Cu ladder structures 
were experimentally realized in SrCu203 and Sr2Cu305 
oxides (Azuma et al., 1994; Hiroi et al., 1991). Couplings 
along the rung and along the leg are similar, of the order 
of 1000 to 2000K. In the two-leg compound, the ladders 
are almost decoupled from each others since the system 
does not show any magnetic ordering down to T=50 mK 
(Kojima et ai, 1995a, b). This is not the case in the three- 
leg compound, which orders at 52 K, probably due to 
weak transverse couplings. A gap is observed only for the 
two-leg system, as shown by macroscopic measurements 
and NMR (Azuma et at, 1994; Ishida et ai, 1994, 1996). 
In a similar ladder cuprate Sri4_2:Caa;Cu2404i, hole dop- 
ing could be achieved through Ca-Sr substitution, leading 
for high enough pressure to a superconducting behavior 
at X = 13.6 (Uehara et al., 1996). 

In these spin ladders, exact solutions arc not avail- 
able, but, for low energies, we can use some of the 
techniques that we discussed above to determine ground 
states and excitation properties (Citro and Orignac, 
2002; Defl'Arringa et ai, 1997; Gogolin et al, 1999; 
Nersesyan et ai, 1998; Orignac and Giamarchi, 1998). 
In addition, DMRG, exact diagonalization and Quan- 
tum Monte-Carlo (no sign problem for spins!) also add 
useful information (Dagotto, 1999; Greven et ai, 1996; 
Johnston et ai, 2000; Poilblanc et ai, 1994c; White and 
Noack, 1994). In the case of the 2- leg ladder, bosoniza- 
tion is easy to perform after introducing acoustic (sym- 
metric under the exchange of the two legs) and optical 
(antisymmetric) modes (Strong and Millis, 1994). Re- 
markably, the form is precisely that obtained by Schulz 
in his analysis of the Haldane spin-1 chain in terms of 
two coupled spin-i legs (see subsection 2.). Accordingly, 
the spin sector is fully gapped as soon as any interchain 
coupling is turned on. Once again (see Eq.(26) and subse- 
quent comments), a Josephson coupling between the two 
chains of the form (jfT~^^ is relevant and provides the 
gap between the valence bond RVB ground state and the 
first triplet excited state, see Fig. 7. When comparing 
the spin-1 chain with a two-leg spin-i ladder, one should 
bear in mind the particular topology (see right panel of 
Fig. 7); it plays an important role, for instance in the 
presence of defects or of open boundary conditions (Ng, 
1994; Qin et ai, 1995). 

Similarly, the analysis of the 3-leg ladder within 
bosonization leads to the result that one sector is gap- 
less. As previously mentioned, a gap was experimen- 
tally observed in the two- leg compound SrCu203 while 
no gap was detected in the three-leg material Sr2Cu305. 
It is possible to analyse the ladder problem with arbitrary 
number of legs (nj), using the non linear sigma model, 
and one finds indeed that an odd number of legs leads to 
a gapless state with power law decay of the spin-spin cor- 
relation functions, while for an even number of legs corre- 
lations fall off exponentially, such that ^lad ^ exp (§«;) 
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(Dell'Arringa et aL, 1997). In the weak coupling limit 
{J rung << >^)) the gap decreases exponentially to zero 
with increasing number of legs: Ajad ^ ■ These re- 
sults are confirmed, and extended to the isotropic limit 
{Jrung = J) using numerical work (Greven et al, 1996): 
Aiad = 0.50,0.16,0.05 for m = 2,4,6. At large dis- 
tances, spin-spin correlation functions behave exactly as 
in the spin-1 case (Eq.(22)) and the susceptibility is given 
by Eq.(24) where one substitutes Aiad for AhuI- Here 
£,'= 3.2,10 for ni = 2,4. 



5. Summary of ID physics in the spin sector. A pic- 
ture emerges from the above description: for (spin) rota- 
tionally and translationally invariant systems, the ground 
state of ID spins chains and ladders is either a degen- 
erate dimer phase where reflection symmetry about a 
site is broken, or a non degenerate valence bond phase, 
a.k.a a Resonating Valence Bond phase (RVB), where 
parity symmetry exists for any site (Anderson, 1997). 
It consists of singlets, which resonate coherently on a 
length scale ^. Since the spectrum is linearized near the 
Fermi points, the binding energy in this RVB state is 
A ^ Vs£,^^. Within the RVB phase, quasiparticles con- 
sist of spinons. These are thus confined within a distance 
^ of each other, forming bound state singlets. An exci- 
tation breaks the singlet, leading to a triplet state. This 
triplet can propagate through the lattice as a magnon 
excitation. For strictly ID chains at T = 0, ^ is infinite 
in the spin-i case, and we have a gapless spin liquid with 
asymptoptically free spinons while ^ is finite in the spin- 
1 case. For ladders, we have a gapped RVB phase and 
Clad ^ exp (f n;) for even numbers (ni) of legs. The limit 
n; ^ oo is singular, since the ground state of a 2D spin-i 
system is believed to sustain long range antifcrromagnetic 
order (Chakravarty et al., 1989). The RVB phase that 
constitutes the ground state of spin-1 and even- leg ladder 
systems is a strongly correlated phase in that it has long 
range phase coherence. This property is quite indepen- 
dent of the strength of the interactions. Can we actually 
see spinons? Simulations are of course one way to do so: 
a ID ring with an odd number of spins shows a spinon 
extending over the entire circumference. Open chains are 
another way of revealing these quasiparticles. Prom an 
experimental perspective, there is a way to realize these 
situations: adding nonmagnetic impurities breaks the ID 
chain. We will discuss this point in the next section. An- 
other way to probe spinons - akin to the odd number 
of spins case - is to consider ladders with odd number 
of legs. In that case, the "unpaired" spinons form a ID 
spin liquid phase, which is revealed in experiments. If we 
associate a Josephson phase (9) with each singlet in the 
RVB phase, this 6 field is either quasi-ordered in the cases 
of ID spin liquid and odd-leg ladders or ordered for even- 
leg ladders: in the former situation, correlation functions 
of the 6 field decay algebraically with distance, while in 
the latter they tend to a constant value for oo. 



D. Defects in ID spin chains 

^From an experimental perspective, intrinsic defects 

will modify the magnetic properties of materials. Exam- 
ples of impurities arc weak- (magnetic coupling) links, 
broken chains, magnetic or nonmagnetic atomic substi- 
tutions (Brunei et al, 1999; Eggert and Affleck, 1995; 
Eggcrt et aL, 1994; Igarashi et aL, 1995; Rommer and 
Eggert, 1999, 2000: Zhang et al, 1997, 1998). They pro- 
duce spectacular effects which are direct consequences of 
the correlated ID states so that we may actually use de- 
fects in a controlled way in order to probe the ground 
state and the excitations in the spin chains that we in- 
troduced in the previous section. Nonmagnetic atoms su- 
prcss the couplings J between neighboring spins. Thus, 
in ID, a finite concentration of such atoms will break the 
chain into disconnected segments. At low temperature, 
each segment is in its ground state. In the case of spin-S' 
Heisenberg chains, Marshall's theorem tells us that the 
total spin is zero for an even number of spins in a chain 
or S for an odd number of spins in the chain. Thus 
odd chains will give a Curie contribution to the uniform 
magnetic susceptibility proportional to the concentration 
of impurities (assuming small concentrations and equal 
probability of getting even or odd chains) (Sirker et aL, 
2006; Wessel and Haas, 2000). This property does not tell 
much about the nature of the correlated state per se. In 
contrast, we show below that boundaries c;ause additional 
phenomena which reveal the nature of the low tempera- 
ture phases; the reason is because the spin density in the 
z direction is the sum of two terms (see Eq.(17)): a uni- 
form part (oc V(/>) which is essentially unaffected by the 
breaking of translational invariance, except in the vicin- 
ity of the edges plus a staggered (Umklapp) part giving 
two types of contributions, one stemming from finite size 
effects and one from the boundaries (Eggert et aL, 2002). 
We first consider spin-^ compounds, then discuss spin-1 
chains and ladders and conclude with quasi- ID systems 
(Eggert and Affleck, 2004). 



1. Spin-i systems 

Let us first discuss the influence of boundaries on 
the magnetic properties of a spin-| isotropic Heisenberg 
chain. Bethe ansatz, conformal theory or bosonization 
techniques allow one to compute all relevant observables 
(Cardy and Lewellen, 1991; Eggert and Affleck, 1995; 
Essler, 1996; Fujimoto, 2000; Fujimoto and Eggert, 2004; 
Lukyanov, 1998). In bosonization language, a bound- 
ary forces a condition on the r = — 1 and r = +1 
modes (Eq.(9)) (see for instance (Giamarchi, 2004), p 
306). Breaking translational invariance affects spin-spin 
correlation functions. Far from the boundary, one re- 
covers the usual behavior; for instance, if one applies 
a uniform magnetic field H, one gets a magnetization 
equal to x{T)H , where the uniform susceptibility is given 
by Eq.(19). Close to the boundary, the susceptibility 
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becomes site dependent; it is the sum of two contribu- 
tions. One of those has a magnitude close to x{T)H, 
such that in Fourier space its components are centered 
around g = (see (Eggert and Affleck, 1995)). In ad- 
dition, because of the boundary, the alternating part of 
the spin-spin correlation function also contributes to the 
uniform susceptibility. This is an impurity induced ef- 
fect. For finite temperature, one gets a site dependent, 
staggered magnetization of the form (Eggert and AfHeck, 
1995) 



M{x) = H{-iy 



,1.16 



^sinh(^) 



(27) 



At low enough temperature T <C J, the envelope of 
the magnetization increases initially with the distance to 
the impurity site, peaks at x = 0A8J/T in cell units, 
then decreases with an exponential asymptote of typical 
extension ^j^p as shown in Fig. 8. From Eq.(27) one 
finds that ^imp ^ T~^. 
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FIG. 8 Alternating magnetization caused by a boundary in 
a spin-i chain according to Monte Carlo simulations at T = 
J/ 15 (Eggert and Affleck, 1995). 



One needs to be a bit careful when taking the T = 
limit. The expression for the magnetization given by 
Eq.(27) with T set to zero is restricted to "small enough" 
values of x. The proper form reads 



M(x)~(-l)y^sin(|^) 



(28) 



and shows that a finite field H produces an incommensu- 
rate oscillating magnetization (on top of the (quasi-) uni- 
form part). On sites close to the boundary, the effect of 
the "impurity" is seen through extra contributions to the 
magnetization which are nonzero even at T = (Rom- 
mer and Eggert, 1999). Recently, (Fujimoto and Eggert, 
2004) showed that one needs to include a term gener- 
ated by marginally irrelevant bulk operators. It gives a 
quasi-Curie contribution to the susceptibility of the form 

Xb ~ (^log {T)^ . Note that this additional piece ex- 
ists both for even and odd chains, since this is a "surface" 
effect. 



In practice, creating a boundary does not require one 
to cut the chain. Indeed, a weak link (a smaller exchange 
coupling between two neighboring spins) is a strong per- 
turbation and ultimately, as the temperature is reduced, 
renormalization and numerical techniques show that it 
produces similar effects on the susceptibility to the case 
of the open chain (Rommer and Eggert, 1999). 

As we stated in the introduction to this section, defects 
break a chain into disconnected segments. The bound- 
aries of each segment give rise to extra "surface" con- 
tributions to the magnetization. In addition, there are 
"volume" contributions associated with finite size effects: 
since the staggered coherence length diverges at T = 
in the thermodynamic limit, one gets a low T finite al- 
ternating moment per site oc for a chain of length 
L (Eggert et al, 2002). This size dependence is a direct 
consequence of the fall-off of the staggered correla- 
tion function. For odd-L chains this magnetization is a 
manifestation of the presence of a spin-1/2 in the groimd 
state, and for even-L chains it accompanies a spin-1 ex- 
citation above the ground state (see Fig. 9). The local 
magnetic susceptibility Xi on site i is the sum of a uniform 
contribution and an alternating contribution. Fig. 9 
shows x" and the staggered part of the alternating term, 
namely Xi = (-l)'(Xi - X")- 
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FIG. 9 Uniform (xT) and staggered (x?) components of 
the envelope of the local susceptibility Xi obtained with the 
DMRG method applied to the S = 1/2 Heisenborg model on 
a chain with L sites and open boundaries, using m states in 
the iterations and studying the subspace with a total spin in 
the z direction equal to 1. (a) corresponds to L = 80 and 
m — 16, and (b) to L = 33 and m = 32 (Laukamp et al., 
1998). 

On the experimental side, the Sr2Cu03 spin chain in- 
troduced in section (III.C.l) is a good realization of a 
truly ID Heisenberg 5=1/2 system, with an antiferro- 
magnetic couphng J = 2200 ± 200K (Ami et al, 1995; 
Keren et al., 1993; Motoyama et al., 1996). Only non- 
magnetic Pd substitution could be achieved on Cu site, 
which leads to a Curie behavior of the macroscopic mag- 
netization (Kojima et al., 2004). This Curie term is a 
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consequence of the chain segments with odd number of 
chains, as argued above. Local probes such as NMR are 
necessary to probe the local variation of the magnetiza- 
tion M{r). Measurements were performed only in the 
case of a non-substituted chain, taking advantage of the 
presence of unavoidable intrinsic defects, probably oxy- 
gen vacancies and interstitials (Boucher and Takigawa, 
2000; Takigawa et at, 1997). The NMR spectrum shows 
a central line and a broad symmetric background with 
sharp edges (see Fig. 10). This is evidence of the pres- 
ence of an alternating polarization induced by the na- 
tive defects, as will be demonstrated in Fig. 14. The 
edges of the AH pattern follow a l/VT shift, as would 
be expected from Eq.(27) for values of x near the peak 
position. 
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FIG. 10 left panel : the NMR spectrum expected for a stag- 
gered magnetization of the spin 1/2 chain displayed in the 
inset, which corresponds to its histogram. Right panel : the 
experimental Cu NMR spectrum in Cu chain Sr2Cu03 dis- 
plays a similar AH pattern. It is interpreted as induced by 
intrinsic defects (Takigawa et al., 1997). 



In the Cu spin 1/2 chain dichlorobis (pyridine) cop- 
per(II) (CUCI22NC5H5), referred to as CPC, with J = 
27-?^, muon spin resonance was performed, where im- 
planted muons were used both as a probe of the local 
susceptibility and as a defect in the chains (Chakhalian 
et al, 2003). Indeed, muons are positively charged par- 
ticules which stop in the bulk of the material and then 
couple to the chains, according to comparisons of Knight 
shift measurements with DMRG computations. While 
this technique is valuable in order to study chains where 
native defects cannot be produced by simple chemistry, 
it is limited by the fact that the muon location is not 
exactly known. 



2. spin-1 chains 

In many respects, it is easier to describe the response 
of a spin-1 system to an impurity. As we pointed out 
above, the ground state of the pure system is a non de- 
generate RVB state, with long range order (of a singlet 
order parameter). Introducing a nonmagnetic impurity 
breaks a valence bond and, in the representation of a 
spin-1 in terms of two spins- ^, implies that a spin-i is 



"freed" at the location of the broken valence bond, i.e it 
is pinned in the vicinity of the impurity, as shown in Fig. 
11 (Kennedy, 1990; Laukamp et ai, 1998; Miyashita and 
Yamamoto, 1993; Ng, 1992). 
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FIG. 11 The introduction of a non magnetic impurity in the 
Spin-1 chain displayed in Fig. 6 breaks two singlet bonds, 
thereby releasing two spins- i. 

Monte Carlo simulations performed on chains of L sites 
(L > 100) show that the polarization oscillates with an 

r (i-1) 1 

exponential envelope < >cx e s ' + e ^ « 
with ^ 6. A comparison of the induced magnetization 
for the S ~ ^ and S — 1 cases is shown in Fig. 12. 
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FIG. 12 Alternating magnetization in the case of spin-| and 
spin-1 systems with a defect. The polarization for S = 1 
is normalized to that for S — ^. At high temperature, the 
two polarizations are roughly identical, exponentially decay- 
ing. At low temperatures, for S = 1, the polarization doesn't 
evolve further because ^ saturates. On the other hand, the 
extension of the polarization for S = | systems continues 
to increase as T is decreasing and exhibits a nonmonotonic 
behavior (see Fig. 8). 

This value of ^ is that of the undoped compound. 
This result suggests that the perturbation which created 
the impurity is localized and decays with the coherence 
length of the pure system see Eq.(23). This fact is con- 
firmed in various DMRG and QMC studies (Alet and 
Sorensen, 2000; Kim et al., 1998a; Miyashita and Ya- 
mamoto, 1993; Polizzi et ai, 1998; S0rensen and Affleck, 
1995). A remarkable feature emerges from these simula- 
tions, namely, the relative sign of the polarization at the 
end of a segment of the chain formed by two impurities. 
Assume that we have a spin-up state at one end of the 
segment. The oscillation of the polarization implies that 
the next site is in the down state and so on. If we ex- 
tend the up-down alternation of the sites all the way to 
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the other end of the segment, we observe that the po- 
larization on the site closest to the second impurity cor- 
responds to that expected from the up-down alternation 
rule, despite the fact that the distance between the impu- 
rities is much larger than 2^. In other words, the singlet 
order parameter propagates the relative phase of the spin 
operators across the entire segment (see (Miyashita and 
Yamamoto, 1993)). 

This observation is confirmed by a study combining 
theory and experiments (Hagiwara et al, 1990) on a spin- 
1 chain with a spin-^ (Cu) impurity. The lowest energy 
state is such that sites located on either side of the Cu are 
BOTH in the up or down state and form a triplet with 
< S"^ ±1, as shown in Fig. 13. DMRG and QMC 

FIG. 13 Spin-1 chain with a Cu spin-1/2 defect. Spin-1/2 
states are formed on either side of the Cu atom and align to 
form a spin-1 state. 



confirm this picture (Roos and Miyashita, 1999; S0rensen 
and Affleck, 1995; Tonegawa and Karubaki, 1995). 

Contrary to the case of the S = ^ spin chain, the 
impurity induced magnetization in the 5* = 1 spin 
chain could be measured experimentally with great pre- 
cision. In the S = 1 antiferromagnetic Ni chain 
[Ni(C2H8N2)2(N02)]CI04 (NENP), EPR measurements 
provide evidence for the existence of a spin state asso- 
ciated with the substitution of Ni by nonmagnetic Zn, 
Cd or Hg atoms, consistent with the VBS picture where 
each impurity frees two S = ^ states (Glarum et at, 
1991). Specific heat measurements confirm the existence 
of a spin state by detecting the corresponding Schottky 
anomaly (Ramirez et ai, 1994). While this study was first 
argued to give a different spin state value than EPR, the 
controversy was finally resolved by a more careful treat- 
ment of long enough spin chains (Batista et ai, 1998). 
More refined experiments were then performed in the Ni 
chain oxide YBa2Ni05, a better prototype of the spin 
1 chain than NENP, with J = 280K (Buttrey et ai, 
1990). When introducing nonmagnetic Zn or Mg im- 
purities, EPR studies reveal the presence of spin states 
accompanied by a Curie contribution to the macroscopic 
susceptibility (Batista et ai, 1998; Das et ai, 2004; Payen 
et ai, 2000; Ramirez et ai, 1994). The Curie constant 
corresponds to about two S = ^ per Zn, as expected. 
Note, however, that any attempt to analyze the macro- 
scopic magnetization quantitatively is hindered by the 
fact that a Curie behavior is already observed in the 
nominally pure compound without any Zn, a feature en- 
countered in most experimental realizations of low di- 
mensional systems. This native Curie term is usually due 
both to the presence of spurious paramagnetic phases and 
local defects within the crystal itself, such as vacancies or 



interstitials of oxygens. As anticipated, any defect which 
is coupled to the spin chain itself will induce in turn a 
local moment, i.e. a Curie law in a magnetization mea- 
surement. This limits any refined quantitative analysis 
of such macroscopic experiments, as well as EPR studies. 
To bypass this difficulty, NMR has proven to be an ideal 
probe, since it allows one to identify the actual effect 
of the Zn substitutional atom in its vicinity and more- 
over to resolve this effect spatially, while other techniques 
only give the total spin value induced by the defect. In 
YBa2Ni05, NMR of the yttrium nucleus probes the local 
magnetization M(r) of the Ni sites, with a site by site 
resolution. Without impurity, the ^^Y NMR shift mea- 
sures the pure uniform susceptibility of the Ni chain and 
displays the usual Haldane gap behavior (Shimizu et at, 
1995). When adding nonmagnetic isovalent impurities 
such as Mg^+ or Zn^+, new NMR satellite lines appear, 
each of them being due to all the Y nuclei at a given 
distance Ri from an impurity as shown in Fig. 14 (Das 
et ai, 2004; Tedoldi et ai, 1999). 




chain site (cell units) 

FIG. 14 NMR reconstruction of an impurity-induced mag- 
netization in a spin chain. The left panel displays an 
NMR spectrum in pure (green) and Zn-substituted (red) 
YBa2Ni05. The various satellites which appear with Zn dop- 
ing allow the reconstruction site by site of the induced stag- 
gered magnetization near the Zn atoms substituting Ni in the 
chain (Das et al, 2004; Tedoldi et al, 1999). 



Their frequency shift relative to the pure case is di- 
rectly proportional to the local average spin {Szii}) along 
the external field direction z at the corresponding Ni site. 
Since the polarization (S^) alternates, the shifts relative 
to the main line produce both positive and negative satel- 
lites, as shown. They decay exponentially from the Ni 
first neighbor of the impurity site. Their temperature 
dependence is governed by that of the Ni closest to the 
impurity, which follows approximately a Curie law, to- 
gether with the T-dependence of the envelope extension 
Cimp plotted in Fig. 15. This spatial behavior is a con- 
sequence of the exponential decay of the spin spin cor- 
relation functions in such Haldane gapped spin chains. 
Furthermore, the extension S.imp measured for different 
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FIG. 15 Comparison of the extension of the staggered in- 
duced magnetism due to Mg (Tedoldi et ai, 1999), Zn or Cu 
(Das et ai, 2004) impurities in YBNiO. The correlation length 
measured by QMC by (Kim et al., 1998b) for the pure system 
is shown for comparison. 



impurities appears almost identical to that of the pure 
spin-spin correlation function S^pure computed by QMC 
techniques, also displayed in Fig. 15. The fact that the 
impurity reveals the intrinsic spin-spin correlations of the 
system is the key result which will give us insight into the 
other more complex and less controlled correlated sys- 
tems such as high Tc superconductors: a nonmagnetic 
impurity can be used as a probe of the intrinsic prop- 
erties of a correlated system, very much like a magnetic 
impurity in a metallic host could be used to probe the 
Fermi surface properties through the RKKY oscillation. 

Inelastic neutron scattering experiments also allow one 
to measure the end-chain excitations, which are peaked 
at the antiferromagnetic wavevector along the chain, and 
show an alternating pattern with extension ^imp = 8 at 
T=0.1K, consistent with the NMR studies (Kenzelmann 
et al, 2003). 

The case of a magnetic impurity in a spin-1 chain has 
been stressed previously: a Cu spin-i substitution on 
a spin-1 site results in the formation of a spin-i effec- 
tive state as probed by EPR in NENP (Hagiwara et al., 
1990) (see Fig. 13). Specific heat measurements con- 
firm this observation through the observation of a Schot- 
tky anomaly as in the case of nonmagnetic Zn (Batista 
et al, 1998). In addition, DMRG and QMC computa- 
tions predict that the magnetic impurity induces an al- 
ternating magnetization close to the defect (Roos and 
Miyashita, 1999; S0rensen and Affleck, 1995; Tonegawa 
and Karubaki, 1995). In the case of YBa2Ni05, ^^Y 
NMR experiments similar to those reported for Zn and 
Mg above show similar results for magnetic Cu impuri- 
ties (Das et al, 2004). The envelope remains exponential 
with the same extension S^imp as shown in Fig. (15). In 
comparison with nonmagnetic Zn defects, one sees a dif- 



ferent shift for these satellites, i.e. a smaller amplitude 
of the induced polarization, by about 20%, for all the Ni 
sites producing the polarization. A QMC computation 
can be performed in this specific case, which reproduces 
the experiment if Cu is coupled antiferromagnetically to 
its neighboring Ni with a coupling J' ~ 0.1 — 0.2J (Das 
et al, 2004), a result qualitatively similar to that found 
in NENP (Hagiwara et al, 1990). 



3. Ladder systems 



In Section III.C.4, we argued that the low energy 
physics of an even-(n) leg spin-^ ladder could be mapped 
onto an integer spin- 2. chain and that of an odd -leg lad- 
der to a spin-i chain. Thus, the influence of an impurity 
in the ladder follows from the considerations pertaining 
to the spin-1 and spin-^ cases. For instance, for a two 
leg ladder, we conclude that when a nonmagnetic impu- 
rity is located on one of the legs, a spin-i state appears 
on the opposite leg; its amplitude is centered around the 
site across from where the impurity sits. We emphasize 
once again (see discussion in Section III.C.4) that it is 
important to keep in mind the particular topology of the 
ladder that the above mapping involves, when discussing 
end-states. For instance, a spin-^ state is localized near 
the edge of an open spin-1 chain whereas no such state 
is found at the open end of the two-leg ladder shown of 
the left panel of Fig. 7. 

Numerical studies of the spin-^ ladder with impurities 
indeed find a magnetic polarization which oscillates from 
site to site along the leg of the ladder, with an amplitude 
which decays exponentially with a characteristic length 
fimp (Laeuchli et al, 2002; Laukamp et al, 1998; Martins 
et al, 1996, 1997; Mikeska et al, 1997; Miyazaki et al, 
1997; Sandvik et al, 1997). But whether or not the value 
of £,imp is equal to the coherence length of the pure system 
(equal to 3-4 lattice spacings) is still controversial. 

These polarized states either produce a Curie law for 
the susceptibility or may even lead to an antiferromag- 
netic transition (Greven and Birgeneau, 1998). For the 
3-leg ladder the physics seems consistent with that of 
the spin-i chain with a nonmagnetic impurity (Sandvik 
et al, 1997). Experimentally, nonmagnetic Zn can be 
substituted in the two-leg ladder SrCu203. Macroscopic 
susceptiblity measurements reveal the appearance of a 
Curie law at high temperature, which corresponds to one 
spin-^ per Zn in the dilute limit (Azuma et al, 1997). A 
concomitant 1 /T-broadening of the Cu NMR line is ob- 
served, proportional to the number of impurities, which 
provides evidence for the existence of an alternating in- 
duced magnetization which extends over the ladder (Fuji- 
wara et al, 1998; Ohsugi et al, 1999). No satellites could 
be resolved, in contrast with the case of YBa2Ni05. This 
limits any quantitative evaluation of the magnetization 
shape and extension. Magnetic S* = 1 Ni substitution 
leads to similar effects. 
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4. Summary of impurities in spin chains 

Impurities create edges in spin chains. For strictly ID 

systems, pairs of edges divide the system into even- or 
odd- length segments. The magnetic response to such 
defects contains two parts, which mirror the two con- 
tributions to the spin density operator of the pure sys- 
tem: (a) the uniform contribution - the gradient of the 
phase - is essentially identical to that of the impurity 
free system, except close to the boundaries and (b) the 
Umklapp contribution near the antiferromagnetic wavec- 
tor, by contrast, is qualitatively affected by the bound- 
ary; this is because the RVB parameter (for spin-1 chains 
and for ladders) or the staggered magnetization (for spin- 
^ chains) show long range or quasi-long range order at 
this wavevector for the pure system. Since impurities 
break translational invariance, momentum is no longer 
a good quantum number so that a uniform perturbation 
can give rise to a spatially oscillating magnetization and, 
conversely, an alternating perturbation can contribute to 
the uniform response. This is precisely what is seen in 
chain segments. The ground state of odd chains is a spin- 
i doublet. An infinitesimal magnetic field will give a uni- 
form Curie contribution which, in turn, will promote an 
alternating magnetization. For even chains, one also gets 
the same effect at finite temperature when transitions to 
the excited (triplet) state are allowed. The reverse effect 
is observed (at least numerically) in the case of a single 
edge, in the presence of a magnetic field. In this case, 
the induced uniform spin polarization causes an alter- 
nating magnetization to develop near the boundary. Its 
spatial extension is of the order of the coherence length 
of the pure system and this mode contributes to the uni- 
form magnetic response. These features are hallmarks of 
strongly correlated systems where several length scales 
are simultaneously relevant (Anfuso and Eggert, 2006). 
We will see in a forthcoming section that similar effects 
occur in two dimensions, when nonmagnetic defects are 
introduced in the copper-oxide planes. 



5. Impurities and long range order 



There cannot be true long range order at finite tem- 
perature for purely ID systems, but adding a small in- 
terchain coupling Jx between spin chains will drive the 
quasi-lD system antiferromagnetic. The Neel tempera- 
ture Tjv can be determined from the exact expression of 
the staggered susceptibility of the chain, using RPA (Eg- 
gert et al, 2002, 2003). One finds that Tn J±. Non 
magnetic impurities do cut the chains into segments, but 
LRO is not destroyed since two segments belonging to a 
given chain can still be magnetically connected through 
the interchain J± . In fact Tjv / Jl is a function of the 
scaling variable [(1/p) — l]Jj_. As the concentration p 
of impurities is increased, Tjv falls. However, in some 
situations, the magnitude of the staggered moment and 



the Neel temperature may increase with p (Eggert and 
Affleck, 2004). Similar effects can be obtained for a large 
enough single ion anisotropy (Sakai and Takahashi, 1990; 
Zhcludev et a/., 2000). 

In other words, the introduction of impurities enhances 
long range order. This is one of the various routes to get- 
ting "order from disorder" (Shender and Kivelson, 1991; 
Villain et al., 1980). For example, in two-leg spin-i lad- 
ders SrCu203, specific heat and susceptibility measure- 
ments show that a concentration x = 1% of Zn at Cu site 
leads to an antiferromagnetic behavior at T<T]v = 3K 
(Azuma et al, 1997). When increasing the number of im- 
purities, Tjv increases to 8K at 4% and decreases again. 
Concentrations as small as x = 0.1% still lead to an an- 
tiferromagnetic order (Manabe et ai, 1998). In spin- 
Peierls chains CuGcOa as well, either magnetic (Ni,Co) 
or nonmagnetic (Mg,Zn) defects at Cu site simultane- 
ously reduce the spin-Peierls transition temperature Tsp 
and induce a long range AF order whose frozen moment 
is spatially inhomogeneous (Anderson et at, 1997; Hase 
et al., 1993; Kadano, 1997; Kojima et ai, 1997; Oseroff 
et al., 1995; Regnault et al., 1995; Renard et al, 1995). 
Similar effects are also observed in some S = I spin 
chains like PbNi2V205 (Uchiyama et al., 1999) oi S = ^ 
chains like Sr2Cu03 (Kojima et al, 2004). Even though 
the physics exhibited in these various substituted sys- 
tems is very rich, it is not directly relevant to the aim of 
this review as it is governed by subtle deviations from a 
purely one- dimensional Heisenberg situation. 



IV. EXPERIMENTS ON POINT DEFECTS IN 2D 
SYSTEMS: THE METALLIC CUPRATES 

A. Overview 

The work mentioned up to now on ID or quasi ID 
systems has been devoted mostly to insulating spin sys- 
tems, as the major available experimental systems are 
non-metallic. The most direct extension to 2D of the 
work sketched above for spin chains would be to study the 
response to impurities in the case of a spin 1 /2 Heisenberg 
plane, for which theoretical computations along similar 
lines are available (see Section II). This situation might 
be encountered in the case of the undoped cuprates, 
e.g. La2Cu04 and YBa2Cu306 which are good Mott- 
Hubbard insulators. We shall explain in section IV. C 
why this simple experimental extension has not been suc- 
cessfully performed so far. On the other hand, in the 
cuprates, hole or electron doping is easy to achieve and 
allows one to span the phase diagram of Fig. 2 displayed 
in the introduction. The most important work done so 
far using defects or impurities to probe 2D systems has 
thus been focused on the doped layered cuprates, among 
which are the HTSC phases. 

Before beginning with a detailed description of the ex- 
periments, it is essential to recall briefly here the impor- 
tant features which characterize the different phases of 
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Fig. 2. Here the loss of AF order docs not give way 
immediately to a metallic and superconducting ground 
state, but rather to an intermediate regime in which a 
disordered magnetic ground state, usually described as a 
spin glass (SG) phase separates the AF domain from a 
superconducting region which takes the form of a dome. 
One denotes respectively the left and right of the dome 
as the underdoped and overdoped regime, and refers to 
the doping at which the maximum Tc occurs as optimal. 
An essential feature which has been demonstrated exper- 
imentally is that the superconducting state does not have 
a uniform gap amplitude in k space but is rather c?-wave, 
i.e. it has a maximum gap in the CuO bond axis di- 
rections and nodes in the diagonal directions (Tsuei and 
Kirtley, 2000). Now all this would have been relatively 
simple if the metallic state above the superconducting 
regime were a simple metal describable within a Fermi 
liquid theory. In fact this limit seems to be achieved only 
for very large doping, on the far right of the phase dia- 
gram, i.e. as far away from the Mott AF state as possible. 
Directly above the superconducting dome the transport 
properties, e.g. a T linear dependence of the resistivity 
with large scattering rates « ksT/h at optimal dop- 
ing, appear to be at odds with a simple Fermi liquid de- 
scription, hence the designation "strange metal" (Varma 
et al., 1989). Most importantly, in the underdoped regime 
the magnetic properties exhibit a quite unusual behavior, 
with a spin susceptibility revealed by NMR Knight shift 
data which exhibits a large drop at low T, and is nearly 
suppressed before the superconducting singlet state is 
even established (AUoul et al., 1989). This drop in the 
susceptibility occurs as if a gap were opening in the ex- 
citation spectrum of the system, justifying the "pseudo- 
gap" designation introduced since 1989, as a small den- 
sity of states still remains at T^. The corresponding gap 
in the spin excitations has indeed been detected by spin 
lattice relaxation data(Yasuoka et al, 1989) and neutron 
scattering experiments (Rossat-Mignod and et al., 1991). 
This pseudogap, which has been observed as well some- 
what later by specific heat, STM, ARPES etc... (see 
(Timusk and Statt, 1999) for a review on the pseudo- 
gap) opens at a temperature T* which drops sharply 
with increasing doping, defining a crossover line in the 
phase diagram of Fig. 2. This pseudogap is most proba- 
bly intimately linked with the correlated nature of these 
systems, and its actual physical origin is still intensely 
debated. One class of interpretations could be the estab- 
lishment of a hidden order disconnected from supercon- 
ductivity such as a spin ordering ( Resonating Valence 
Bond (RVB) (Anderson, 1997) or d-density wave (ddW) 
(Chakravarty et al., 2001)), or involving orbital currents 
(Varma, 2006), or a charge segregation into stripe or- 
der (Carlson et al., 2002). Another interpretation is a 
precursor pairing state, the superconducting phase be- 
ing only established at when the pairs achieve long 
range phase coherence (Emery and Kivelson, 1995). At 
the present writing the understanding of the pseudogap 
state remains controversial. 



In this section we shall describe the experiments on 
point defects performed in the normal state of the HTSC 
systems, which cover many phases of the phase diagram 
introduced hereabove. We shall start with a discussion 
of the defects which can be used (part IV. B). Experi- 
ments which allowed one to probe the magnetic response 
of the Cu02 planes both on a local scale and on the 
macroscopic scale will be reviewed then in part IV. C. 
Since the systems in question are metallic, the defects 
also cause significant modifications of the transport prop- 
erties, which differ dramatically in the various parts of 
the phase diagram, as will be discussed in part IV. D. 
Overall, these data have some profound implications on 
the normal state properties of the cuprates. The most 
prominent qualitative features will be discussed in part 
IV. E, while comparison with theoretical approaches will 
be postponed to Section V. 



B. Controlled defects In and out of the Cu02 planes 

Impurity substitutions are a priori the simplest de- 
fects allowing one to induce well-defined perturbations 
and to "tune" the nature (spin or/and charge) of the 
perturbing potential. The ideal experimental situation 
would be to choose both the impurity and the atomic 
substitution site in the synthesis process. Furthermore, 
one would like to achieve this for any doping level in 
the planes, in order to compare the impurity properties 
in the various phases of the phase diagram. Unfortu- 
nately physical chemistry considerations play a large role 
since competing solid state phases combining the impu- 
rity atom and the constituents of the cuprate may exist. 
Even in the absence of such phenomena, various substi- 
tution sites may occur depending on the charge reservoir 
layers of the cuprate material, and interference between 
the intended substitution and the ions responsible for the 
doping of the Cu02 planes may occur. Efforts have been 
made to control both the doping process and the associ- 
ated disorder effects in the pure cuprates, as well as the 
impurity substitutions and/or defect content. 



1. Hole doping of the cuprates 

Doping the parent insulating cuprates to span the 
phase diagram is usually done by a modification of the 
layers separating the Cu02 planes, either by a heterova- 
lent substitution or by changing their oxygen content. In 
both cases, the ions introduced into the structure create 
a modification of the Coulomb potential which disrupts 
the lattice periodicity and which will be felt as a scat- 
tering potential by the carriers in the Cu02 layers. The 
infiuence of this source of disorder on the properties of 
the cuprates was somewhat neglected in the early studies 
of the cuprates. 

It is clear that dopant disorder is particularly strong if 
the ionic dopant is in the separating layer of multilayer 



20 



cuprates or in the reservoir layer near neighbor of the 
Cu02 planes. For example, Sr^+ /La'^+ in La2Cu04 and 
Ca^+/Y'^+ substitutions in YBCOg induce hole doping 
of the Cu02 planes. These heterovalent substitutions are 
easily achieved, as no other substitution site seems to 
be available, and enable one to span the phase diagram, 
which is found to be quite similar for the two families 
quoted above, but with rather low optimal Tc w 4QK. For 
Yi^yC&y Ba2Cu306, the solubility limit is apparently 
smaller than y « 0.35 required to reach the optimum T^. 
(Casalta et al., 1993). These systems are quite important 
to study the electronic properties of the lightly doped AF 
charge-transfer insulator. 

In other families of cuprates, the doping is produced by 
introduction of extra oxygens in the charge reservoirs, on 
sites located far from the Cu02 planes. The optimal T^s 
found in such cases, i.e. in single layer Hg-1201, Tl-2201 
or bilayer YBCOe+i; and Bi-2212, are higher ( > 90K) 
than those found for the heterovalent substitutions men- 
tioned above. So it has been understood for some time 
that the dopant disorder certainly plays a role in limit- 
ing Tc over all the superconducting dome. This source 
of disorder is indeed revealed in NMR experiments either 
on the NMR widths as a distribution of magnetic hy- 
perfine effects (Bobroff et ai, 1997a) or as a distribution 
of Ti relaxation times (Singer et al, 2002). Recent ex- 
periments in which the various possible substitution sites 
have been systematically used to dope Bi2212 have con- 
firmed elegantly this correlation between the optimum Tc 
and the proximity of the substitution site to the Cu02 
plane (Eisaki et ai, 2004). The existence of dopant dis- 
order in Bi2212 and its influence on the local supercon- 
ducting properties have been recently revealed by STM 
experiments (McElroy et al., 2005). 

Finally, let us point out that the disorder of chemical 
dopants is a well identified source of disorder, but that 
other sources of disorder do occur in the cuprate families. 
For instance, the structural modulations due to the misfit 
between the lattice parameters of the Cu02 plane and 
BiO plane in BSCCO, or the tilt of the oxygen pyramids 
and the stripe structure around x=0.12 in the Ba doped 
lanthanum compound. This type of "disorder "is even 
less controllable and dependent on the family. Because 
we are dealing here with specific impurity effects, we are 
forced to consider these aspects and also to select systems 
in which "intrinsic disorder" is minimal. Here we include 
in the definition of intrinsic disorder both any structural 
disorder arising in the charge reservoir or Cu02 layers, as 
well as that due to the random distribution of the dopant 
atoms. 

We highlight here the particular case of YBCOg+a; for 
which the hole doping of the parent compound YBCOe, 
which has no oxygen in the intermediate pure Cu(l) 
layer, is obtained by inserting oxygens in this layer. As 
can be seen in Figure 16, the charge transfer is only ini- 
tiated when X is large enough to give near neighbor oxy- 
gens. With the growth of chain segments the hole con- 
tent Uh of the planes increases, but does not scale with 



X. Thus, the phase diagram versus x differs markedly in 
appearance from that obtained by heterovalent substitu- 
tions for which rih is controlled directly by the atomic 
substitution. At fixed values of x, ordered oxygen struc- 
tures are formed, the most prominent one being for x = 1, 
with filled Cu-0 chains, which corresponds to slight over- 
doping of the plane (with Tc about 1°K lower than op- 
timal). For a; = 0.5 the ortho-II phase with alternating 
filled and empty chains can be formed and corresponds 
to a good underdoped situation with a well defined large 
pseudogap. But again for x > 0.5 the added oxygens be- 
tween chains (see Fig. 16) do not dope the planes as long 
as two oxygens arc not near neighbors, that is ideally 
in an ordered structure up to x = 0.75, which explains 
the Tc plateau at 60K around a; ~ 0.6 . For these simple 
reasons, x = 1 and x ~ 0.6 are important compositions 
characteristic of the near-optimal and pseudogap phases, 
which can be easily reproduced. Furthermore, at x ~ 0.6, 
the ordering of the chains, and the reduced effect of the 
Coulomb potential from the Cu^+ - 0^~- Cu^+ units, 
with an overall charge equivalent to 2 Cu"'", promotes a 
situation where intrinsic disorder is much less important 
than in most other cuprate systems, as has been high- 
lighted by NMR measurements (Bobroff et ai, 2002). 
The fact that the ordered chain phases are important 
fixed points of the phase diagram has been thoroughly 
investigated by the Vancouver group (Liang et al, 2000). 
YBCO is therefore quite an interesting family for impu- 
rity substitution studies, since once a substitution has 
been performed, the undoped AF, the pseudogap phase 
and the near optimal phase arc accessible with simple 
heat treatments and correspond to phases with weak in- 
trinsic disorder. Y1248 (Karpinski et al., 1988), a vari- 
ant of YBCO with a well ordered two chain layer, is also 
a very interesting case of underdoped pseudogap com- 
pound with little disorder. 



2. Impurity substitutions in the cuprates 

Since the early 1990s, impurity substitutions have been 
thoroughly studied in the two preferred cuprate fami- 
lies YBCO and LSCO, and later, after the development 
of ARPES and STM experiments, data on substituted 
BSCCO have been taken as well. The easiest path to 
ensure that an impurity substitutes on a given site is to 
assume a substitution in the ceramic synthesis and check 
by x-rays the absence of impurity phases which definitely 
occur if the substituent does not enter in phase, or sub- 
stitutes on an alternate site. This allows one to perform 
Rietveld analysis of the powder patterns (sec for instance 
(Maeda et al., 1989)), to measure the variation of lattice 
parameters and to determine the solubility limit in solid 
state. Synthesis of single crystals docs not allow similar 
checks but can be calibrated with respect to results on 
ceramic samples. 

It is important to note that the physical chemistry 
complications linked with substitutions raises serious 
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FIG. 16 An oxygen inserted between two copper sites in a 
vacant chain row only changes the (initially) Cu"*" into Cu^^ . 
This does not transfer any charge to the Cu02 plane. When 
a chain segment with at least two oxygens is formed charge 
transfer to the plane occurs. Full chains as that shown on the 
right of the figure might then form an ordered pattern, with 
specific values of the hole content in the planes, and minimal 
disorder. 



problems if one deals with a system which has some in- 
trinsic disorder, as correlation between the impurity lo- 
cations and the intrinsic disorder do occur naturally, and 
therefore these effects are not additive. So studying im- 
purity effects directly will be easier in families of cuprates 
with minimal intrinsic disorder. 



a. Homovalent substitutions in the separating layer of mul- 
tilayer cuprates. Substitution on the sites by rare 
earths, allowing magnetic probes in the blocking layers, 
have been thoroughly performed in the YBCO family. 
Most rare earth moments apparently interact very weakly 
with the Cu02 planes, however, since superconductivity 
is nearly unaffected even for the pure compounds such as 
DyBaCuO or GdBaCuO (Hor et aL, 1987). Apart from 
applications purposes, the main interest of these substi- 
tutions is to provide new spin probes for the experimental 
study of the magnetic properties of the Cu02 layers. A 
few per cent of Gd substituted allowed thorough studies 
by Gd3+ ESR (Janossy et aL, 1994), while Yb substitu- 
tions allows one to use Yb'^"'" as a Mossbauer effect probe 
(Hodges et aL, 1995). These techniques provide hyperfine 
probes similar to the ^^Y nuclear spin probe, but allow 
the exploration of different T ranges as well as different 
electron spin relaxation time windows. 

Pr acts differently from the other rare earths as it re- 
duces the charge carrier content so that PrBaCuOr is 
antiferromagnetic and resembles the parent compound 
YBCOg. At small Pr content, it has been demonstrated 
that Pr hybridizes more with the O of the Cu02 planes 
than the other rare earths, and that its effect is then to 
reduce the hole content in planes by filling holes on the 
neighboring oxygen (MacFarlane et aL, 2002). In this 
picture, the effective charge of Pr is then closer to 4-I-. 



However, the corresponding perturbation in the plane is 
found to display a magnetic character as well. 



b. Substitutions on the Copper sites The most important 
substitutions are the Zn and Ni substitutions on the 
Cu(2) sites of the Cu02 planes which allows one to re- 
place the Cu 3d^ S — 1/2 site by nominal S — Q and 
5=1 impurities, respectively. These substitutions have 
been easily performed in LSCO as they only take place 
on the single Cu site of the structure. But as this system 
contains significant intrinsic disorder due to Sr^+/La^+ 
substitution, correlated disorder effects might severely 
modify the influence of the impurities. 

In YBCO, the essential difficulty comes from the possi- 
ble substitutions on either the chain or plane Cu sites. It 
is easy to ensure that an impurity substitutes for Cu, by 
checking with x-rays the absence of impurity phases. But 
once it is ensured that up to 4% Zn and 3% Ni do substi- 
tute for Cu, thorough studies are required to determine 
on which site this substitution occurs. Various experi- 
mental investigations using EXAFS (Maeda et aL, 1989) 
or XANES (Yang et aL, 1990) suggest that Zn substitu- 
tion preferentially occurs on Cu(2) in the Cu02 planes 
while Ni substitutes on both the chain Cu(l) and pla- 
nar Cu(2) sites. In the case of Zn the observation by 
NMR of the ^^Y near neighbors of Zn allowed one to 
ensure, from the integrated intensity of this line, that 
more than 90% of the nominal Zn content substitute on 
the Cu(2) site, within experimental accuracy (Mahajan 
et aL, 1994). Since such an experiment is impossible in 
the case of Ni, the equal probability of atomic substi- 
tutions on the Cu(l) and Cu(2) sites given by structural 
studies could not be controlled independently Finally, Li 
has also been found to substitute solely on the Cu(2) site 
(Sauv et aL, 1996) in YBCO, with the actual Li content 
being again determined from ^^Y near neighbor intensity 
measurements (BobrofF et aL, 1999). 

Many impurities have been found to substitute prefer- 
entially on the chain Cu(l) site, such as Co^+, Fe^+ or 
Ga^+((Bridges et aL, 1989), and references therein). The 
impurities on the chain site donate their extra electrons 
to the planes and also interrupt the chains, which further 
reduces the charge transfer to the planes (Miceli et aL, 
1989). These 3-|- impurities often do not like the square 
oxygen coordination of the chain site and accomodate 
extra oxygens, which partly compensate the reduction of 
hole doping and maintains a tetragonal structure of the 
compound whatever the oxygen content, as for instance 
for Ga^+(Xiao et aL, 1988). The dominant overall effect 
of such substitutions on the chain site is a reduction of 
the hole doping of the Cu02 plane, making the optimally 
doped composition impossible to reach. In many cases, 
even for Ni (Adachi et aL, 2000), it has been found that 
heat treatments in controlled atmosphere do allow par- 
tial transfers of impurities between the Cu(2) and Cu(l) 
sites. 
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c. Defects induced by electron irradiation. To control de- 
fect creation in a solid, one may use irradiation by ener- 
getic particles, preferably light ones which have a large 
enough penetration length to damage bulk samples. 

Ions or even protons produce a large density of defects 
which recombine into clusters of defects. In particular 
high energy heavy ions yield linear tracks of heavily dam- 
aged material which can be used to pin vortices in the SC 
state (Civale et ai, 1991; Konczykowski et ai, 1991). To 
produce isolated point defects, the most suitable method 
is to irradiate with high energy (MeV) electrons at low 
T, optimally in liquid H2 at 20K, which avoids recombi- 
nation of the point defects created. In cuprates, this al- 
lows the production of Cu and O vacancies in the planes, 
the ejected Cu and O being trapped in interstitial lo- 
cations out of the Cu02 planes. Upon heating the sam- 
ple, some of the created defects recombine, the oxygens 
first as their mobility is larger (Giapintzakis et al., 1994; 
Legris et al., 1993; Tolpygo, 1994). The complete set of 
studies performed so far indicate that the damage is very 
homogeneous on samples of thickness < 50 /im, and that 
the vacancies in the Cu02 planes dominantly modify the 
transport and superconducting properties, these effects 
being even quantitatively identical to those produced 
by Zn substitutions (Rullier-Albenque et al., 2000). If 
this method applies ideally for transport studies in single 
crystals of small thickness, it unfortunately requires too 
much electron beam time for the large amount of sample 
required for bulk or NMR experiments. 



C. Magnetic properties induced by in-plane impurities 

1. Impurities in the AF phases 

The parent compounds of the cuprates appear to be 
good Mott insulators, with a very large Heisenberg ex- 
change interaction J ~ 1500K (Lyons et al., 1988). This 
is revealed by neutron scattering experiments which dis- 
play an inelastic scattering peak at the AF wave vector 
(tt, tt) . The q-width of this peak gives a quantitative es- 
timate of the dynamic AF correlation length S,(T), which 
fits perfectly at high T in La2Cu04 the exp(J/fcBT) 
variation expected for the Heisenberg 2D square lattice 
(Birgenau and Shirane, 1989). However, while a spin 1/2 
Heisenberg 2D system should only order magnetically at 
T — 0, the undoped cuprates display a well established 
3D AF order with Neel temperatures Tjv , which should 
be largely dependent on the magnitude J± of the inter- 
planar couplings. As this is not the case (AUoul et ai, 
1990; Mendels et al., 1991) it has become clear that the 
ordering is governed instead by a deviation of the plane 
Hamiltonian from the ideal Heisenberg situation. Indeed 
a very small XY anisotropy of the exchange interaction 
drives a Kosterlitz-Thouless transition in the 2D system 
at a temperature TffT which is slightly modified by 3D 
interactions (AUoul, 1991; Keimer et al., 1992). 



a. Substitution effects in the undoped phase. What hap- 
pens when substituting a Cu by a nonmagnetic or mag- 
netic impurity? In YBCOg magnetic impurities substi- 
tuted on the chain Cu sites are only found to change 
the coupling between bilayers in such a way that planes 
of neighboring bilayers switch from an AF ordering to 
a ferromagnetic ordering, still retaining the bilayer AF 
order (Dooglav et al., 1996; Lutgemeier, 1988). But Tjv, 
which is fixed mainly by Tkt, is not modified signifi- 
cantly by this change of 3D magnetic structure of the 
ground state. For a nonmagnetic impurity such as Zn, 
which substitutes on the planes, the dominant effect is 
the dilution of the magnetic lattice, which has been has 
been studied (Vajk et al., 2002) and is found to reduce 
Tjv moderately at low concentration as a result of the re- 
duction of ground state energy and of Tkt- Overall, the 
influence of substituted impurities on the macroscopic 
properties of the ground state do not tell us much about 
the single plane physics. 

In this AF state, local measurements of the modifi- 
cations induced by a nonmagnetic defect are somewhat 
difficult. The main observation done so far has been the 
reduction of intensity of the ''"'Cu NMR signal detected in 
the internal magnetic field without external applied field 
(Zero Field NMR: ZFNMR) (Mendels and Alloul, 1988). 
This gives an indication of the spatial range over which 
the impurity induces a sizable modification of the AF or- 
der (Mendels et ai, 1990). So far these experiments did 
not enable a detailed study of the local modifications of 
the magnetic order, as their analysis depend on the modi- 
fications of the various parameters mentioned above and 
are not solely associated with the single plane physics 
and its spin liquid state. Experiments above Tjv in the 
paramagnetic state allowed one to detect modifications of 
the spin dynamics from NMR spin lattice relaxation data 
(Carretta et al., 1997), but again no local study of the 
magnetic properties around impurities could be achieved. 



b. Impurities in the doped antiferromagnetic state. Let us 
first recall the effect of small hole dopings ( below 10%) 
in the "pure" systems. One important fact learned early 
on from studies on LSCO is that carriers abruptly re- 
duce Tn, less for electron doping than for hole dop- 
ing (Birgenau and Shirane, 1989), but much more than 
spinless homovalent substitutions, such as Zn^+. A sim- 
ilar decrease of Tjv occurs for hole doping by Ca in 
Yi_yCayBa2Cu306 (Casalta et ai, 1993). In any case 
these carriers are mobile at high T as the planar resis- 
tivity increases above room T as in a metal (Ando et al., 
2004) . They should localize at low T in the Coulomb po- 
tential of the dopant, so that non generic modifications 
of the AF state are expected. Although detailed studies 
have not been performed so far, suppression is how- 
ever achieved for a comparable doping of Uh — 0.03 for 
Sr and Li, which correspond respectively to out of plane 
and in plane coulomb potentials (Suh et ai, 1998). It has 
been generally seen that a new "phase" appears below 
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T/v , which contains some amount of static disorder, and 
which has some characteristics of a mixed AF/spin glass 
phase (Chou et ai, 1993). 

For hole dopings larger than Uh — 0.03, the ordered 
magnetic state disappears and gives place to a static dis- 
ordered magnetic phase detected by magnetization mea- 
surements, /iSR or NMR (Cho et al., 1992; Chou et ai, 
1995; Niedermayer et ai, 1998). This latter phase has 
many characteristics of a true spin glass (SG) phase. In 
some systems like the lanthanum family, it extends over 
a large doping range, between 0.03 and 0.07 but still ex- 
ists far into the SC state nearly up to optimal doping 
X = Uh = 0.15, albeit with a low spin glass temperature 
Tg (Niedermayer et ai, 1998). Whether the tendency to 
stripe formation at a; = 0.12 in the Ba doped case favors 
a coexistence of the SG and SG phases on a microscopic 
scale, or leads to a phase separation is not fully clari- 
fied. In YBGOg+a; the SG phase has been occasionally 
detected for intermediate oxygen contents x ~ 0.45 but 
has totally disappeared for x ~ 0.6, that is rih — 0.10. So 
it is clear that the SG phase is not generic in the phase di- 
agram of the "pure systems" , its extension being linked 
with the intrinsic disorder present in the systems. In a 
specific family of bilayer cuprates CaBaLaGuO, system- 
atic studies reveal a clear relation between the Tg and 
Tc values in the phase diagram, which suggests that the 
influence of the intrinsic disorder in these systems affects 
these temperatures similarly (Keren and Kanigel, 2003). 

The introduction of spinless impurities adds to the ex- 
isting disorder in these phases. It is remarkable that Zn 
substitution in YBCO extends significantly the SG phase 
range in this system (Mendels et ai, 1994a)(see Fig. 17), 
so that the phase diagram of YBGO:Zno.o4 resembles that 
observed in nominally pure LSCO, or in Yi-j^CayBCOg. 
This suggests again that extrinsic or intrinsic disorder 
may have similar influences on the properties of these 
systems at least in the low hole doping range. 

Regarding the influence of controlled impurities in the 
intermediate doping range, let us once again emphasize 
that we are not dealing with single impurity effects, and 
that the interferences with the intrinsic disorder of the 
"pure" systems complicate the situation markedly. As 
will be seen in Section VI. E, Zn impurities in LSGO 
around n^^O.lO obviously nucleate a long range mag- 
netic order detected by neutron scattering experiments 
which is not present in the pure system. This kind of 
order induced by disorder bears some analogy with the 
situation encountered in spin chains for which magnetic 
order can be enhanced by impurities (see Section III.D.5). 
Although many qualitative and materials aspects are re- 
vealed by the experiments involving impurity substitu- 
tions, such observations are even harder to analyze the- 
oretically than those done in the undoped AF case (see 
however Section V). 
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FIG. 17 The phase diagram of pure YUCOe+x (full sym- 
bols) only exhibits a small range of oxygen contents, near the 
tetragonal to orthorhombic transition at a; = 0.4 for which 
disordered static magnetism occurs at low T. With Zn sub- 
stitution, at a concentration of 4% (open symbols), the su- 
perconducting range is highly depressed, the local magnetism 
induced by Zn drives the appearance of a large spin glass (SG) 
regime. The phase diagram becomes then similar to that of 
Fig. 2. Notice that the onset of pseudogap opening T* is not 
modified by Zn substitution (Figure taken from (AUoul et ai, 
2000)). 



2. Impurities in the pseudogap phase above Tc 

When the doping is large enough to destroy the static 
AF phase and to establish a metallic and superconduct- 
ing state, the magnetic properties display, in the under- 
doped regime, a "pseudogap" (see Fig. 2). It opens up 
at a very high temperature T* (« 350 K in YBCOe.e) 
and is signalled, e.g., by a decrease in the spin suscepti- 
bility (see Fig. 18). Here we are in an interesting doping 
range where the system is a metal with magnetic corre- 
lations, which is a hard theoretical problem. This is the 
regime for which the use of impurities to probe the physi- 
cal properties was expected to be the most fruitful. Such 
an approach was started in the early 1990's and provided 
first basic qualitative information. Improvement of the 
samples and experimental techniques have since allowed 
one to obtain more refined quantitative information on 
the perturbations introduced by the impurities. A major 
effort has been devoted to the case of Zn, Ni and Li im- 
purities substituted for Cu in the YBaGu Oq+x system 
specifically for x ~ 0.6 for which the intrinsic disorder is 
minimal. 



a. Influence of impurities on the pseudogap The first point 
which was investigated in the early experiments was the 
influence of the substituted impurities on the tempera- 
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FIG. 18 The ''O NMR shift versus T for YBCOe.e pure and 
substituted with 2% Zn and 2%Ni. The NMR shift for the sites 
far from the Zn is found totally independent of the substitu- 
tion, with the onset of opening of the pseudogap remaining 
at «350K for all samples. 
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FIG. 19 The 

tuted YBCOe, 

also seen in the inelastic neutron scattering data taken at the 
AF wave vector (7r,7r) shown in the inset. Both data are dom- 
inated at low T by the contribution of sites near the Zn, but 
still display an evidence for the onset of the pseudogap as a 
high T maximum of the dynamic susceptibility around 150K, 
seen in the pure samples (Julien et ai, 2000). 



ture T* at which the pseudogap opens. As can be seen in 
Fig. 18, the NMR shift of the nuclei located far from the 
impurities is not affected by the impurities Zn or Ni, even 
when the impurity content is large enough to completely 
suppress superconductivity (AUoul et ai, 1991; BobrofF 
et ai, 1997b; Riseman et at, 1994). 

Most probes sensitive to the normal state pseudogap 
far from the defects confirmed these finding later on, and 
showed that they applied to any kind disorder, such as 
irradiation defects or Li impurities. Some early inelastic 
neutron scattering (INS) (Kakurai et ai, 1993), and ^^Cu 
Ti experiments (qing Zheng et al, 1993) claimed a filling 
of the pseudogap in the spin excitations at the wave vec- 
tor (tt, tt), in the presence of impurities. However more 
detailed data taken later have confirmed that the Zn im- 
purities create low energy midgap excitations (Bourges 
et ai, 1996; Harashina et ai, 1993). These excitations, 
which will be described in detail in part IV. C. 4, are as- 
sociated with the fluctuations of the local moment AF 
staggered response spatially localized near the Zn. They 
partially fill in the pseudogap which characterizes the be- 
haviour far from the Zn; this is apparent in both the ^^Cu 
Ti data and the INS, as can be seen in Fig. 19. Indeed, 
NQR investigations (Itoh et al., 2003) in which the reso- 
nances of the sites far from the impurity could be singled 
out clearly reveal the persistence of the pseudogap for 
the Ti data of these sites (see Fig. 20). 

So this body of experiments give consistent evidence 
that the pseudogap is not modified far from the Zn site, 
while the local magnetic properties arc strongly modified 
near the Zn location. This is also confirmed by other 
experiments which also allow one to detect the onset of 
the pseudogap from macroscopic data, such as specific 
heat (Loram et ai, 1994), transport for Zn substitutions 
(Mizuhashi et ai, 1995; Walker et ai, 1995), and also 
irradiation defects (RuUier-Albenque et ai, 2003). 
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FIG. 20 The Cu nuclei far from the Zn impurities display 
in YBC01248 the same reduction of (TiT)"^ as that found in 
the pure compound below the pseudogap temperature (Itoh 
et al, 2003). 



b. Local moments induced by Zn. The first question 
which arises is whether a nonmagnetic site induces a free 
paramagnetic moment in a metallic correlated system, as 
was seen in the case of undoped spin chains and ladders. 
Ideally, this might be probed by macroscopic susceptibil- 
ity measurements, or by electron spin resonance (ESR) 
experiments. A theoretical proposal for the occurrence 
of a local moment with an ESR experiment in LSCO:Zn 
to support this view was given (Finkelstein et at, 1990). 
Although the theoretical idea was quite valuable, the de- 
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tected signal could not be associated with isolated Zn im- 
purities, as will be shown below in part IV. C. 4. Indeed, 
the early experiments were plagued by materials prob- 
lems. For instance, ESR signals detected in YBCO sam- 
ples were exclusively due to the so-called "green phase" 
Y2CUO5, which also induced Curie terms in the suscepti- 
bility of "pure" samples. So, although some initial suscep- 
tibility data on Zn substituted YBCO samples did sug- 
gest evidence of Curie contributions, they could not be 
considered as proofs of the existence of impurity induced 
paramagnetism (Cooper, 1991). The actual contribution 
of substituted impurities to the magnetic susceptibility 
has only been determined since, in very carefully con- 
trolled samples, as we shall see below. 

The first indirect but unambiguous evidence that Zn 
induces a dynamic paramagnetic moment in an under- 
doped cuprate was obtained by monitoring the ^^Y NMR 
linewidth in YBC06.6:Zn.y (AUoul et at, 1991). A non- 
magnetic impurity in a metallic host normally only yields 
a T independent broadening of the unshifted host NMR 
signal due to Friedel oscillations of the density of states. 
By contrast, the unshifted ^^Y NMR has a linewidth 
which is found to increase markedly at low T, while the 
susceptibility of the host decreases due to the pseudogap. 
This reveals that the spatially dependent spin susceptibil- 
ity induced by the Zn impurities displays a T dependent 
paramagnetism which increases at low T and is therefore 
"Curie-like", as illustrated in the case of a spin chain 
(see Fig. 14). The only difference is that here the signals 
of the different near neighbor shells of the impurity are 
overlapping and just yield a broadening of the line. The 
actual T dependence of the susceptibility could not be as- 
certained by such a measurement, since the relationship 
between the magnitude of the induced spin polarization 
and the magnetic susceptibility is not straightforward, as 
will be clear from the later work described below. That 
a Curie like moment developed in the vicinity of the Zn 
was clearly confirmed by resolving with dilute samples 
the satellite NMR signal of the ^®Y near neighbor (n.n.) 
sites of the substituted Zn (Mahajan et ai, 1994) (see 
Fig. 21). This provided the first local detection of the 
field induced paramagnetism near the Zn, well before the 
equivalent information could be monitored in the case of 
spin chains. These data imply that the spin polarization 
of the Cu n.n. to Zn is already at 100 K more than 
ten times larger than that of the pure host, so this is not 
a mere minor modification of the host density of states, 
but a strong effect of the electronic correlations of the 
system, similar to that observed in spin chains. The very 
short spin lattice relaxation Ti of the ^®Y n.n. NMR at 
low-T is a further striking evidence which establishes the 
occurrence of a local moment induced by the spinlcss Zn 
impurity. Quite analogous ^^Y n.n. NMR data were ob- 
tained later on Zn substituted YBC01248, which is also 
a low disorder underdoped cuprate, thus confirming the 
generality of these observations (Williams et ai, 1998). 

This Zn induced paramagnetic contribution to the 
macroscopic susceptibility Xc could then be determined 
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FIG. 21 The NMR spectrum of YBC06.6;Zn (or equiva- 
lently YBCOe.s^Li) exhibits a main line and two satellite res- 
onances as shown in the inset. Their shifts ^^K are plotted 
versus T as empty and full symbols respectively for Zn and 
Li substitution. While the main line shifts (circles) retain 
the characteristic pseudogap T variation of the pure sample 
(crosses), the n.n. satellite shifts depart markedly at low T, 
the difference of shift with the main line displaying a Curie 
variation. This results from the large Curie like spin po- 
larization which appears on the four Cu sites around the Zn 
or Li, two of them being directly sensed by the *^Y n.n. to the 
Zn (see Fig.22). The ''Li NMR shift (cross in a circle, right 
scale) directly probes as well this Curie polarization of its 
four Cu neighbours (figure constructed from (Bobroff et al., 
1999)). 



from SQUID magnetization data taken on samples in 
which the fraction of parasitic impurity phases could be 

minimized and remained only as traces (Mcndels et al., 
1994b; Zagoulaev et al, 1995). For x = 0.66, a Curie- 
Weiss law 

Xc = f^lff/[ikB{T + e)] 

could be associated with Zn induced magnetism down to 
T wlOK, with a smah value < 4K for a; = 4% and an 
effective moment = l.O/is (Mendels et al., 1999). 

With the known hyperfine couplings for the ^^Y nu- 
cleus, the spin polarization on the 4 Cu n.n. to the Zn 
would correspond to a total susceptibility of the same 
magnitude as Xc- This was initially taken as an indica- 
tion that the moment was mainly localized on the 4 Cu 
first n.n. to the Zn (Mahajan et ai, 1994). But the data 
only allows one to conclude that the sum of the staggered 
magnetic susceptibility on more distant sites roughly can- 
cel out and contribute little to Xc (Mahajan et al., 2000). 

The ^^Cu nuclear spin itself, which is more sensitive to 
the on site magnetization, allows one to probe this long 
distance staggered magnetism (Julien et ai, 2000; Walst- 
edt et al, 1993). The large increase of the ^^Cu linewidth 
induced by Zn established the existence of long distance 
spin polarization effects but did not allow to provide 
enough quantitative information to correlate them with 
the short distance information gained from ^^Y NMR 
data. 
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FIG. 22 The Cu 3d^2_j^2 and 02p^, 02py orbitals of the 
Cu02 plane are represented, with a Li or Zn impurity substi- 
tuted on one Cu site. The magnetization induced on each Cu 
site by the applied field is staggered and decreases in magni- 
tude to reach the very small uniform value found in the pure 
compound. In the inset we display the hyperfine coupling 
paths linking the different nuclear spin sites with their near 
neighbor Cu orbitals (also with the on site orbitals for Cu). 
The Y sites are below the Cu02 plane (and above the other 
Cu02 plane of the bilayer). 



FIG. 23 The difference A;/ of ^^O NMR linewidth between 
the substituted and pure YBCOe.s samples gives a measure 
of the width of the histogram of local fields induced by Zn or 
Ni on the ^'^O nuclei. The product TAj/ increases markedly 
at low r, while it should remain constant if the spin polar- 
ization had a mere Curie l/T variation. This is evidence that 
both the magnitude of the local fields and the width of their 
histogram, that is the decay length ^imp of the spin polar- 
ization, increase sizably at low T (data taken from (Bobroff 
et al, 1997b) and (Ouazi et al, 2004)). 



c. The Ni impurity: a weakly coupled local moment and the 
Rrst evidence for a large ^imp- The Ni impurity, which is 
in a 3d® state, was expected to display significantly dif- 
ferent magnetic properties owing to its on site local mo- 
ment. In this case, the macroscopic spin susceptibility 
Xc per Ni site was found to be indeed larger than for Zn, 
with /ie// = 1-7/iB for x = 0.66 nearer from that ex- 
pected for a spin S* = 1/2 than for a spin 5* = 1, so that 
the local moment is somewhat reduced with respect to 
that expected for a free ion (Mendels et al., 1999). A 
crude explanation would be that the 3dx2_y2 component 
of the Ni moment couples to the Cu02 band and would 
roughly replace the Cu spin, while the 3d^2_r2 orbitals 
would contribute to a spin 1/2 local moment weakly cou- 
pled with the Cu02 band. This conjecture is supported 
by the evidence that the broadening of the NMR 
(Riseman et ai, 1994), is dominated by the dipole field 
of the Ni local moment, and is only slightly affected by 
the staggered magnetism on the Cu sites, no n.n. 
line being even detected in that case. 

On this basis, it was clear that data on the ^^O NMR, 
which has larger hyperfine couplings with the Cu(2) mag- 
netization, would be less sensitive to the direct dipolar 
local fields induced by the local moment and would give 
more information on the staggered paramagnetism. The 
interesting result obtained was that the ^^O linewidth 
did not increase at low T as the Curie l/T behavior 
displayed by Xc but much faster (see Fig. 23). This was 
considered a proof that in contrast with expectations for 
RKKY oscillations in a non-correlated electron host (Bo- 
broff et ai, 1997b), the T variation of the linewidth is not 
solely due to that of the magnitude of the staggered sus- 
ceptibility but also that its decay length ^imp around the 



impurity has a large T variation (Morr et ai, 1998), (Bo- 
broff et ai, 1998). This has been found to occur as well for 
Zn impurities as can be seen in Fig. 23, so that it appears 
more specific to the host than to the impurity response. 
To obtain quantitative determinations of^imp, the Li im- 
purity which is a spinless impurity, as Zn, appeared more 
practical technically, as will be shown below. 



d. The Li spinless impurity: comparison with Zn. Li has 
been found to substitute on the Cu(2) sites (Sauv et at, 
1996) in YBCO, but detailed comparisons with Zn could 
only be performed when it was established that the ®^Y 
and -'^^O NMR data in YBCOe.e^Li samples were nearly 
identical to those found for YBCOe.e^Zn, as shown in 
Fig.21 (Bobroff et ai, 1999). This allowed one to under- 
stand that a Li"*" spinless impurity which introduces a 
hole in the plane induces the same local magnetic prop- 
erties as the spinless isovalent substitution of a Zn^+ on 
a Cu^+ site. The fact that the ®^Y n.n. sites of the Li 
and the Zn isolated impurities display almost the same 
shift is evidence that the charge difference is not screened 
on the atomic scale. The hole donated by Li to the band 
is rather delocalized, and of course changes the doping 
slightly, but this is experimentally a negligible effect as 
long as one considers dilute concentrations of Li. So the 
magnetic response of the Cu02 plane to the local pertur- 
bation is dominated by the spinless character of the im- 
purity which disrupts the AF interactions with the neigh- 
boring sites, while the extra Coulomb potential of the Li 
impurity is negligible in that respect. 

The great advantage in using Li as an impurity is that 
the large gyromagnetic ratio of the ^Li nuclear spin and 
its small quadrupole moment allows one to detect its 
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NMR signal, its shift being then, as seen in Fig. 21, a 
direct measure of the magnetization of its four Cu n.n. 
sites. This has then ahowed three important technical 
advantages beyond what was permitted by the study of 
the n.n. ^^Y NMR: 

• the very high sensitivity of the ^Li NMR permitted 
accurate measurements of x(^) of the Cu n.n. of 
the spinless impurity. In the underdoped samples, 
this variation is found to display a perfect Curie T 
dependence (see Fig. 21), which confirmed the ob- 
servation made from the NMR that the impu- 
rity induced state behaves as a nearly free param- 
agnetic moment (Bobroff et ai, 1999). 

• the dynamics of the local moment could be stud- 
ied through the ''Li nuclear spin relaxation time Ti 
data (MacFarlane et ai, 2000). 

• the possibility to use in the same sample the ^Li, 
^^Y , ^^O and ^^Cu nuclear spin probes has en- 
abled the quantitative determination of the spatial 
structure of the induced polarization, that is its 
magnitude and £,imp (T) (Ouazi et ai, 2004). 

This latter point will be now discussed in detail, both 
in the underdoped and optimally doped parts of the 
cuprate phase diagram. 



e. Spatial extent of the staggered spin polarization. While 
ideally one would like to measure directly the local mag- 
netization on the nuclear spins of the magnetic sites 
themselves, i.e. Ni sites in the S = 1 chains or Cu 
in the ID or 2D systems, this happens to be in most 
cases hampered by technical limitations linked with the 
very large hyperfine couplings on those sites. We have 
shown in the previous paragraphs that other neighbor- 
ing nuclear probes such as can help to extract some 
information in a less direct way. While in ID YBaNiO 
chains the ^^Y nuclei do probe the magnetization of a 
single Ni moment in a chain, the information is not as 
straightforward with a nucleus which has many different 
Cu neighbors, such as ^^O which feels the response of 
two nearest neighbor Cu sites in YBCO. This results in a 
form factor which connects the field induced on a nuclear 
site to the local magnetization in its vicinity. Although 
this appears as a complication at first sight, it is on the 
contrary an advantage since the availability of different 
nuclear spins allows one to probe various spatial ranges 
of the induced spin polarization. In particular, when the 
magnetization has a staggered spatial dependence it is 
clear that, for nuclei which have z Cu near neighbors in 
a plane, the cancellation of the local field sensed is more 
efficient for larger z (see Fig. 22). One can then un- 
derstand that the range of the spin polarization sensed 
is larger for ^'^Cu, which feels the on site polarization, 
as compared to ^^0(z = 2), further compared to ^^Y 
(z = 4). 
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FIG. 24 The spatial variation of the model spin densities 
which fit the NMR data on '^Li, **^Y and ^'^O nuclei in 
YBCOe.e are represented as a function of distance to the im- 
purity site along the (10) direction in the Cu02 plane. It 
can be seen that the local magnetization magnitude and de- 
cay length ^imp both decrease with increasing T. For a large 
number of sites the local magnetization is much larger than 
that of the pure sample which is the limiting large r value at 
160K, and which is even smaller at 80K ( constructed from 
results of (Ouazi et al, 2004)). 

The data given by the NMR shifts of '^Li , of the n.n. 
satellite resonances of ®^Y, and the lineshapes of the ®^Y 
and ^^O NMR main lines could be compared with nu- 
merical simulations made for model spin polarization ra- 
dial dependences (Ouazi et ai, 2004). This analysis gave 
evidence that the staggered spin density decays exponen- 
tially a few lattice constants away from the impurity, and 
that the induced spin magnetism is not singularly high on 
the near neighbors of the Li (or Zn ) impurities. Similarly 
to the chain case, the local moment perturbation has a 
spatial extent which is characterized by a single range 
parameter ^imp and the macroscopic moment appears as 
the sum of this oscillating polarization. Performing such 
a data analysis at different T allowed one to confirm that 
S^imp and the magnitude of the spin polarization have 
large T variations, as shown in Fig. 24. 



3. From the underdoped to the optimally doped case 

The same approach can in principle be applied for 
increasing hole doping in the Cu02 plane. In the pure 
compounds, the pseudogap temperature decreases and 
approaches Tc. What happens then to the magnetic re- 
sponse? Does it change drastically above T*? This ques- 
tion has of course a direct relationship with the important 
question still pending regarding the phase diagram of the 
cuprates: does the pseudogap line cross the Tc line , giv- 
ing rise to a quantum critical point near optimal doping, 
or does it merge with the line? We shall examine this 
point first. 
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a. Pseudogap line and QCP? Experimentally, the pseudo- 
gap temperature T* has been found to drop sharply with 
increasing doping and approaches the superconducting 
dome around optimal doping. As we have seen above, the 
introduction of impurities reduces Tc and allows one to 
probe the normal state properties at temperatures below 
Tc of the pure compound. As impurities do not modify 
the pseudogap in the underdoped part of the phase di- 
agram, a reduction of Tc induced by impurities should 
open a window allowing the study of the continuation of 
the T* line in the optimally doped regime. NMR stud- 
ies performed so far on substituted samples of YBCO7 
do not show any evidence for a discontinuous change of 
the NMR shift of any nuclear spin species of the material 
between the Tc of the pure compound and that of the 
substituted one. Even with 4% Zn, the ^^Y NMR shift is 
found to remain constant when Tc has been reduced to 
55K (AUoul et ai, 1991). So the existing experiments do 
not allow to conclude whether the pseudogap line crosses 
the Tc line, although some authors report data points in 
their phase diagram which would suggest the opposite 
(WiUiams et ai, 1997). 



b. Spinless impurity magnetism. Experimentally, al- 
though no nearest neighbor ^^Y NMR resonance could 
be resolved in YBCOyiZn, the early NMR data (AUoul 
et ai, 1991) provided evidence for a low T increase of the 
^^Y NMR linewidth which could not be associated with 
impurity phases. In view of the analysis done above, 
this was already clearly evidence that magnetic correla- 
tions persist, since in a nonmagnetic metallic system one 
would expect a T independent linewidth when the spin 
susceptibility is Pauli-like as in YBCO7. Indeed, the non- 
Korringa T dependence of the Cu(2) Ti was also a direct 
proof that magnetic correlations remain important (Taki- 
gawa et ai, 1991). Most evidence in the SC state to be 
discussed in Section VI below also indicate that below Tc 
some paramagnetic moment character appears. But this 
qualitative evidence cannot be connected easily to the 
more quantitative analysis done above in the underdoped 
normal state. The first more quantitative point came 
from SQUID data on very clean samples of YBCOyiZn 
in which the impurity induced Xc was found much smaller 
at low T than in the underdoped compound, which in- 
dicated a weaker magnetism, and certainly not a free 
S = 1/2 moment (Mendels et ai, 1999). 

Subsequently, ^Li NMR data allowed one to measure 
x(T) directly on the 4 n. 71. Cu, which is shown in the left 
panel of Fig. 25 to deviate from a Curie law. The Curie- 
Weiss behavior found in YBCOywith a sizable 8 only 
weakly dependent on Li content ensures that this results 
from the isolated impurity response (Bobroff et ai, 1999). 
Again, comparison of the ^^O and *^Y NMR shapes and 
widths data leads one to conclude that there is no ap- 
preciable difference between the Li and Zn induced mag- 
netism, while Ni has a larger magnetic susceptibility Xc 
but, as in the underdoped case, induces a comparatively 
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FIG. 25 Left: ^Li NMR shift plotted versus T"^ for a series of 
YBCOy samples as well as on an overdoped Ca^"'"/Y^"'" substi- 
tuted sample for y = 7. The data are fitted with Curie Weiss 
laws. The corresponding Weiss O values are plotted on the 
right bottom panel versus oxygen content y, and are seen to 
increase sharply near optimal doping. There the full curve rep- 
resents the variation of Tc in the pure samples. Right top: ef- 
fective moment parameter A^e// of the Curie Weiss fit. (Bobroff 
et ai, 1999). 



smaller in-plane perturbation (Riseman et ai, 1994). The 
careful study of the T dependence of the ^Li NMR shift 
versus hole doping displayed in Fig. 25 provides evidence 
that the Curie- Weiss temperature Q varies abruptly with 
hole doping and becomes nearly zero in the underdoped 
regime. As the effective moment found appears indepen- 
dent of doping, one can conclude that the low T reduc- 
tion of susceptibility in the optimally doped case is due 
to the onset of the energy scale fc^ 8 in analogy with 
the Kondo reduction of local moments in classical metal- 
lic systems. We shall see later on, after performing the 
analysis of ^imp(T) that this analogy also applies to the 
dynamic properties of the local moment. 

The occurrence of a local moment induced by substitu- 
tion of Cu^+ by nonmagnetic AP+ has been established 
as well in optimally doped La2-£i;Sra;Cu04 (Ishida et at, 
1996). They showed that Al^+ exhibits a local moment 
behavior in optimally doped La2_a;Sra;Cu04. Again, the 
^''Al NMR shift exhibits a Curie- Weiss T dependence, 
with a sizable Weiss temperature (8 « 50K). It is not 
clear whether this value of 8 could have been reduced 
with respect to the single impurity limit by the large im- 
purity concentration (3% ) used in this experiment. 



c. Magnetic correlations and ^impiT) The weakness of 
the induced susceptibility, that is the large value of 8 
measured for optimal doping, is quite compatible with 
the absence of detectable ^^Y NMR nearest neighbor res- 
onances within the experimental resolution ( limited by 
the ^^Y linewidth) for Li or Zn in YBCO7. This does not 
prevent the multi-nuclei analysis of the staggered mag- 
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FIG. 26 Temperature dependence of the spatial extension 
^imp of the polarization induced by Zn or Li impurities in un- 
derdoped YBaCuOe.e and optimally doped YBaCuOr. The 
T variation is larger for the underdoped case. (Ouazi et al., 
2004). 



netism, as the hyperfine couplings determined for the n. n. 

in the underdoped case can still be used. Indeed, it 
was shown some time ago that the hyperfine couplings, 
which are based on atomic orbital physics(Mila and Rice, 

1989) , do not depend markedly on the hole content in the 
cuprates (AUoul et ai, 1989; Pennington and Slichter, 

1990) . 

The analysis (Ouazi et ai, 2004) of the T dependences 
of the ^Li shift and of the ^^Y and ^^O NMR spectral 
shapes suggest that the magnitude of the staggered mag- 
netism still decays exponentially a few lattice constants 
away from the impurity. The local susceptibility is still 
about four times larger at lOOK on the Cu n.n. to Zn 
than that of the pure optimally doped compound. This 
establishes the similarity with the staggered magnetism 
induced in the underdoped case, despite claims to the 
contrary ((Tallon et ai, 2002) and see reply by (BobrofF 
et al, 2002)). Although ^imp has a similar magnitude at 
room temperature in the normal and in the underdoped 
cases, £,imp(T) increases much less at lower T in the for- 
mer situation, as Fig. 26 shows. These results demon- 
strated that the local moment susceptibility and ^imp 
have much weaker T variations in the optimally doped 
case than in the underdoped one. The energy scale O 
may control the T variations of both quantities, how- 
ever. Since overdoping corresponds to an increase of 
well beyond the value found for optimal doping (Fig. 25), 
such a scheme would allow a smooth crossover towards 
the Fermi liquid limit for large overdoping. 



4. Spin Dynamics 

We have only considered so far the static response 
of the electronic states to a magnetic field. We might 
wonder what the microscopic process is which drives the 
dynamics of the local moment induced by spinless impu- 



rities? In undoped spin systems, the local moment fluc- 
tuations should be driven by the magnetic excitations of 
the system, if not by a different degree of freedom such as 
a spin phonon coupling. In a 2D metallic system, charge 
degrees of freedom will play a role as well. 

The nuclear spin lattice relaxation Ti of the nuclear 
spins are governed by the fluctuations of the local fields 
induced by the dynamic susceptibility (Pennington and 
Slichter, 1990). Modifications of the Ti rates of the pure 
system due to local moment fluctuations are seen on nu- 
clei even at large distances from the local moment, but 
are usually very hard to analyze due to the distribution 
of couphng constants (Alloul and Bernier, 1973). The 
most reliable experimental situation is that for which a 
well defined nuclear spin site with respect to the impurity 
can be singled out (Alloul, 1977). For instance, the reso- 
lution of a shell of ^^Cu neighbors of Zn by NQR (Itoh 
et al, 2003), showed that the nuclei near the impurity 
display the fastest relaxation rates, as could be expected 
in the local moment picture, contrary to previous claims 
(Ishida et ai, 1993). 

This is the case as well for the n.n. neighbors of 
the impurity (Mahajan et ai, 2000), or for the impurity 
^Li NMR itself in the case of Li substitution (MacFarlane 
et ai, 2000). In these simple cases the contribution to 
1/Ti of the local moment fiuctuations overwhelms that of 
the pure system (Mahajan et ai, 2000), and was studied 
thoroughly for the Li case (MacFarlane et ai, 2000). 

In a simple relaxational model for which the local mo- 
ment is assumed to fluctuate as a single entity at a rate 
T~^, the nuclear spin lattice relaxation rate is related to 
and to the local moment susceptibility Xc by 



(29) 



where /i„ and fiB are the nuclear and electronic magnetic 
moments. Here, A is the effective hyperfine coupling of 
the nuclear spin to the local moment obtained directly 
from the shift of the nucleus by K = Axc/ij-b- Therefore 
{TiTK)~^ is a direct measure of t^. as reported for the 
^Li nuclei in in Fig. 27. They are fully compatible with 
those obtained from the ^^Y n.n. Ti data which were 
found identical for Li and Zn substitutions. 

The low energy excitations detected by inelastic neu- 
tron scattering at (tt, tt) in the normal state of under- 
doped YBCO:Zn have energy widths which agree per- 
fectly with h/Ts (Sidis et ai, 2000, 1996). This quantita- 
tive consistency among confirms that the dynamics is not 
that of a purely atomic local moment, but concerns the 
extended staggered AF state, which is the assumption 
made above in Eq. (29). . 

The Ts values are so short that the expected ESR 
linewidths are too large to be observable with standard 
lOGHz ESR, which explains the absence of any such di- 
rect observation in cuprates and implies that the ESR sig- 
nals detected in the early LSCOiZn experiments (Finkel- 
stein et at, 1990) were not linked with the Zn induced 
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FIG. 27 Electronic relaxation rate obtained from T\TK, 
as measured for the ^Li NMR for different hole dopings and 
Li contents ( empty symbols 1%, full 2%) in YBaCuO. After 
Ref. (MacFarlane et al, 2000). 

local moments (Mahajan et ai, 2000). 

The other evidence which can be drawn from the data 
of Fig. 27, is that 1/ts has a behavior quite similar to 
that found for 1/xc, as both increase quasilinearly with 
T for the underdoped cases and are less T dependent for 
the optimally doped case. In fact, a nearly perfect linear 
variation of (TiTK)^^ versus K applies, with a slope 
S which corresponds to the universal Korringa constant 
S — AwkB/fi X (/x„//Xe)^. This provides direct evidence 
that Xc a-nd r~^are governed by the same energy scale. 
In other words, the applicability of a Curie- Weiss law 
for Xc, with a Weiss temperature Q is associated with 
a limiting low T behavior t^^ — ksQ/h for the local 
moment fluctuation rate. This is strikingly reminiscent 
of the Kondo effect for a local moment in a noble metal 
host, the Kondo energy being then the width of the res- 
onant state which governs at low T the fluctuations in 
the Kondo state. These results will be discussed in more 
detail in Section V. 



D. Transport properties 

Transport properties are among the most unusual 
properties of the cuprates, and may be strongly influ- 
enced by the disorder-induced magnetic phenomena we 
have been discussing. Scattering of the charge carriers 
from inelastic magnetic fluctuations is often assumed to 
be responsible for, e.g. the strange metal behavior ex- 
emplified by the linear T dependence of the resistivity 
above Tc, although there is no complete theory available. 
The influence of impurities and defects on the transport 
properties should provide information on the properties 
of the carriers themselves, and on the coupling between 
charge and spin degrees of freedom. 

In contrast with simple metals where experiments 
probing impurity effects can be performed down to low 
T, there are severe limitations in the present case, owing 
to onset of superconductivity at a high Tc- On the other 




FIG. 28 (a) Temperature dependence of the in plane resistiv- 
ity averaged over a and b for a pure single crystal of YBCO7 
and for various Zn concentrations in substituted single crys- 
tals . (b) T variation of cot 9h, where 6h is the Hall angle, 
for the same samples (Chien et al, 1991). 

hand, the suppression of Tc by defects is an additional 
element of characterization available which serves as a 
reference in the experiments although it does not mea- 
sure a priori only the properties of the normal metallic 
state. 

Let us recall, first, that reliable in-plane transport data 
require small high quality single crystals, and that the 
synthesis of impurity-substituted single crystals does not 
necessarily benefit from the experience gained with ce- 
ramics, which are reacted in quite different thermal con- 
ditions. Thus, the impurity contents, fractional site occu- 
pancies and actual carrier content of single crystals which 
are quoted in the literature are not always accurately 
determined. As already pointed out in Section IV. B. 2 
above, electron irradiations are then quite helpful, as one 
can span a series of defect concentrations using only one 
single crystal by increasing the irradiation dose progres- 
sively. 

We shall first discuss below the basic transport mea- 
surements which allow one to demonstrate the simple 
characteristics of the scattering induced by specific de- 
fects: residual resistivities. Hall effect, thermopower, 
which will be correlated with the induced reduction of 
Tc- This will then allow us to discuss more sophisticated 
effects such as low T resistivity upturns, and magnetore- 
sistance. 



1. High T in-plane transport and scattering rates 

An important question in transport processes is the 
applicability of Mathiessen's rule, which in classical met- 
als states that the impurity scattering rate adds incoher- 
ently with the inelastic scattering rate in the pure ma- 
terial, usually due to phonons. In the cuprates, the ap- 
proximate additivity of the impurity scattering and that 
by magnetic excitations of the pure samples was pointed 
out for Zn substitutions in YBCO7 as shown in Fig. 28 
(Chien et ai, 1991). 
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FIG. 29 Resistivity versus T for a single crystal of YBCOe.e 
irradiated at 20K by 2.5 MeV electrons for increasing irradia- 
tion doses. The resistivity curves are seen to display at high T 
a parallel vertical shift associated with the increase of the im- 
purity induced residual resistivity. Therefore the defect con- 
tent does not modify the high T inflection point, linked with 
the pseudogap T* . For high defect contents, when Tc is suf- 
ficiently reduced, low T upturns are evident, and ultimately 
the resistivity curve diverges exponentially (RuUier-Albenque 
et al, 2001). 



Data taken later with electron irradiation on a single 
crystal with increasing defect content are even more con- 
vincing, as can be seen in Fig. 29. At high T, well above 
Tc the p{T) curves parallel each other, which establishes 
that the hole content is not modified by the defects in- 
troduced by electron irradiation and that a residual re- 
sistivity term simply adds to the T dependent part of the 
scattering. However, at low T and for high enough defect 
content, low T upturns of the resistivity appear. This 
phenomenon, to be discussed later, has often been inter- 
preted in terms of a metal insulator transition (MIT). 

An estimation of the residual resistivity due to the im- 
purities can be obtained either by extrapolation of the 
high T behavior to T = 0, or in the best cases as that 
displayed in Fig. 29 by measuring the parallel shift of 
the p{T) curves at high T. The residual resistivity per 
plane Ap2D is a measure of the scattering rate of the 
conduction electrons by the impurity potential. In clas- 
sical metals, the analysis of the scattering is done in terms 
of phase shifts. Such analyses are very successful as long 
as the scatterers are atomic scale defects and if one takes 
into account a decomposition of the Bloch wave functions 
into partial waves of different symmetries s, p, d... which 
undergo different phase shifts. In particular cases, as for 
transition metal impurities, one scattering channel dom- 
inates and a single phase shift is sufficient to characterize 
the scattering. While such a simple limiting situation has 
no a priori reason to be valid for impurities in correlated 
electrons systems, the data have mostly been analyzed in 
cuprates by assuming an s wave scattering with a unique 



Ap2D = -IT — sin S. 



(30) 



where Ud and n are respectively the impurity and hole 
concentrations. 

In the case of Zn substitutions the residual resistivity 
associated with the Zn scatterers yields a large value of 
the scattering rate, with sin'^ (5 sa 1/3 if Eq. (30) above 
is applied (Chien et al, 1991). This is more than two 
orders of magnitude larger than that of the Zn impu- 
rity in copper. This unexpected large scattering by a ho- 
movalent impurity is a striking illustration of the strange 
metal properties of the cuprates which may only be ex- 
plained by the correlated nature of the electron states. 
The residual resistivity associated with Cu and O vacan- 
cies in irradiated samples is large as well (Legris et al., 
1993). 

The Hall effect was also found to display an anoma- 
lously large T dependence in the cuprates, as was pointed 
out in (Chien et al., 1991). By comparing pure and Zn 
substituted samples, it was found that the cotangent of 
the Hall angle varies as 



cotOH = Pxx/Pxy = ana + /3T^, 



(31) 



as can be seen on Fig. 28. This supports the idea ad- 
vanced by Anderson (Anderson, 1991) that cot Oh reflects 
a second scattering rate , with a variation in the 
pure samples, which is independent of the resistivity scat- 
tering rate r^^. The first term in Eq. (31) implies that 
the Zn impurities contribute independently to the two 
scattering channels. 



2. Variation of the high T transport properties with hole 
doping. 

These experiments have been extended to the under- 
doped regime in order to try to quantify the variation 
of the scattering rates with hole doping. As seen in Fig. 
29 the in-plane resistivity of underdoped pure samples 
acquires a nonlinear T variation which signals a sup- 
pression of the scattering when the pseudogap opens at 
T*. Experiments on substituted YBCO or LSCO sin- 
gle crystals with Zn (Fukuzumi et al., 1996; Mizuhashi 
et al., 1995), or on electron irradiated samples (RuUier- 
Albenque et al., 2000) demonstrate that this T depen- 
dence is not altered by disorder, as can be seen in Fig. 
29, which indicates again that the hole doping and T* 
are not modified by these defects. 



a. Residual resistivity The applicability of Eq. (30) for 
different hole dopings has been checked as well. As we 
have shown above, determinations of Ap2D^ nd and n 
are done with an accuracy hardly better than 20%, so 
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that estimates of sin^ S are possible only within a fac- 
tor of 2. However (Fukuzumi et ai, 1996) have assumed 
that, within the approach described above, the scattering 
rate is large enough to correspond to unitary scattering 
{6 « 7''/2) whatever the hole doping. They suggested 
then that one should identify the number of carriers n 
with the number of doped holes for the underdoped 
samples and with I — rih for overdoped ones, to get a 
consistent explanation of the data. 

Alternatively, experiments using point defects created 
by electron irradiation allow more accurate determina- 
tions of Ap2D, as shown in Fig. 29. Furthermore, in 
YBCOe-i-x samples, the concentration Ud of defects cre- 
ated in situ at low T is independent of oxygen content 
and solely determined by the electron fluence. The ac- 
curate comparisons of Ap2_D for different hole contents 
show then that unitary scattering cannot apply for both 
YBCOe.e and YBCO7 (RuUier-Albenque et ai, 2000). 
Rather, keeping n = Uh, one finds that the scattering 
is stronger for the underdoped case, which suggests that 
the scattering rate has some relation with the range ^imp 
of the magnetic perturbation induced by spinless impu- 
rities. 



b. ATc and residual resistivity. As accurate determina- 
tions of the scattering rate are not easy to obtain directly 
from resistivity data, an alternative approach is to seek 
information from the decrease of Tc induced by the scat- 
tering centers. Although this is rather a measure of the 
pairbreaking effect associated with the defects, it gives us 
information on the scattering rates within a d-wave sce- 
nario for the superconducting state, which is supported 
by many experiments in cuprates, as will be discussed in 
Section VI. In such a case, within a BCS theory of a 
d-wave superconductor, the initial reduction of Tc is de- 
termined solely by the normal state scattering rate (see 
Sec. VI.C), 



AT, 



TT H 

Aks t' 



(32) 



Here the scattering rate should be given directly by the 
residual resistivity as = ne'^Ap2D/'m*. Thus, as long 
as the effective mass of the carriers is not sample depen- 
dent, ATc and n Ap2D should be independent measures 
of the scattering rate, and therefore scale with one an- 
other, regardless of the concentration and microscopic 
structure of the point defects. This universal scaling has 
been found to apply perfectly for Zn and electron irradi- 
ation defects in many cuprates (RuUier-Albenque et al., 
2000) (see Fig. 30) provided that the number of carriers 
is taken as n = Uh. 

The fact that dTc/dud generally decreases with increas- 
ing nji is therefore mainly due to a reduction of when 
electron correlations become less relevant. In that sense, 
one does expect to eventually recover a weak scattering 
for Zn impurities in the overdoped Fermi liquid limit. 
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FIG. 30 The decrease of Tc is plotted versus nhAp2D for var- 
ious families of cuprates, for Zn impurities or electron irradi- 
ation defects. The universal relation which applies whatever 
the hole doping rih, except for the LSCO family, establishes 
that the Tc depression is entirely governed by the scattering 
rate (RuUier-Albenque et ai, 2000). 



This is qualitatively supported by Tc data taken on Bi- 
2212 for various Zn and doping contents (Kluge et ai, 
1995) but has not been studied quantitatively in the far 
overdoped regime, because Zn substitutions in overdoped 
systems like T12201 have not been successfully performed 
so far to our knowledge. Direct determination of the scat- 
tering rate in such systems cannot be done by electron 
irradiation as the actual number of defects cannot be de- 
termined accurately either. 

Finally, this universal scaling of ATc and n A/92D does 
not apply for the LSCO system. This is particularly clear 
for Zn substitutions for which the scaling is not even valid 
at optimal doping. Fractional homovalent substitution of 
Nd on the La sites at optimal doping also yields drastic 
reductions of Tc(Fujita et ai, 2005). But this specific sys- 
tem is well known to favor stripe ordering in the Cu02 
plane ( the La i.eNd o.4o substituted samples exhibiting 
a fully ordered stripe structure even at optimal doping 
(Tranquada et al, 1997)). So, quite generally, point de- 
fects appear to favor the existence of structural and /or 
electronic instabilities in LSCO, which may explain why 
this compound behaves differently from the other HTSC 
materials. 

There are several aspects of the observed Tc suppres- 
sion which are poorly understood at this writing. At op- 
timal doping, the measured initial suppression has been 
compared to Abrikosov-Gorkov (AG) pairbreaking the- 
ory and deviations explained by anisotropic scattering 
effects (Balatsky et a/., 2006) (see however S. Graser et 
al.(Graser et ai, 2007)). In electron irradiation exper- 
iments the concentration of metastable defects created 
is not limited by the equilibrium thermodynamics which 
governs the solubility limits of impurities in a high T solid 
state synthesis. This allows to follow the suppresion of 
Tc until Tc vanishes, and the initial linear slope continues 
and eventually curves up , in contrast to the AG theory 
for which it would curve down (RuUier-Albenque et at, 
2003). This suggests the potential importance of in- 



33 



eluding phase fluetuations when the superfluid stiffness 
becomes small. 

Finally, let us recall that, as we emphasized in Sec- 
tion IV.B.l, the dopant and intrinsic structural disor- 
ders certainly play a large role in controlling the actual 
observed Tc values in many "lower Tc cupratc families" 
(Bobroff et ai, 1997a; Eisaki et ai, 2004; Fujita et al, 
2005). Obviously, if the situation for in plane disorder 
seems under experimental control, the influence of other 
types of disorders such as those discussed above in the 
case of LSCO are only qualitatively understood (Graser 
et al., 2007). We believe that quantifying the correlation 
between the decrease of Tc induced by various types of 
disorder, the associated residual resistivities and other 
physical properties as well still requires further exper- 
imental studies, which will require overcoming difficult 
materials challenges. 



c. Hall effect and thermoelectric power Hall effect mea- 
surements have been extended as well to the overdoped 
and underdopcd regime (Abe et al, 1999) and the analy- 
sis performed for optimal doping seems to apply as well, 
which highlights the existence of two scattering rates. 
However quantitative comparisons between samples are 
even harder to perform than for the residual resistivity, 
so that the data did not allow the determination of the 
physical origin of the actual difference between the two 
scattering rates. The modifications of the specific heat 
induced by Zn or Ni substitutions has been investigated, 
but the normal state electronic contributions of Zn or 
Ni could not be separated from the dominant phonon 
contribution (Loram et ai, 1990a; Mirza et ai, 1994). 
(Tallon et ai, 1995) have also found that the electronic 
thermoelectric power (TEP) is independent of Zn con- 
tent at room T, which is an independent indication that 
hole doping is not modified, as the TEP value at room 
T is found to vary markedly with hole doping in the pure 
cuprates (Obertelli et al, 1992). 



3. Low T upturns of the resistivity 

Since the parent AF state of the cuprates is insulating, 
a generic phase diagram of the cuprates should display a 
MIT with decreasing hole or electron doping. Of course 
the introduction of in plane impurities in a metallic state 
of these 2D systems is also expected to yield a MIT as- 
sociated with localization effects which should occur in 
any case if the residual resistivity exceeds the Mott-Ioffe- 
Regel limit. Both in pure and substituted cuprates, low 
T upturns of the resistivity as T is decreased have often 
been observed and have been usually associated with the 
occurrence of a MIT. 

Due to the occurrence of the SC state, the most ef- 
fective way to explore these effects is to suppress SC 
in large applied fields. It has for instance been found 
that LSCO exhibits at optimal doping a — logT increase 



of P2D in a 60 Tesla applied field (Ando et al, 1995). 
This was initially taken as an indication that the MIT 
already occurs for optimally doped cuprates, but it has 
been shown since that resistivity upturns only occur in 
the underdoped regime in BiLa2201 (Ono et al, 2000), 
suggesting that the location of the MIT is greatly influ- 
enced by the intrinsic disorder of the considered family 
as demonstrated in (Rullier-Albenque et al., 2007a). 

In the cleanest cuprates such as YBCOe.e or T12201, 
the influence of controlled defects on the conductivity can 
be studied more directly. The possibility to induce a MIT 
by introducing enough Zn in the Cu02 planes of YBCO 
was also noted in the early experiments (Fukuzumi et al., 
1996). As seen in Figure 29, resistivity upturns are also 
easily detected for electron irradiation defects. In the 
overdoped regime of T12201, the upturns of p2D only ap- 
pear for large defect content, and the magnitude of the 
— logT terms scale with , as expected for weak lo- 
calization contributions (Rullier-Albenque et al, 2001), 
which probe coherence effects in multiple impurity scat- 
tering (Lee and Ramakrishnan, 1985), in good metallic 
Fermi liquids. Therefore these experiments are consis- 
tent with the expectation that the overdoped cuprates 
are good Fermi liquid metals. 

For Zn substituted YBCOe.e the upturns of p2D were 
also detected by partly suppressing SC with an applied 
field (Segawa and Ando, 1999) and it was concluded 
that they could not be attributed to weak localization 
effects. The data accuracy for electron irradiation de- 
fects in YBCOg.e allowed the detection of upturns in p2D 
above Tc, even in zero field. A careful subtraction of the 
resistivity of the sample before irradiation allowed one 
to demonstrate that the resistivity upturns scale initially 
linearly with na so that these upturns are associated with 
single impurity scattering (Rullier-Albenque et al., 2001). 
It is therefore quite tempting to correlate these experi- 
mental evidences with the Kondo-like behavior observed 
for the impurity induced magnetism. 

If these resistivity upturns were solely explained in 
terms of spin-fiip scattering, they should be suppressed 
by a field. However, prior to any conclusion on the nor- 
mal state properties in presence of impurities, one needs 
to ensure that SC is fully suppressed by the field. In- 
deed a large SC fiuctuation regime above Tc has been 
initially detected by Ncrnst and diamagnetic measure- 
ments (Wang et al., 2006). In the presence of disorder, 
this regime has been shown to survive even after a large 
reduction of the 3D Tc (Rullier-Albenque et al., 2007a, 
2006), a fact which remains to be understood. This dif- 
ficulty to separate out the various contributions to the 
normal state scattering properties constitutes a serious 
complication when one studies the dilute impurity limit. 

For instance, a positive "magnetoresitance" has been 
reported initially in Zn substituted YBCO (Segawa and 
Ando, 1999). The survival of the SC fluctuation contri- 
bution to the conductivity which only saturates at very 
large applied fields could explain this behavior (Rullier- 
Albenque et al., 2007b). By contrast, in underdoped 
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Zn substituted BiLa2201, for which Tc could be reduced 
to zero, a negative magnetoresistance has been detected 
(Hanaki et ai, 2001). Recent experiments for irradia- 
tion defects, where high fields suppress superconductiv- 
ity, show evidence that the upturns saturate at low T 
(Rullier-Albenque et al., 2007a), also consistent with the 
Kondo analogy. 

Theoretical arguments allowing to establish the oc- 
currence of the local moment and to reproduce at least 
partly the observed behavior will be given in Section V 
below. The arguments for and against this analogy with 
Kondo effect will be discussed in some detail there. 

For increasing defect content P2d{T) becomes more 
singular than — logT and diverges at low T, signalling 
a real MIT in which localization of the carriers occurs. 
Optical experiments on Zn substituted cuprates do reveal 
a reduction of the Drude contribution to the conductiv- 
ity and the appearance of a high frequency absorption 
mode for large Zn content which have been attributed to 
such localized carriers (D.N.Basov et al, 1998). However 
this type of experiment does not give information on the 
length scale of these localized states. 



E. Significance for normal state physics 

We have described in some detail in the previous sec- 
tion the available experimental information on the influ- 
ence of impurities on the magnetic and transport prop- 
erties in the normal state of the cuprates. In this section 
we shall try to collect them and to present the main con- 
chisions which can be drawn from the qualitative aspects 
of the data. We shall try to focus on the implications for 
the physical properties of the pure system, and on the 
questions which arise then for the theory of this normal 
state. 



1. What about the phase diagram? 

Let us first return to the properties of the pure cuprates 
and summarize the implications of the impurity investi- 
gations on the phase diagram. In the low doping range 
we have seen that the controlled disorder associated 
with impurity substitutions or irradiation defects affects 
the phase diagram and expands the SG regime. So in 
this range the phase diagram might be dependent as 
well on the existing intrinsic disorder in the "lower Tc 
cuprate" families. On the other hand the pseudogap 
crossover temperature T* is unaffected by disorder, which 
explains why T* is nearly identical for most families of 
cuprates whatever their optimum Tc (Williams et al., 
1997) . However we have recalled why the impurity ex- 
periments do not allow one so far to extrapolate the be- 
havior of T* at, and above optimal doping. 

The magnetic behavior of correlated systems with im- 
purities changes abruptly near optimal doping (Fig. 25) . 
It is then tempting to associate this evolution with a 



crossover line between the strange metal and the corre- 
lated Fermi liquid regimes. NMR experiments yield a 
characteristic energy scale O for the magnetic suscepti- 
bility which would be associated with such a crossover. 6 
increases markedly as one overdopes the system and, for 
temperatures far below Q/ks, it appears that ^imp be- 
comes T independent. Better experimental accuracy is 
required before a generic phase diagram can been drawn 
using the above criteria, but the overall trend of the data 
is not inconsistent with this picture. The fate of 6 in the 
highly overdoped limit when the system should evolve 
toward a standard Fermi liquid is still an open question. 



2. Importance of electronic correlations 

Another point which has been clearly established is 
that, as in spin chains, a nonmagnetic impurity induces 
a local moment in its vicinity. This moment results from 
the sum of the staggered response of the near neighbor Cu 
sites which extends over a distance which can be charac- 
terized by a single spatial length scale ^,;,„p. This antifer- 
romagnetically correlated object is not static but highly 
dynamic, and a macroscopic magnetization can only be 
detected in a an applied field, or with a probe with a 
higher time resolution than the fluctuation rate of the 
correlated cloud. These isolated local moments cannot 
be detected by zero field /zSR experiments, which only 
probe static moments. 

These observations of magnetic responses of a metallic 
state to a nonmagnetic impurity can only occur in (and 
are therefore indisputable evidences for) the presence of 
electronic correlations in the cuprates. Similarly the large 
scattering cross section of conduction electrons by these 
nonmagnetic impurities cannot occur in a simple broad s- 
p band. It results here from the inability of the Cu3d-02p 
hybridized quasiparticles to penetrate the Zn or Li sites 
which translates into an extended scattering potential 
caused by correlations. . This effect is as well completely 
driven by the magnetic correlations in this narrow band. 

Although the importance of magnetic correlations is 
revealed by some properties of the pure systems such as 
the peak of the dynamic susceptibility at the AF wave 
vector, these experiments on impurities allow one to bet- 
ter follow the evolution with hole doping of the magnetic 
correlations and to measure both the magnitude and the 
range ^imp of the induced staggered magnetization. In 
addition, they have provided evidence that the correlated 
state is still present in optimal and lightly overdoped sys- 
tems, although with a ^imp which only exceeds slightly 
the unit cell length, and with a weak susceptibility. 

Finally static moments are only seen for large enough 
defect content, when interactions between moments pins 
a disordered spin glass like state at low T. Also in the spe- 
cific LSCO system for which the intrinsic order can yield 
the formation of stripes, impurities and stripe structures 
might interact and yield then frozen magnetic states. But 
such a situation is not encountered in cleaner cuprates. 
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3. Influence of the charge degrees of freedom on the magnetic 

effects 

In regular metallic systems, such as noble metals, a 
charged impurity induces an oscillatory Friedel charge 
oscillation which constitutes the response of the system, 
and combines the screening of the charge and the associ- 
ated spin density oscillations. In cuprates, where spin and 
charge might be distinct degrees of freedom, one might 
imagine that different length scales would characterize 
the screening and the staggered magnetism. It is some- 
times considered that a charge depletion might occur in 
the vicinity of the Zn^+ with a corresponding increase 
of the magnetic correlation in this range. Conversely if 
Li+ was to bind a hole on its n.n. oxygens, that would 
prevent a Zhang Rice singlet from visiting the Cu sites 
n.n. to Li+ and would yield quite different magnetic ef- 
fects on these sites than near Zn^+ . However the absence 
of microscopic difference over several lattice spacings for 
Li"*" and Zn^^ is a strong evidence that such effects are 
unimportant, and that the spinless character of the im- 
purities in the correlated host is the major perturbation 
which yields the appearance of a local moment. 

The scattering of the carriers also does not much de- 
pend significantly on the impurity charge, as seen from 
the analogous Tc reductions by Li+ and Zn^+. Although 
the occurrence of superconductivity severely limits the 
analysis of the transport data, we have seen that the 
scattering by spinless impurities is not always unitary, 
but that the scattering cross section increases from the 
overdoped to the underdoped regime. 



4. Quantitative information on the physical properties of the 
pure compounds 

The length scale ^imp is a signature of the response of 
the system, which is apparently only weakly connected to 
the properties of the impurity itself, as even a magnetic 
Ni induces a staggered magnetism which has a similar 
spatial extent. However the magnitude of the magnetic 
response still depends on the nature of the perturbation. 
Then ^imp is likely related to the magnetic correlation 
length ^ of the pure cuprate, reflecting the q dependence 
of x(q;'^)- In the spin 1 chains, the occurrence of a sin- 
gle degree of freedom (the spin) left obviously no other 
alternative than ^imp = ^, as confirmed by theory and 
experiment. As the impurity charge has little influence 
on ^imp in the cuprates, it is quite natural to consider 
that ^imp is at least a lower estimate of ^. 

The determination of ^ in the pure cuprates has been 
at the source of a longstanding controversy. Inelastic 
neutron scattering experiments (INS) reveal that peak 
widths of x(q,ti)) at the AF wave vector q^F- are T in- 
dependent and, if interpreted as a measure of ^"^would 
correspond to ^ w 2.5a and ^ Rs a respectively in 
the underdoped and optimally doped YBCO. In the 
slightly overdoped YBCO7 no INS peak has been de- 



tected (Bourgcs, 1998). However the INS experimental 
q width might be broadened by incommensurability ef- 
fects. Also taking into account the actual band struc- 
ture could yield a different shape of the q dependence of 
x(q,w) thereby increasing the q width of the AF peak 
with respect to ^"^ (Markiewicz, 2003, 2006). Another 
approach to determine ^ has been to analyze the longitu- 
dinal and transverse relaxation times Ti and T2G of the 
^■'Cu NMR within a phenomenological approach, assum- 
ing that xCl)"^) has a Lorentzian shape peaked at q^F 
(Millis et al., 1990). For YBCO7 this analysis gives an es- 
timate of ^ « 2a at room T, with a slight increase down 
to Tc- The extension of such a phenomenological analysis 
to the underdoped compounds in presence of the pseudo- 
gap is highly questionable and totally inconclusive (Goto 
et al, 2002). 



In summary, the length scale S,imp{T) deduced from 
the impurity experiments agrees quantitatively with ^(T) 
obtained from the analysis of the Ti and T2G ^^Cu NMR 
data in the optimally doped compound, and is always 
larger than the determinations done by INS experiments. 
So all this strongly suggests that ^imp , which increases 
markedly at low T and in the pseudogap phase, might be 
one of the best estimates of ^ of the pure system. 



We have also shown that the incidence of impurity 
scattering on transport properties gives information on 
the properties of the pure material. For instance, the 
universality of the relation between ATc and H/t induced 
by impurities can be taken as an evidence that the carrier 
density in the pure compound is indeed that of the doped 
holes. Furthermore, the Hall effect data in cuprates with 
impurity substitutions has allowed one to ascertain that 
two independent relaxation rates characterize the trans- 
port in these 2D electronic systems and that impurities 
contribute independently to both. 



Impurity experiments have revealed many important 
aspects of these correlated magnetic systems, including 
the importance of the electronic correlations and mag- 
netic effects in the perturbations induced by in-plane 
point defects in these systems. We have stressed the need 
to produce, if possible, systems free of intrinsic disorder 
to avoid non generic conclusions about the phase dia- 
gram. The actual properties of the specific response to 
isolated spinless impurities have also been found to re- 
veal quantities which are believed to be specific to the 
pure systems. Any theory aiming to describe the physi- 
cal properties of the high Tc superconductors should be 
able to explain both the existence of the local moments 
induced by spinless defects, the residual resistivities asso- 
ciated with these, the magnitude of the correlation length 
^imp, the energy scale 9 and the occurrence of the Kondo 
like behavior in magnetic and transport properties. 
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V. THEORIES OF IMPURITIES IN CORRELATED 

HOSTS 

The theory of moment formation in uncorrelated met- 
als was begun by Friedel(Friedel and Cimento, 1958) and 
Anderson (Anderson, 1961), who pointed out that a suf- 
ficiently strong local Coulomb interaction U could effec- 
tively project out of the Hilbert space corresponding to 
a weakly hybridized atomic level the state with double 
(spin-up and spin-down) occupation. This leaves in the 
effective low-energy Hamiltonian for a metal with a sin- 
gle impurity an essentially localized spin- 1/2 magnetic 
degree of freedom. The subsequent screening of this mo- 
ment by the surrounding noninteracUng electron gas as 
temperatures or energies are further reduced is a remark- 
able many-body phenomenon known as the Kondo effect, 
the subject of intense theoretical activity in the 1960s and 
70s. This literature and many further developments have 
been nicely reviewed by Hewson (Hewson, 1993). 

While wc will not discuss the Kondo problem in a non- 
interacting metal in detail here, it is important to note 
that the energy scales for the processes of magnetic mo- 
ment formation and for moment screening by the elec- 
tron gas are typically different by orders of magnitude, 
and may be discussed independently. By contrast, the 
physics of vacancies and other impurity states in strongly 
interacting chains and ladders, as discussed in Sec. Ill, 
suggests that the two phenomena of "moment forma- 
tion" (impurity-induced magnetization) and screening in- 
evitably occur together, even for weak interactions in the 
ID or quasi-lD host material. This naturally raises the 
question of the effect of a nominally nonmagnetic poten- 
tial in an interacting host in two or higher dimensions. 
In these cases, special exact analytical methods available 
to treat the single-impurity Kondo problem, or nonmag- 
netic impurities in interacting chains or ladders, are not 
applicable. However, many authors have attempted to 
treat the homogeneous strongly correlated electron host 
system using the approximate methods discussed briefly 
in Sec. II, and subsequently added an impurity-type per- 
turbation. In this section we will try to identify the in- 
gredients of a minimal tractable model which can be used 
to study the impurity in a correlated host. 



of d-band states at the pure Pd Fermi level. The Ni site 
is then simply assumed to have its own, larger Hubbard 
interaction for doubly occupied sites, Uq + AU, which is 
above the threshold to create a magnetic state according 
to the Stoner model. The spatial extent or character of 
the local magnetic state does not appear to have been 
studied extensively, as the concern was primarily to ac- 
count for bulk thermodynamic and transport properties. 

Questions regarding interaction-induced magnetic mo- 
ments arose again in studies of the metal-insulator tran- 
sition in P-doped Si and similar systems (for a review, 
see (Lee and Ramakrishnan, 1985)). Scaling treatments 
of the disordered interacting metal (Finkel'shtein, 1983) 
near the transition included interactions effectively to all 
orders, but disorder only perturbativcly, to leading order 
in (e/T-r)"^. Within such theories, the spin susceptibility 
was found to be strongly renormalized near the transi- 
tion (Castellani et at, 1986), suggesting the formation 
of clusters of magnetic moments with anomalously slow 
dynamics. This notion was confirmed generally by ex- 
periment (AUoul and Dellouve, 1987; Paalanen et al., 
1986), but results on the metallic side were found to 
be more consistent with a picture of weakly interact- 
ing diffusive quasiparticles and random local moments. 
First-principles theories of moment formation in such 
random systems were then investigated, in which inter- 
actions were treated perturbativcly (Milovanovic et al., 
1989; Paalanen et al., 1988; Sachdev, 1989), and dis- 
order treated exactly, within numerical finite-size stud- 
ies. A Hartree-Fock factorization of the interaction term, 
equivalent to that performed by (Anderson, 1961) for 
the single-impurity model, leads to a an effective single- 
particle problem with random hopping. For sufficiently 
large on-site Coulomb interaction U, it was found that 
some of the magnetic susceptibility eigenvalues turned 
negative, signalling an instability towards local moment 
formation. (Lakner et al, 1994) then pointed out that 
qualitatively very similar results could again be obtained 
by studying the isolated scattering potential in an oth- 
erwise homogeneous interacting medium. These early 
works on the Si:P problem also do not appear to have 
been concerned with the details of the magnetic state 
thereby created. 



A. Early results 

Early interest in the problem of impurity induced mag- 
netic moments in interacting systems arose in the context 
of studies of transition metal alloys such as Pd:Ni. In 
the dilute Ni limit, the appropriate problem is to study a 
single Ni embedded in the nearly ferromagnetic Pd host 
metal. A good deal of progress was obtained using so- 
called "local paramagnon theories" (Lederer and Mills, 
1968), in which the nearly magnetic host was treated 
within Stoner theory, whereby the uniform susceptibility 
was enhanced by the factor [1 — Uf)Po]~^, where Uo is the 
Hubbard interaction on Pd sites, and po the bare density 



B. Weak coupling approaches 

In these early works, correlations in the host material 
were modelled with a Hubbard Hamiltonian (Eq. 1), 
and moment formation was found to be enhanced with 
increasing repulsive U. The basic physics can be obtained 
via the usual mean field factorization: 



(33) 
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where is the on-site impurity strength, im is the impu- 
rity site, jjL is the chemical potential, and Tii^; = {c\^Cjcr)- 
The local magnetization in this decoupling where spin 
rotation invariance has been explicitly broken is simply 
Mj = (cJ|Ci-|- — c\^Cii). If the impurity potential on site 
im is taken proportional to the local occupation of op- 
posite spin CT, = e«iCT, this model is identical to the 
local paramagnon models at the mean field level, and can 
be used to describe a magnetic impurity. On the other 
hand, if Cq- = e independent of spin, the model describes 
a short-ranged screened Coulomb potential which may 
be more appropriate for an impurity like Zn, which has a 
closed shell. This latter case manifests no magnetic mo- 
ment in zero applied magnetic field if U is below the crit- 
ical value for long-ranged "background" antifcrromag- 
netic order, consistent with the experiments described 
in Section IV which one wishes to explain. This is the 
approach adopted by (Bulut, 2000, 2001) and (Ohashi. 
2001, 2002) in an attempt to model Zn impurities, and 
which we outline here. Note also that Bulut analyzed 
Hamiltonian (33) for tij limited to nearest neighbor hop- 
pings, while Ohashi included next nearest neighbor hop- 
pings. 

The goal of such calculations is therefore to calculate 
the spin susceptibility x(q = 0, R) which determines the 
response of the inhomogeneous system in the presence of 
the impurity to an applied homogeneous magnetic field, 
and compare its temperature dependence to that of ex- 
periment (here the Fourier transform is with respect to 
the relative electron coordinate r — r', leaving the depen- 
dence on the center of mass coordinate R = (r -I- r')/2). 
Once this quantity is known, the Knight shift and, in 
principle, the distribution of local fields around the im- 
purity as probed by NMR can be obtained. The physical 
idea put forward by these authors is that the metallic 
state in which the impurity is embedded manifests strong 
but noncritical antiferromagnctic spin correlations maxi- 
mized near q = (tt, tt), which are enhanced locally by the 
impurity. The related problem at criticality was studied 
by (Millis et at, 2001). 

The impurity potential then couples the enhanced re- 
sponse of the system at (tt, tt) and that near q = which 
is then measured in the NMR experiment. This can be 
most easily seen by examining the susceptibility deter- 
mined within the random phase approximation (RPA) 

X(r,r') = xo(r,r') + C/^xo(r,r")x(r",r'), (34) 

r" 

where xo is the susceptibility determined in the non- 
interacting system {7 = 0, given exactly by the non- 
interacting Green's function Go(r,r') in the presence 
of a single short-range impurity described by a single- 
impurity T-matrix. In the homogeneous case x(r, r') = 
x(r — r')j the result in Fourier space is the familiar 
x(q) = Xo(<l)(l — C^Xo(q))~^, i-e. all momentum compo- 
nents are decoupled. On the other hand, in the general 
case (34), Fourier transforming with respect to the rel- 
ative coordinate cannot disentangle the various momen- 



tum components, leading to a general integral equation 
relating x(q, R) to x(q',R')- Similar results for the nu- 
merical evaluation of these equations were obtained in 
(Bulut, 2000, 2001) and (Ohashi, 2001, 2002), and are 
plotted e.g. in Fig. 31 as a function of temperature. 
While the enhancements at q = are smaller than at 
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FIG. 31 Temperature dependence of magnetic susceptibility 
in RPA at q = for different sites R in the presence of a 
strong potential Vi = 20t at (0,0) and a weaker one on nearest 
neighbor sites of strength V2. x is a weighted average of x to 
account for the fact that R may coincide with a lattice site, 
the center of a bond or of a plaquette (Ohashi, 2001). 

(tt, tt) as expected, they are still quite significant at low 
temperatures and appear to follow a ~ 1/T form down 
to the lowest temperatures on the next nearest-neighbor 
sites (1,1) when the extended impurity potential on the 
nearest neighbor sites is large. For smaller nearest neigh- 
bor coupling, the enhancement at g = weakens. At in- 
termediate values of the extended potential strength, the 
susceptibility T dependence resembles the experimental 
NMR results discussed in Section IV where in the metal- 
lic phase this temperature dependence is found to follow 
a behavior x ~ (T-f-O)"^. Note that the mean field form 
of the Hamiltonian solved, whereby the spin rotation in- 
variance is explicitly broken, precludes these theoretical 
results originating from Kondo-type physics. Indeed, to 
describe Kondo physics one needs to allow spin-flip pro- 
cesses yielding transitions between the two spin eigen- 
states locally, and these large amplitude fluctuations are 
not included within RPA. 

Using an essentially identical formalism, Bulut explic- 
itly attempted to fit the Knight shift data on Li and Zn 
impurities in optimally doped YBCO samples, with con- 
siderable success (Bulut, 2001). In his work, however, the 
potential on nearest-neighbor sites is taken to be attrac- 
tive, leading to a magnetic moment residing primarily on 
these sites. Hyperfine parameters, band parameters, and 
the Coulomb repulsion U were chosen from fits to NMR 
data on the pure system. In Fig. 32 we reproduce the 
Knight shift data of (Bobroff et al., 1999) together with 
Bulut's calculation. It is interesting to note that similar 
quality fits to data on a different impurity, Zn, were ob- 
tained with nearly identical values of the nearest neighbor 
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impurity potential, suggesting a universality of the mag- 
netic response to spinless in-plane defects as discussed 
already in Section IV. The same model also reproduces 
the staggered magnetization pattern seen by experiment 
on sites near the impurity (although a quantitative com- 
parison was not attempted). 

The physical meaning of the nearest neighbor potential 
which appears so vital to explain the impurity-induced 
enhanced magnetic correlations is not completely trans- 
parent. Band calculations in fact suggest that such im- 
purity potentials have a range of an A or less (Wang 
et al, 2005). This term may therefore play the phc- 
nomenological role of a different impurity-induced cor- 
relation effect in a more complete theory. The on-site 
impurity potential breaks translational invariance such 
that the response to a uniform magnetic field contains 
a contribution related to the staggered susceptibility of 
the pure system. This quantity is quite large, but its 
weight in q-space is small. Adding a nearest neighbor 
potential V2 boosts the weight at the antifcrromagnctic 
momentum. In addition, Monte Carlo calculations (Bu- 
lut, 2003) begin with an on-site potential only, and find 
that an effective nearest neighbor potential of roughly 
the right order to fit the Knight shift data in the RPA 
approach is dynamically generated. This phenomenon of 
creation of dynamical extended potentials by nonmag- 
netic pointlike perturbations in the presence of Coulomb 
correlations had been pointed out earlier (Poilblanc et al., 
1994a) and discussed more generally in the weak cou- 
pling context (Ziegler et ai, 1996b). A similarly en- 
hanced range of the effective impurity potential due to 
renormalized electron-phonon coupling was predicted by 
(Kulic and Oudovenko, 1997). So in this sense, a vacancy 
or strong potential in a Hubbard model may indeed pro- 
vide a minimal model to describe a Cu substituent in the 
Cu02 plane of the cuprates around optimal doping, if the 
extended effective potential is generated systematically. 
However, the RPA theory in the form given in the work 
of Ohashi and Bulut requires a phenomenological finite 
range potential. Furthermore, it is not systematically ex- 
tendable to the underdoped side of the phase diagram, 
although plausible results can be obtained by putting in 
the pseudogap in an ad hoc way (Bulut, 2001). 



C. Strong coupling approaches 

It is generally believed that the transition to the Mott 
antiferromagnet at low doping in the cuprates can be de- 
scribed theoretically only with models which invoke local 
Coulomb interactions in the strong coupling limit. His- 
torically, the problem of a nonmagnetic defect in a cor- 
related system was first proposed in the cuprate context 
by (Finkelstein et at, 1990) to account for the depression 
of superconductivity in La-based compounds (see next 
section). For the normal state case, (Poilblanc et al., 
1994a, b) were the first to calculate the electronic struc- 
ture of the t — J model (Eq.(2)) when a single spin was 
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FIG. 32 Temperature dependence of Knight shift (Bobroff 
et al, 1999) for Li impurities in YBCO within RPA approach 
in the presence of a strong potential at Vi = — lOOt at (0,0) 
and a weaker one on nearest neighbor sites of strength given 
in the figure. Prom (Bulut, 2001). 



frozen, i.e. all exchange couplings and hoppings con- 
nected to a given site on a square lattice were removed. 
They found that the rest of the system had a net mag- 
netic moment corresponding to net spin-1/2 distributed 
over the sites near the impurity. A hole added to the sys- 
tem was found to bind to the impurity site in one of four 
symmetry channels s, Px,y and d corresponding to the ir- 
reducible representations of the tetragonal group. These 
bound states occurred for couplings above a critical value 
of J of order ~ O.lt, and it was speculated that for a 
thermodynamically finite density of holes the critical ex- 
change would be enhanced due to the additional allowed 
decay channels for the bomid states. Further finite size 
numerical calculations ont—J impurities using Hubbard 
operator techniques (Odashima and Matsumoto, 1997; 
Odashima et al., 2000) indeed showed the weakening of 
bound states at finite doping, but did not investigate the 
transition to spontaneous local magnetization. 

While the numerical calculations of (Poilblanc et al., 
1994a,b) indicated a tendency to spontaneous moment 
formation around a nonmagnetic potential in a strongly 
correlated host, they left open many questions regard- 
ing the physical origin of the phenomenon, as well as its 
dependence on doping and correlation strength. Along 
with the experiments on Zn-doped cuprates reviewed in 
Section IV at about the same time, this stimulated con- 
siderable interest in the search for tractable analytical ap- 
proaches to the problem. It is worth recalling that there 
are many plausible schemes to obtain results the strong- 
coupling regime in 2D, none of which are fully controlled. 
(Gabay, 1994) studied the doping dependence of an im- 
purity moment in a slave boson approach to the t — J 
model, as discussed in Section II, and tried to relate the 
size of the moment associated with Zn in different cuprate 
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materials and for different dopings to the screening length 
of the spinon fields, particularly in the spin gap phase. He 
also predicted a local destructive effect of the impurity- 
induced magnetism on the spin gap itself. (Nagaosa and 
Ng, 1995) considered an impurity in a tt-Hux phase of 
the Heisenberg model using a slave boson representation 
and predicted the existence of a spin- 1/2 bound state. 
Similar results were obtained using a slave fermion rep- 
resentation. These authors noted explicitly the analogy 
with ID end-chain defects, and pointed out that although 
the effective potential for the spin was weaker in 2D, a 
Curie- Weiss term in the susceptibility would be the likely 
result. 

(KhaliuUin et at, 1997) and (Kihan et ai, 1999) pro- 
duced detailed calculations in a similar spirit which re- 
sulted in specific analytical predictions for impurity prop- 
erties in the spin gap state. The spin- or pseudogap 
homogeneous state was modelled by a "drone" -fermion 
(Khaliullin et al., 1997) or more conventional slave bo- 
son approach (Kilian et al, 1999) which yielded similar 
results. In both cases auxiliary bond field phases are 
chosen within a homogeneous mean field description to 
give a linear 'pseudogap" in the spinon sector of the the- 
ory, with density of states p''\uj) = \oj\/D'^, where D is 
the spinon half bandwidth, see Fig. 33b. Formally, this 
choice of mean field parameters corresponds to the flux 
phase of (Affleck and Marston, 1988). To model a non- 
magnetic impurity in the pseudogap state, a site is then 
removed from the system by means of a local potential 
term Aq/q/o, with Aq — > oo, and calculation of the im- 
purity T-matrix then yields a spinon bound state at zero 
energy Pimp{^^) = S{lv). The local susceptibility due to 
the bound state varies as 



1 <^(R) 



27r i?2 TlnD/T' 



(35) 



where the phase factor $(R) gives the spatial dependence 
of the spin polarization in a magnetic field, 



$(R) 



(36) 



with i?± = + Ry and = {R^ ± iRy)/R. This is 
a staggered pattern of spins on one sublattice only , in 
contradiction with the results of the weak coupling cal- 
culations described above, and considered by the authors 
(Kilian et al., 1999) as an artifact of the mean field the- 
ory which severely underestimates long-range antifcrro- 
magnetic correlations, as evidenced by the fact that the 
antiferromagnetic correlation length is absent from the 
expression (35). To improve the accuracy of the calcula- 
tion and include these effects, a phenomenological RPA 
summation was performed along the lines of (Morr et ai, 
1998), x(q) = X"(q)/(1 + ^qX°(q)), where xo is the ho- 
mogeneous susceptibility and Jq = 2J{cosqx+cosqy)), as 
in the theory of the nearly antiferromagnetic Fermi liquid 
(NAFL) (Monthoux and Pines, 1994). For the inhomoge- 
neous susceptibility, approximate analytical expressions 



can be obtained for the long distance behavior, 

.x(R,T).-cos(Q.R)i-i,^i§!^, (37) 

where Q =_(^,7r), and ^(T) = [Jxo(Q,r)/(l - 
4 Jxo(Q, T))] is the NAFL antiferromagnetic correla- 
tion length for the pure system,. The final result Eq.(37) 
shows the commensurate staggered behavior of the mo- 
ment as well as the (TlnT)"^ divergence of the suscep- 
tiblity, reflecting the marginal bound state found already 
in the mean field theory. 

(Kilian et al., 1999) compared their results with data 
from (Bobroff et ai, 1997b) on imderdopcd YBCO, and 
concluded that a good fit could be obtained if the tem- 
perature dependence of the correlation length in Eq.(37) 
was assumed to follow (,{T) ^ l/(a + bT), with a and b 
adjustable parameters. With these parameters fixed, the 
theory can account for the broadening of the NMR line 
with increasing Zn concentration. It is seen that, while 
not all parameters could be computed directly from mi- 
croscopic theory, a physically reasonable picture of the 
impurity in the pseudogap state is ultimately obtained. 
It should be noted however that the spatial decay of the 
local susceptibility as one moves away from the impurity 
has a different functional form in (Ouazi et al., 2004) 
compared with Eq.(37). Indeed, NMR lineshapes cannot 
be accounted for using R~'^ or R~^ spatial dependences 
for Sx- 




FIG. 33 Schematic plot of spinon density of states for a) spin- 
1 impurity; b) spin-0 impurity in model of (Kilian et al, 1999). 
Dashed line: density of states for pure flux phase. Solid line: 
with impurity. 

The strong coupling approaches discussed above as- 
sumed the existence of a pseudogap in the one-particle 
density of states, which is ultimately required for the 
formation of the spinon bound state. One would there- 
fore expect the singular magnetic response to disappear 
when the pseudogap does; indeed, this is consistent with 
the fact that no static moment is observed in zero field 
experiments on clean optimally doped cuprates and that 
the field induced magnetic response is not a simple Curie 
(1/T) form. (Section IV). However a theory of this 
type should be able to describe the paramagnetic mo- 
ments observed by NMR experiments in finite field on an 
equal footing. This was accomplished recently by (Gabay 
et ai, 2007), who considered the optimally doped state 
of cuprates within at — t' — J model (Eq. 8); the strong 
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impurity potential implies a site from which spinons and 
holons are forbidden. The Hamiltonian was treated in a 
slave boson mean field approach, following (Ubbens and 
Lee, 1992). This model remains metallic at low T and 
is thus suitable for the study of the normal state. One 
introduces the local order parameters 
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where fi and bi create spinon and holon fields, respec- 
tively, and i,j are nearest neighbors. Within such a 
theory, the spinon bandwidth = tQij + {J/2)xij is 
strongly suppressed close to half-filling i5 — > 0. When the 
effect of the impurity is accounted for, the bandwidth is 
found to decrease locally as well, as shown in Fig. 34c. 
Such effects are quite important for a quantitative treat- 
ment of local correlations and electronic structure around 
the impurity, and are missed in the weak-coupling ap- 
proaches discussed above. A semianalytic solution of the 
mean field equations reveals that the local decrease of 
the spinon bandwidth in the vicinity of the impurity has 
two effects. One is to produce an extended scattering 
potential which enhances the staggered Fourier compo- 
nent of the local magnetizations mi. The other is to in- 
crease the magnetic response, since - in a Stoner-like pic- 
ture - a larger value of J/tf(i) brings the system locally 
closer to a magnetic phase. In order to treat correctly the 
non double occupancy constraint for the holons and the 
spinons, a numerical solution of the mean field equations 
is required (Gabay et ai, 2007). It yields quite plausible 
results for both the spatial decay of the staggered mag- 
netization induced in a field, and the temperature depen- 
dence, which is found to correspond to a Curie- Weiss law 
similar to that found in experiment on overdoped sam- 
ples. Indeed, the results displayed in Fig. 34 are in rea- 
sonable agreement with their experimental counterparts, 
shown in Figs. 25,26. It is intriguing that the results of 
the slave boson mean field theory for the strong-coupling 
t — t' — J model in the spin liquid state (optimally doped 
case) are so qualitatively similar to the weak coupling re- 
sults presented in Sec.V.B. This suggests indeed that for 
the optimally doped materials a Fermi liquid approach 
captures the essential details of the physics. On the other 
hand, in the strong coupling case there is no analog of 
the requirement of an extended bare impurity potential 
to fit experimental data as found by (Bulut, 2000) and 
(Ohashi, 2001). This may indeed reflect the failure of 
the weak coupling approaches to treat the renormalized 
band structure in the vicinity of the impurity, e.g. the 
spinon bandwidth shown in Fig. 34. At present, no mi- 
croscopic approach to the problem of an impurity in a 
correlated host has been able to reproduce systematically 
the crossover between underdoped and optimally doped 
regimes observed in the cuprates, and described in Sec- 
tion IV. In principle, slave boson and similar renormal- 
ized mean field theories are a promising avenue to pursue. 
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FIG. 34 Results for fully self-consistent evaluation of mag- 
netization from slave boson equations (Gabay et al., 2007) . 
a) Normalized magnetization si on nearest-neighbor site as 
function of T for values of doping 5 from 0.14 to 0.30. Solid 
lines show fits to si = C/(r -|- 0). b) Effective moment con- 
stant C and Curie- Weiss constant Q extracted from fits in a) 
vs. 5. c) Spinon bandwidth in fully self-consistent evalua- 
tion as a function of distance from the impurity site at filling 
5 = 0.15 and T = lOOK (filled squares) and 300K (fiUed cir- 
cles). Bandwidth from semianalytic calculation at J = 0.3 
and lOOK is shown (open circles) for comparison, d) Temper- 
ature dependence of extracted theoretical correlation lengths 
^ vs. T (compare Fig. 26) . 



In addition, a phenomenological approach by (Prelovsek 
and Sega, 2004) has suggested a way to parameterize 
this crossover in terms of the doping dependent correla- 
tion length ^ of the pure system. These authors showed 
that a crossover between the Curie and Curie- Weiss be- 
havior of the local susceptibility observed in experiments 
as the system is doped is to be expected generically in 
a homogeneous doped antiferromagnet. They made the 
ansatz that the impurity suppresses only the spin cor- 
relations (S^Sj) involving the impurity site 0, with all 
others taken to correspond to the pure system. While 
this seems to contradict other microscopic analyses dis- 
cussed above, the results indeed resemble experiments 
qualitatively and suggest that the pure ^ determines the 
T-dependence of the system. 



D. Magnetic impurities in correlated systems 

One of the attractive features of the (Kilian et al, 
1999) approach is that it accounts naturally for differ- 
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cnccs observed between Zn and Ni substitutions for Cu 
in undcrdopcd YBCO. In Fig. (35), we show a schematic 
picture of an RVB type singlet state in the presence of an 
impurity of various types. In the case of a nonmagnetic 
impurity hke Zn, a spin singlet is broken while a spin- 1/2 
is freed from a given singlet bond. By contrast, if one 
models the Ni perturbation as a spin-1 object, the spin- 
1/2 on the Cu site nearest the impurity screens the Ni 
spin, leading again to a free spin-1/2 degree of freedom. 
This picture is very analogous to that which has been 
established rigorously for impurities in ID spin chains 
in Sec. III. Of course, as we have seen, more detailed 
calculations lead to the "free" spin 1/2 being spread out 
over a healing length in both cases. To model a Ni spin 
So of magnitude 1 coupled to the surrounding Cu spins- 
1/2, (Kilian et al, 1999) add an impurity potential con- 
sisting of nearest neighbor exchanges with the impurity 
So • (Js to the system. The spin 1 is decom- 
posed So = Sa -|- Sfe, where Sa,b are spin-1/2 operators, 
and the spin 1 algebra is retained by adding a ferromag- 
netic exchange term — JcSa • S(, to the Hamiltonian, then 
allowing Jg — > oo. This gives 



analysis as above is found to be 



imp 



Y^i^'sfLha + h.C.) - JeSeffao, (39) 



5a 



where 



^crifsafoa) IS the local flux phase am- 



plitude on bonds near the impurity, Sgjy is a localized 
spin-1/2 operator, and uo is the spinon spin density on 
site 0. Thus this case reduces to an impurity which lo- 





(b) 



FIG. 35 Schematic picture of spin-1/2 impurity produced in 
case of a) Ni and b) Zn. 

cally renormalizes the spinon bandwidth and adds a local 
exchange with a free spin-1/2 in a fcrmionic bath with 
pseudogap. Note that such a treatment clearly neglects 
the transverse spin fluctuations which lead to the Kondo 
eflect. Nevertheless the pole in the T- matrix is found 
to have a resonance in the pseudogap at finite energy 
= (7r/4) J'/ ln(£'/J'), leading to a contribution to 
the density of states 



1 col 



2 InD/iVK (w2 - u;^)2 + (77^2^/4 In D/wk)^ ' 

(40) 

whose form is shown schematically in Fig. (33a). The im- 
purity susceptibility after performing the identical RPA 



-7;^cos(Q.R) — 



, (41) 



where Jc(0) interpolates between D for J' = J and 
2ojk for J' <^ J. This expression is quite similar to 
the nonmagnetic case (37), but decays more rapidly, as 
~ and has a more singular l/T dependence. The 

reduction of both spin-0 and spin-1 impurity models to 
Hamiltonians with effective localized spin-1/2 exchange 
terms raises the question of the fate of the Kondo eff'ect, 
the screening of a single localized spin exchange coupled 
to a free Fermi gas, as electronic correlations increase. 
This effect was neglected in the mean field treatments 
above. Historically, the question was studied first by 
(KhahuUin and Fulde, 1995; Schork and Fulde, 1994). 
These authors began with the slave boson mean field the- 
ory of (Read and Newns, 1983) for the Anderson impu- 
rity model, which is known, in the noninteracting case, 
to correctly reproduce the low-temperature Kondo reso- 
nance in the /-electron density of states at the correct 
scale Tk- Perturbative corrections in the Hubbard U 
were then found to simply renormalizc the various Read- 
Newns parameters, including the Kondo scale Tk- 



Tk 
Tk 



— — = exp 



po-'^(l + a) 



(42) 



where a =- (3/2) In 2C/ po, ^ is the exchange of the Ander- 
son impurity with the conduction electrons in the Kondo 
limit, and po the density of conduction electron states 
at the Fermi level. Thus for the usual repulsive-f case 
the Kondo scale is enhanced due to the suppression of 
the charge degrees of freedom. Similar conclusions were 
reached by (Hofstetter et al, 2000) in a more sophisti- 
cated calculation combining dynamical mean field theory 
with numerical renormalization group methods. In addi- 
tion to the eff'ects of interactions on impurity properties, 
there are interaction-dependent renormalizations of con- 
duction electron properties by the Kondo impurity which 
are predicted to be visible in thermodynamics. (Neef 
et al, 2003). 



E. Transport 

There is currently relatively little theoretical work at- 
tempting to explain the apparent close relation between 
the disorder-induced magnetic properties of cuprate sys- 
tems and transport properties, e.g. the resistivity "up- 
turns" at low T observed in cuprate materials when they 
become very disordered or are subjected to high fields. 
Based on the general properties of gauge theory descrip- 
tions of the HTS phase diagram, (Nagaosa and Lee, 1997) 
attempted to understand the residual resistivity measure- 
ments of (Fukuzumi et al., 1996) in terms of a Boltzmann 
transport description with assumed unitarity limit scat- 
terers and carrier concentration proportional to hole dop- 
ing Uh in the underdoped cuprates, crossing over to l — rih 
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for optimal doping. Some criticisms of this picture have 
already been mentioned in Section IV. Recently, (Kon- 
tani and Ohno, 2006) calculated the local susceptibilities 
in a Hubbard model with strong impurities within FLEX- 
type approximations, and fed the resulting self-energies 
and vertex functions into a disorder-averaged transport 
calculation. The enhanced scattering due to the creation 
of local moments led in some cases to upturns in the resid- 
ual resistivity at low temperatures similar to those seen 
in experiments. While some details do not correspond to 
observation (upturns in YBCO were found to be stronger 
than in LSCO, for example), this is a promising step 
allowing one possibly to approach the "metal-insulator 
transition" from the metallic side. 



F. 2D versus ID 

Since there is to date no consensus on a theoretical 
description of the homogeneous system in the strongly 
correlated limit, further insight into the 2D case may be 
gained by considering, once again, the ID situation. The 
literature devoted to the study of impurities in doped ID 
systems is quite significant, but we discuss here a very 
small number of results which are directly relevant to 
the 2D case. In ID, a weak link, i.e a slight local mod- 
ification of the kinetic energy and of the magnetic cou- 
pling or a weak potential modelling a charge imbalance 
renormalizes to a very large perturbation as the temper- 
ature is reduced, in the case of repulsive interactions. As 
T ^ 0, the transmission in the spin and charge chan- 
nels approaches zero (Furusaki and Nagaosa, 1993; Kane 
and Fisher, 1992a,b). Applying a uniform magnetic field 
H yields a local alternating magnetic polarization with 
an amplitude proportional to H and a spatial extension 
generically of the order of the correlation length of the 
pure system. The crossover between the weak pertur- 
bation and the strong perturbation regimes is frequently 
described in Kondo language, and the characteristic tem- 
perature Tk depends exponentially on the difference be- 
tween the kinetic (or potential or coupling) energy of the 
pure system and that of the modified link (Rommer and 
Eggert, 1999, 2000). 

The impact of the impurity on transport and on mag- 
netism can be considered separately, because of spin- 
charge separation at low energies, but, in both channels 
backscattering is very strong at low temperature. The 
Luttinger liquid parameter K varies from i to 1 upon 
increasing the hole doping, showing a drastic reduction 
of the enhancement of the staggered polarization as one 
moves toward the Fermi liquid limit. A similar behavior 
is observed in numerical studies of ladder systems (White 
et al., 2002). Calculations of transport properties also 
produce low T upturns in the resistivity which reflect 
the pinning of charge density waves by the impurities. 



G. Summary 

At this stage, we may attempt to interpret some of the 

results described in Section IV in the framework of the 
theoretical models that we just discussed. Experiments 
reveal that a spinless impurity introduced in a correlated 
host produces in its vicinity a very large and spatially ex- 
tended alternating magnetic polarization, similar to ID. 
On the basis of NMR and /kSR spectra it has been demon- 
strated that this magnetization is not associated with a 
static moment. The maximum amplitude of the stag- 
gered polarization is found to be of the form C H/{T+Q), 
where C is a constant and the extension ^j^p of this po- 
larization away from the defect seems to follow a similar 
T dependence. For a fixed T, the ( spatial) extent of 
the staggered state increases as one underdopes the sys- 
tem (0^0). A finite density of impurities produces a 
magnetization which scales linearly with the concentra- 
tion. Moving to the optimally doped and, beyond, to the 
overdoped regime, both the amplitude of the alternat- 
ing polarization and its spatial extent decrease, bringing 
these gradually back in line with values expected for a 
Fermi liquid metal. 

Furthermore both NMR spin lattice relaxation and in- 
elastic neutron scattering data suggest that this entity 
should be viewed in the optimally doped case as a reso- 
nant state rather than a free moment. By this we mean 
that, in zero field, the total spin of the system in the 
ground state is S* rather than S = 1/2. The fluctua- 
tion rate h/r « ks® corresponds then to about lOOK or 
lOmeV. 

For a wide range of hole dopings, the transport scatter- 
ing rate in Cu02 materials displays a — log T upturn at 
low T . In the overdoped case, this behavior is perfectly 
explained by weak localization effects in a good metallic 
system. In the underdoped case, the upturns are much 
larger than predicted by weak localization theory, and 
may be due to impurity induced magnetism. 

The above features are quite suggestive of Kondo 
physics, with the picture of a screened moment formed 
around a nonmagnetic impurity. However, as we have 
seen, a number of these magnetic and transport signa- 
tures are qualitatively captured by the theoretical mod- 
els reviewed in the present section which ignore processes 
leading to Kondo screening. Typically, the impurity is 
represented as a strong or infinite nonmagnetic poten- 
tial, which creates a nearby cloud of staggered polariza- 
tion and net paramagnetic moment by coupling the sys- 
tem's antiferromagnetic correlations to a homogeneous 
external magnetic field. The apparent universal mag- 
netic response of any impurity substituting for Cu in the 
Cu02 plane found by NMR is strong evidence in favor of 
such models. Remarkably, treatments which start either 
from a weak-coupled Fermi liquid or from a strong cou- 
pling description of the host are able to reproduce the 
key features seen in experiments; a characteristic energy 
© is foimd, the temperature dependence of the Knight 
shift, the resonant nature of the staggered magnetic po- 
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larization, and general trends in the observed doping de- 
pendence are obtained roughly correctly. Weak coupling 
approaches have been mostly devoted to the study of the 
optimally doped regime; they appear to require tuning 
the Hubbard ?7 to a value close to an AF instabilility 
and using extended impurity potentials. In the strong 
coupling limit, mean field approximations give a qualita- 
tively correct carrier doping dependence for 6 (and for 
C, which changes only slightly with the hole content). 
A current drawback, however, is the inability of any sin- 
gle model to span the entire cuprate phase diagram and 
make predictions for the full doping dependence of the 
response. 

Further considerations suggest caution against taking 

the Kondo analogy too literally. In the Kondo/ Anderson 
picture there are actually two energy scales. The upper 
one is related to the formation of the local spin, i.e to 
a (local) bound state, and the lower is related to the 
screening of the moment by conduction electrons. In 
cuprates, moment formation and screening are driven by 
the same set of conduction electrons. Furthermore, it 
is unclear how different an extended itinerant staggered 
magnetization cloud with a net moment which is ob- 
tained when a non magnetic impurity is introduced in a 
correlated system - is, from the single localized impurity 
spin 1/2 that forms at high energies, in the Kondo case. 
In the Kondo scenario of a magnetic impurity in a no- 
ble metal host, the characteristic energy scale is set by 
the carrier density and the exchange interaction through 
6 ~ Jexp— [Jp{Ep)]~'^ . For underdoped cuprates, the 
non superconducting regime corresponds to temperatures 
T >> 6, since Q is found to approach zero; in a Kondo 
description, T >> 6 would correspond to the weak 
coupling limit J p{Ep)<< 1: yet spin dynamics, seen in 
NMR experiments, continue to show the fluctuations of 
a strongly correlated extended magnetic entity around 
the impurity. This might be due to the existence of the 
pseudogap. 

In the overdoped regime, as the physics gets closer to 
the Fermi liquid limit, one does not expect to observe 
Kondo-like signatures any longer. Indeed, a nonmagnetic 
impurity is not expected to produce a Kondo resonance in 
conventional metals. Thus, understanding the reduction 
of the amplitude of the staggered moments, of the spatial 
extent of the polarization and the evolution of C and G 
for large dopings may require a different description. 



VI. IMPURITIES IIM THE SUPERCONDUCTING STATE 
A. Overview 

The subject of impurities in the superconducting state 

of the cuprates, and the theory of impurity states in un- 
conventional superconductors generally, is already quite 
vast. We benefit from the availability of recent reviews of 
this field by (Balatsky et al., 2006) on local aspects of im- 
purity states, by (Hussey, 2002) on effects of disorder on 



bulk quasiparticlc properties, by (Hirschfeld and Atkin- 
son, 2002) on nonperturbative disorder effects in two di- 
mensions, and by (Flatte and Byers, 1999) on general 
impurity properties in superconductors. Here we review 
briefly the simplest aspects of the effects of nonmagnetic 
impurities on d-wave superconductivity before discussing 
theoretical approaches to the inclusion of electronic cor- 
relations in the superconducting state impurity problem, 
as well as relevant experiments. 

It was recognized shortly after the publication of the 
BCS theory that while nonmagnetic impurities would not 
affect thermodynamic properties of a conventional super- 
conductor (Anderson, 1959), they might do so in hypo- 
thetical p-wave pairing systems (Balian and Werthamer, 
1963). When a high density of low-energy excitations 
in superconducting heavy fermion systems, such "un- 
conventional" states were considered as possible ground 
states of compounds like UBcia and UPta. (Gorkov 
and Kalugin, 1985) and, independently, (Ucda and Rice, 

1985) pointed out that linear nodal regions of the or- 
der parameter on the 3D Fermi surface would lead to a 
nonzero residual density of states at zero energy, p(0). 
This self-consistent perturbative treatment for weak im- 
purities, averaged over disorder, was similar to that used 
by (Abrikosov and Gorkov, 1960) in their discussion of 
gapless superconductivity induced by magnetic poten- 
tials. (Pethick and Pines, 1986) observed that weak scat- 
tering models of this kind could not describe transport 
properties measured on heavy fermion materials, and 
proposed that large effective impurity potentials could be 
treated in a T-matrix approximation. (Hirschfeld et al., 

1986) and (Schmitt-Rink et al., 1986) then proposed the 
self-consistent T-matrix approximation (SCTMA) for im- 
purities of arbitrary strength, and predicted a disorder- 
induced "plateau" in p(w) at low energies. This plateau 
consists of disorder induced " impurity band" states con- 
ceptually analogous to the mid-gap localized states pre- 
dicted earlier for dirty isotropic p-wave superconductors 
by (Buchholtz and Zwicknagl, 1981). The nodal quasi- 
particles which contribute to this "impurity band" give 
rise to thermodynamic properties similar to a normal 
metal with reduced DOS. (Stamp, 1987) then investi- 
gated the character of a local single-impurity resonance 
within this model. 

In a similar way, when experiments on good samples 

of cuprate superconductors in the early 90s began in- 
dicating the existence of low-energy excitations in the 
superconducting state (Hussey, 2002), the importance of 
understanding the effect of disorder for identifying the 
symmetry of the order parameter was quickly realized. 
The effect of disorder on Tc itself was controversial early 
on. For example, irradiation with heavy ions was found 
to suppress Tc only when the scattering rate, as deduced 
from resistivity, became of order the Fermi energy (Valles 
et al., 1989), as expected for an isotropic s-wave super- 
conductor. On the other hands, impurities substituting 
for Cu in the Cu02 plane tended to depress Tc rapidly, as 
expected for an unconventional superconductor, and con- 
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trolled electron irradiation experiments confirmed this 
(Giapintzakis et al., 1994), (Legris et at, 1993) (see Sec- 
tion IV.B.2). 

On the theoretical side, the self-consistent T-matrix 
approximation (SCTMA) was employed to provide a 
qualitative explanation of the effects of disorder on d- 
wave superconductors, offering important evidence for 
the symmetry identification. These arguments belong 
in the same category as the explicitly phase sensi- 
tive measurements like the tricrystal experiments (Tsuei 
and Kirtley, 2000), as pointed out by (Borkowski and 
Hirschfeld, 1994). This is because the fact that nonmag- 
netic impurities break pairs is directly associated with 
the sign change of the d-wave order parameter (see Sec. 
VLB). Within the SCTMA, the strength of the impurity 
potential is an important parameter, and the recognition 
that planar impurities like Zn acted as strong scatter- 
ers which, from the standpoint of theoretical modelling 
could be considered as unitary scatterers, was an im- 
portant step forward. For example, observations of a 
crossover of the magnetic penetration depth from a T to 
a behavior with a small concentration of Zn impuri- 
ties were not initially midcastood to be a natural con- 
sequence of d symmetry, because it was expected that 
in such experiments Tc would be strongly suppressed 
in a d-wave system. In the unitarity limit, however, 
strong modifications of low-energy properties occur over 
an energy scale 7 (residual scattering rate or "impurity 
bandwidth") which can be much larger than the nor- 
mal state impurity scattering rate F which controls Tc 
(Hirschfeld and Goldcnfeld, 1993). A similar analysis by 
(Hotta, 1993) allowed for an explanation of the depen- 
dence of the residual T Knight shift of YBCO on 
Zn concentration. Calculations of transport properties 
within this framework also provided evidence for near- 
unitarity scattering, and helped rule out candidate ex- 
tended s-wave states (Borkowski et ai, 1995). 

While quite successful in explaining the effects of 
strong in-plane scatterers in optimally doped cuprates, 
the basis of the SCTMA phenomenology was called into 
question in a work by (Nersesyan et ai, 1995), who 
pointed out that the neglect of multiple-impurity scatter- 
ing processes ("crossed diagrams"), which are of higher 
order in l/kpi oi the density of impurity sites rii in three 
dimensions, was not justified at asymptotically low en- 
ergies since these diagrams have an additional log E sin- 
gularity in 2D coming from the line nodes of the d^^-y^ 
gap. This leads to a breakdown of single-site perturba- 
tion theory and the need to resum at least a second class 
of crossed diagrams. This paper led to several attempts 
to solve the problem nonperturbatively (Nersesyan et al, 
1995; Pepin and Lee, 1998; Senthil and Fisher, 1999; 
Ziegler et al, 1996a), each claiming to calculate exactly 
the density of states of a disordered 2D rf-wavc siipercon- 
ductor, but arriving at dramatically different results. A 
combination of numerical work (Atkinson et al., 2000b) 
and analytical weak localization calculations (Yashenkin 
et al., 2001), reviewed in (Hirschfeld and Atkinson, 2002), 



led to the conclusion that the low-energy DOS of the 2D 
d-wavc state was quite sensitive to the symmetries of the 
disorder model and the normal state one-electron band 
employed in the calculations, in contrast to expectations 
based on standard localization calculations for normal 
metals. In addition, it was argued that the corrections 
to the SCTMA were unlikely to be relevant to experi- 
ments on real cuprates. 

While the SCTMA appears to be a reasonably suc- 
cessful approximation for properties of optimally to over- 
doped cuprate superconductors with nonmagnetic disor- 
der, it clearly cannot account for the magnetic states 
induced by potential scatterers as described in Section 
IV, nor other phenomena which must arise when pairing 
takes place in the underdoped phase near the Mott transi- 
tion. Here we investigate how the physics of impurities in 
correlated metals discussed in Sections IV and V is mod- 
ified by superconductivity. To establish the groundwork, 
we first give a brief account of the formation of low-energy 
impurity resonances (VLB) in a rf-wavc superconductor, 
and associated impurity bands by nonmagnetic poten- 
tials. We next discuss in Section VI. C the interference of 
these states in the presence of finite disorder, and the for- 
mation of the "impurity band" in the rf-wavc state. We 
then introduce electronic correlations in Section VI.D, fo- 
cussing primarily on results from the weak-coupling the- 
ory of impurity- induced moment formation as discussed 
in the previous section, since it can easily be extended 
to the superconducting state. We then compare a gen- 
eral picture which emerges from these treatments with 
experiments, both those which probe disorder-averaged 
properties, i.e. primarily states far from the impurities 
themselves (Section VI. E), and those which probe states 
near the impurity (Section VI.F). In Section VI.G we 
summarize ideas regarding the effects of disorder in the 
correlated d-wave state gleaned from both theory and 
experiment. 

B. Single impurity in a d-wave superconductor 

The BCS Hamiltonian for a pure singlet superconduc- 
tor can be written as: 

Ho^Yl ^k(ekT3 + Aferi)$k, (43) 

k 

where $k = (ckj^,cl^|), is a Nambu spinor. Analj^tic re- 
sults are presented for a parabolic band et = k'^ /2m, with 
corresponding d-wave order parameter Ak = Aocos2^ 
{(f) is the angle in momentum space which k makes with 
the 100 axis). The matrices are the Pauli matrices 
spanning particle-hole space. To the homogeneous su- 
perconductor we add a single impurity perturbation 

^fint = ^<l4k'$k' (44) 

kk' 

where Vk.k' = / dr e^'-^^^ ^''^ V{r) is the scattering po- 
tential Fourier component. For a density-like perturba- 



45 



tion which couples to the chemical potential locally, the 
Nambu form is Vk,k' = Vk,k'''3- The full Green's function 
for a single impurity is 



G(k,k') = Go(k)+Go(k)Fk,k'Go(k') (45) 
; Go(k)Fk,k"Go(k")Vk»,k'Go(k') + . . . . 



k" 



Go(k) + Go(k)fk,k'Go(k') 



(46) 



where the T-matrix is given by summing all single-site 
scattering processes, 



Tk.k' = Vk,k' + M.,k"Go(k")T4",k' + . . . (47) 

k" 

= Vk,k' +^Vk,k"Go(k")fk",k', (48) 



and all symbols A's are matrices in particle-hole space. If 
the impurity scattering is purely pointlike, V{r) cx S{r), 
the scattering is isotropic, Vk,k' = V, and T is seen to 
be independent of momentum. For the special case of a 
pointlike nonmagnetic impurity potential with strength 
Vq, it is convenient to define the cotangent of the s-wave 
phase shift ^o, with c = cot^o = {''^VqPo)'^ , with po the 
density of states at the Fermi level. The Green's function 
including the impurity is therefore 

G(k, k'; uj) = Go(k, c^)4k' + Go(k, c^)f (a;)Go(k', cu) (49) 

with Go(k, tj) — (uj — Hq)^^, and the T-matrix takes the 
form 



T = TqTq + T3T3 
Tq = go/{S+S^) 
T3 = (c-.g3)/(5+5_). 



(50) 



where go and 53 are the tq and Nambu com- 
ponents of the integrated bare Green's function 
(I/tt/Oo) X^k ^o(k, ij-'). This expression has resonances 
when 



S± = c-{g3T go) = 0. 



(51) 



Note that the off-diagonal components of the T-matrix 
do not occur because gi = 52 = due to the d-wave 
symmetry. Furthermore, in the special case of a particle- 
hole symmetric system, 33=0 and the resonance energy 
is determined entirely by go, which is given in the case 
of a circular Fermi surface by go{'-^) = —i / l^'^I'^^ ~ 
A^]^^/^ which for low energies lo ^ Ao takes the form 
Go[ijj) — — (7rcj/Ao)(log4Ao/a; -1- i). One may then solve 
Re S±{lo -\- iO"*") = and estimate the resonance width F 
on the real axis. In the case of strong scattering c <C 1, 
the resonance energy fijj^ and width F are 

±7rcAo 



F = 



2 log(8/7rc) 

TT^cAo 

41og2(8/7rc) 



(52a) 
(52b) 



0.1 

0.08 
„0.06 
-■0.04 

0.02 




-0.2 -0.1 



a) 




„ V 0) 


-0.2 -0 


1 0.1 0.2 




b) 




FIG. 36 Summary of LDOS results for 1-impurity problem 
on a tight-binding lattice, band given in reference (Norman 
et al., 1995), bond order parameter = 0.1, p = 0, Vb = 10, 
jSTiJl ~ 0.013: a) LDOS vs. u on the impurity site; b) LDOS 
vs. u on nearest neighbor site; c),d) LDOS map at resonance 
u) = Qq" and Lu = Q.'^ . Color scales in c) and d) are relative 
to the nearest neighbor peak heights shown in b). From (Zhu 
et al, 2003) 



This result was first obtained by (Balatsky et al., 1995) 
for the d-wave case, following earlier work on the simi- 
lar p-wave case by Stamp (Stamp, 1987). Note that the 
resonance becomes an undamped bound state only ex- 
actly at the Fermi level il = 0, when c — Q; for finite c 
there are two resonances whose energies — —^0 are 
symmetric in this approximation. As seen from (51), in 
general particle-hole asymmetric systems (G3 7^ 0), the 
resonance is tuned to sit at the Fermi energy for some 
value of the impurity potential Vq which is not infinite, so 
the term "unitarity" (as used here, Q.o = 0) and "strong 
potential" (Vq 00) are no longer synonymous (Atkin- 
son et al, 2000b; Joynt, 1997). 



1. Local density of states near impurity 

Although the energies of the particle and hole resonant 
states are symmetric, their spectral weights on a given 
site are usually quite different (Balatsky et al, 1995), 
such that locally only one resonance dominates. Define 
the local density of state (LDOS) as 

p[r,uj) = --lniGii{v,Y;uj + (53) 

TT 

with the total Green's function in the presence of the im- 
purity G(r,r';c^) = Go(r-r',c^) + Go(r,u;)f (w)Go(r',c^), 
where the 11 index refers to the normal electron part of 
the Nambu Green's function. Now it is easy to see that 
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any finite impurity potential acts as a local breaker of 
particle-hole symmetry which can lead to a pronounced 
peak in the LDOS. At the impurity site r = 0, for, e.g. 
a repulsive impurity potential, the resonant peak at neg- 
ative energy (holelike states) is large, and the positive 
energy (electronlike states) peak is relatively small, as 
seen in Fig. 36. The situation is reversed for the nearest 
neighbor site. Note, however, that for the large value of 
Vq chosen, the LDOS peak at resonance is much smaller 
on the impurity site than on the nearest neighbor site 
since the large potential excludes electrons from the im- 
purity site. 

Impurity bound states of this type have immediate ob- 
servable consequences for the local density of states, mea- 
surable by STM, and for the distributions of local mag- 
netizations, measured by NMR. The impurity modulates 
the local density of states of the homogeneous system 
p(r, uj) by 

Pimp(r,a;) = --Im [Go(r, a;)f(a;)Go(-r, a;)]ii (54) 

TT 

The impurity-induced LDOS falls off as along the 
nodal directions, and exponentially along the antinodal 
ones (Balatsky et al., 1995). The LDOS in the near field 
is more complicated, however: the nearest neighbor sites 
have peaks at ±f2o) with the larger spectral weight at 
-I- rip- In the crossover regime r ~ the LDOS is en- 
hanced along the node direction for holelike states, but 
is spread perpendicular to the node direction for elec- 
tronlike states. These spatially extended LDOS patterns 
are the fingerprint of the impurity-induced virtual bound 
states. In Fig. 36, we illustrate the LDOS pattern ex- 
pected for both particles and holes for a resonant state 
close to the Fermi level in this simple model. 

The existence of impurity bound states in the super- 
conductor will be shown below to be related to strong 
T-dependent upturns in the local susceptibility, so it is of 
interest to establish what happens in the current model 
as the temperature is increased and superconductivity 
destroyed. If the normal state is metallic, as assumed in 
the current model, weak Friedel-like oscillations in the 
density of states remain, but the resonant character of 
the impurity state is destroyed. If however the normal 
state also possesses a spectral gap, resonant states may 
remain. This is why impurity bound states found, e.g. 
in models of the spin gap phase of the cuprates(Gabay, 
1994; Khaliuhin et al, 1997; Nagaosa and Ng, 1995), with 
linear density of states analogous to the d-wave supercon- 
ductor, are quite similar to those in the superconducting 
state. Impurity measurements have in fact been proposed 
as a probe of the pseudogap state in the cuprates(Kruis 
et al., 2001). 

Because nonmagnetic impurities are pairbreaking in 
d-wave superconductors, the most important qualitative 
distinction with magnetic impurities does not exist. Nev- 
ertheless there are some interesting phenomena associ- 
ated with impurity spins in unconventional superconduc- 
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FIG. 37 Magnetization along (1,0) direction through impu- 
rity with d-wave superconductor described by tight-binding 
band t'/t = -.2, Vimp = WOt, g^BB/2 = 0.004t, U/t = 
0,1,1.5. 



tors, which have been amply reviewed in (Balatsky et al., 
2006). 



2. Effect on field-induced magnetization 

In the pure d-wave state with no residual quasiparticle 
interactions, the paired gas has a Friedel-type response to 
the impurity perturbation which leads to a spatially os- 
cillating LDOS p(r, w = 0) at the Fermi level with wave- 
length 2kF and envelope decaying over a length scale 
of ^0 = i^F/(7rAo). Consider a situation where the im- 
purity is nonmagnetic; there are then no corresponding 
spin density waves induced in zero field by correlation 
effects. However, since the local susceptibility in the su- 
perconducting state x(r, T) is proportional to the density 
of states at the Fermi level, such oscillations translate di- 
rectly into magnetization modulations in the presence of 
a nonzero applied field. Since the LDOS is always posi- 
tive, the magnetization is always in the same direction as 
the applied field. This can be seen in the U = curve of 
Fig. 37. In the presence of electronic correlations, these 
modulations are greatly enhanced, and a true alternat- 
ing magnetization appears, as discussed further below in 
Sec. VI.D. 



C. Effects of disorder on d-wave state 

1. Quasiparticle interference of impurity bound states 

As seen in Fig. 36, the strong nonmagnetic impu- 
rity in a d-wave superconductor in this approximation 

is characterized by a bound state with long range tails 
extending along the nodal direction. If many impurities 
are present, these states overlap and interfere, leading 
to a splitting of bound state energies and an accumula- 
tion of low-energy impurity-induced energy eigenvalues 
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which are spread out over a so-called "impurity band" 
(Micheluchi et ai, 2002). The formation of this band, 
which appears clearly in the disorder-averaged theory, is 
actually a rather subtle phenomenon compared to, say, 
the analogous impurity band in disordered semiconduc- 
tors, where a good deal can be understood by postulat- 
ing transport via overlap of spherically symmetric hydro- 
genic orbitals. In the c?-wave case the formation of the 
impurity band and the corresponding quasiparticle local- 
ization problem are strongly influenced by the fact that 
significant overlaps between two impurity states can take 
place only if the impurities are "oriented" with respect to 
one another such that the nodal wavefunction tails over- 
lap (Balatsky and Salkola, 1996). In fact, analyses of 
the two-impurity problem show that interference effects 
can take place over many tens of lattice spacings between 
optimally oriented impurities (Andersen and Hedegaard, 
2003; Morr and Stavropoulos, 2002; Zhu et ai, 2003). In 
Fig. 38, we show the local density of states due to many 
strong impurities using a realistic band structure, aver- 
aged over a few eigenstates within a narrow energy range 
abound the 1-impurity resonance energy (Atkinson et al, 
2003). The interference tails between some impurities are 
visible, while others have been split outside the "observa- 
tion bandwidth" due to the details of the local disorder 
environment. It is this complex system that the SCTMA 
seeks to replace by a disorder-averaged, translationally- 
invariant medium where the quasiparticle states of the 
pure system are broadened. 



2. Self-consistent T-matrix approach 

Since the SCTMA approach is already well-described 
in several original papers and reviews (Balatsky et at, 
2006; Buchholtz and Zwicknagl, 1981; Hirschfeld et ai, 
1986; Hussey, 2002; Joynt, 1997; Schmitt-Rink et ai, 
1986), we will only briefly summarize it here. We be- 
gin with the exact equations for the single impurity 
T-matrix (48). The self-consistent T-matrix approach 
(SCTMA) now extends these equations to an approxi- 
mate treatment where the disordered system is replaced 
by a translationally invariant dissipative medium de- 
scribed by a self energy, E(k, w) = riimpTi^.-k, where T 
is now a disorder-averaged quantity At the same time 
the full Green's function (^(k, k',w) is replaced by a dis- 
order averaged G{h,uj), the propagator for an electron 
moving in the translationally-invariant effective medium, 
given by 



(55) 



where now E is the sum of all disorder-averaged T- 
matrix diagrams, evaluated with the fully self-consistent 
disorder-averaged G on each Green's function line. The 
calculation is similar in spirit to the approximation of 
(Abrikosov and Gorkov, 1960), but valid for arbitrary 
impurity strength. 
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FIG. 38 Left: Density of states p(0) vs. position in a 70x 70 
lattice, showing interference of many impurities in band with 
no particle-hole symmetry. White dots represent strong po- 
tential scatterers. After W.A. Atkinson, unpublished. Right: 
schematic density of states of a d-wave superconductor, and 
impurity band of width 7. 



Rather than give details of the calculation of E, which 
have been amply reviewed elsewhere, we focus on the sim- 
plest example which displays the basic physics of a disor- 
dered d-wave superconductor, a finite density of identical 
strong pointlike scatterers. In this case the self-energy 
reduces to the analog of (50), 



E(k,t^) = F 



(c - g3)T3 + goTQ 



9l 



(56) 



where the impurity scattering rate parameter F = 
nimpi (tt/Oq) and 30^3 are now functionals of G rather than 

Gq. Let us furthermore assume particle-hole symmetry 
and unitary limit scattering c ^ 0. Then the self-energy 
reduces to the particularly simple result E — F/go- 
At the Fermi level, u) — > 0, if we assume the retarded 
self-energy Eg reduces to a small imaginary part —17, so 
that the self-consistency equation becomes, for a circular 
Fermi surface and d-wave order parameter Aq cos 2^, 



-17 



50(0) 



'tt jo 



(57) 



The elliptic integral can be performed and expanded for 
7 <C Ao to give the following self-consistency equation 
for 7 



7 

Ao 



. F /4A0 
2 Ao^ V 7 



(58) 



a transcendental equation whose solution for small F 
can be approximated to log accuracy as 7 ~ 0.63-\/FAo 
(Hirschfeld and Goldenfeld, 1993). The DOS at zero fre- 
quency is then given by 



Po 



-lmgo(O) = F/7 



(59) 



where the first equality is general but the second obtains 
only in the unitarity limit. 
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Thus the SCTMA generically predicts a finite den- 
sity of states at the Fermi level in the presence of dis- 
order for a (i-wavc superconductor. The dependence of 
this residual DOS on the concentration of defects ranges 
from the po ~ y/nimp behavior in the unitarity limit 
just exhibited, to a much weaker exponential dependence, 
p(0) exp(7rpo'^o(l + <^)/nimp) in the weak scattering 
(Born, c ^ 1) limit. The general form for the DOS 
expected over the whole frequency range in the strong 
scattering limit is shown in Fig. 38b. If the scattering is 
weaker, the weight in the impurity band moves to higher 
energies as the position of the single-impurity resonance 
itself moves away from the Fermi level, see Eq. (52a), 
and the Fermi level DOS p(0) decreases accordingly. 

The simplest measurement showing the strong influ- 
ence of impurities is the suppression of the critical tem- 
perature, which has been discussed to some extent al- 
ready in paragraph IV.D.2.b in connection with transport 
experiments. Within the SCTMA, the remarkable pre- 
diction is that for a rf-wavc superconductor with pointlikc 
scatterers the Tc suppression follows a form identical to 
the prediction of Abrikosov and Gor'kov (AG)(Abrikosov 
and Gorkov, 1960) 

where T^o is the transition temperature in the absence 
of disorder, and iIj{x) is the digamma fimction, and 
1/r is the scattering rate at the critical temperature 
due to impurity scattering alone, assuming that the im- 
purity potential is nonmagnetic and of zero range, for 
any scattering phase shift. For small concentrations, 
the initial slope of the Tc suppression is then Tc/T^a — 
1 — 7r/(4Trco). With increasing concentrations, super- 
conductivity is eventually destroyed (T^ = 0) when the 
scattering rate reaches the critical value 

= 7rTeo/27, (61) 

where 7 « 1.78. If Matthiessen's rule is obeyed (Section 
IV. D), the scattering rate is proportional to the change 
in resistivity due to disorder Ap and one expects an AG- 
type form for Tc vs. Ap. 

3. Extended impurities 

As discussed in Section IV. B. 2, impurities which sub- 
stitute for Cu in the CUO2 planes turn out to be strong 
scatterers of electronic states in the plane. While the 
effective bare potential induced by planar impurities is 
expected to be atomic scale (Wang et ai, 2005), it has 
been suggested that the correlations could give rise to 
a larger effective range (Section V.B). This is a possi- 
ble reason why the initial slope of the Tc suppression vs. 
residual resistivity curve in the case of Zn and Ni is found 
to be somewhat smaller (roughly a factor of 2) than pre- 
dicted by the SCTMA (see Section VI.E) (Franz et al, 



1996; Graser et ai, 2007; Haran and Nagi, 1996, 1998; 
Kuhc and Oudovenko, 1997; Kulic and Dolgov, 1999). 

On the other hand, there are many other sorts of de- 
fects in the cuprates, and in fact all cuprate materials 
are doped by out-of-plane ions which in most cases oc- 
cupy random positions in the crystal lattice or intersti- 
tial positions. It is to be expected that these "intrinsic" 
impurities are poorly screened, and therefore that they 
should have a large effective range of the scattering po- 
tential experienced by holes moving in the planes. In 
addition, they tend to be present in rather large con- 
centrations, indeed of order the doping itself. With the 
notable exceptions of YBa2Cu307+^ and YBa2Cu408, 
cuprates at optimal doping must therefore be consid- 
ered to have 0(10%) concentrations of "intrinsic" im- 
purities with extended potentials. Such defects will tend 
to scatter electrons preferentially in the forward direc- 
tion. The consequences for normal state one-electron 
properties were highlighted by (Abrahams and Varma, 
2000), for the critical temperature by (Kec, 2001), and 
for quasiparticle properties in the superconducting state 
at low temperatures by (Dahm et al., 2005a,b; Durst 
and Lee, 2000; Kampf and Devereaux, 1997; Nunner and 
Hirschfeld, 2005; Xiang and Wheatley, 1995; Zhu et ai, 
2004). Impurities of this type are assumed to be sufB- 
ciently weak such that they do not give rise to resonant 
effects at low energies, and are therefore not so easy to 
image with STM techniques. In addition, weaker poten- 
tials are not expected to induce significant local mag- 
netism (see Section VI.F), but there do not appear to be 
theoretical investigations along these lines available. 

(Zhu et al., 2004) and (Dahm et al., 2005a,b) concluded 
that forward scatterers have a tendency to sharpen spec- 
tral features near the antinodc of the d-wave gap, while 
broadening those near the node. This counterintuitive 
result can be understood by noting that if the scattering 
is sufficiently forward, scattering processes which begin 
with quasiparticles near the antinode cannot involve a 
sign change of the order parameter and are therefore not 
pairbreaking; this is the analog of Anderson's theorem 
for forward scattering in the d-wave state. On the other 
hand, if a quasiparticle is sufficiently close to the node, 
phase space for scattering involving sign change of Ak 
is available and spectral features will be broadened sim- 
ilarly to the normal state. 

Forward scattering impurities might be expected to 
have negligible effect on transport, since they cannot re- 
lax the momentum current. On the other hand, (Nun- 
ner and Hirschfeld, 2005) showed that scatterers of inter- 
mediate range could produce qualitatively different tem- 
perature dependences of the superconducting state mi- 
crowave conductivity, providing a possible explanation 
to the quite different results obtained for crystals of op- 
timally doped YBa2Cu307+5and Bi2Sr2CaCu208 • 
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4. Order parameter suppression and pairing potential impurities 

It is worth noting that the theoretical models above 
for the most part do not account for the local variation 
of the superconducting order parameter near the impu- 
rity site. The order parameter has a response to a lo- 
cal Coulomb potential, and the correct way to include 
it within a fully self-consistent Hartree-Fock theory was 
discussed by (Fetter, 1965), (Rusinov, 1968) and others 
for ordinary superconductors. In the context of uncon- 
ventional superconductors order parameter suppression 
around a single impurity was discussed by several au- 
thors (Atkinson et ai, 2000a; Choi and Muzikar, 1990; 
Franz et ai, 1996; Friesen and Muzikar, 1997; Shnirman 
et ai, 1999). These effects were neglected early on prob- 
ably because, in terms of the resonance effects on one- 
electron properties discussed in Section VLB, the self- 
consistent treatment does not lead to major differences. 
In the disorder-averaged theory it can have some impor- 
tant effects on transport as a new source of "off-diagonal 
scattering" (Hettler and Hirschfeld, 1999), however. 

There is a further, more exotic possibility, which is 
just beginning to be taken seriously at this writing. In 
the self-consistent approaches above, the BCS pair in- 
teraction constant A which appears in the gap equation 
is assumed to be homogeneous, and one considers the 
response of the order parameter to a diagonal poten- 
tial with T3 symmetry in Nambu space, i.e. a response 
(5Ak(r)Ti is produced. In principle, it has been recog- 
nized that an impurity can alter the local electronic struc- 
ture such that the pairing interaction itself is modified lo- 
cally (Larkin, 1970; Suhl et ai, 1962), leading to a A(r). 
Since coherence lengths of conventional superconduc- 
tors are much larger than microscopic distances a, these 
effects are typically averaged over in observable quanti- 
ties. In cuprates, however, the effects can be important, 
as shown by (Nunner et ai, 2005), who considered the 
effects of such impurities on the LDOS. Strong pairing 
impurities can produce resonant effects of their own with 
quite different characteristic LDOS patterns than con- 
ventional impurities, as shown by (Andersen et a/., 2006; 
Chattopadhyay et ai, 2002). 



D. Effect of correlations 

1. Single impurity 

The preceding discussion relates to the simplest the- 
ory of localized potential scatterers in a c?-wave super- 
conductor, and to some simple generalizations within 
the framework of BCS theory. While the approach has 
had considerable success at explaining properties of dis- 
ordered cuprates (see Sections VI.E,VI.F), it is expected 
that strong residual quasiparticle interactions neglected 
in this approach are still present even in the supercon- 
ducting state, and should influence particularly the un- 
derdoped state in proximity to the antiferromagnetic in- 
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FIG. 39 Spatial variations of the staggered magnetization 
Ms [(a), (c)], and LDOS spectra [(b), (d)] calculated for a 
32x 32 system. The impurity site is at (16,16). The solid 
lines in panels (b) and (d) are for (16,16), dashed lines for 
(17,16), dashed-dotted lines for (17,17), and dotted lines for 
(16,1). The upper [(a),(b)] and the lower panels [(c), (d)] are 
for {U,Vimp)^ (2,3), and (2.35,100), respectively. After (Chen 
and Ting, 2004). 



sulator. Many of the approaches to the problem of an 
impurity in a correlated system discussed in Section V 
have therefore been generalized to include pair correla- 
tions, and we present some of these efforts here. 

From the weak-coupling point of view, the question 
posed by (Ohashi, 2002) is, does the very substantial re- 
duction of the susceptibility in a d-wave superconduc- 
tor at low temperatures (x ~ pqT/Aq) prevent the en- 
hanced magnetization near the impurity site due to the 
coupling of spin fluctuations at q = and (tt, tt)? In 
fact, the enhancement of the susceptibility by a strong 
impurity relative to the system far from the impurity 
site is found to be considerably more important in the 
superconducting case, due to the resonant enhancement 
of the low-energy density of states already present even 
in the absence of correlations (Section VI. C. 2). The in- 
terplay of resonant density of states enhancements and 
correlation-induced magnetic moments is one of the inter- 
esting questions which arises only in the superconducting 
(and possibly pseudogap) states, to be addressed here. 

A phase diagram for the single-site impurity treated 
in the weak-coupling mean field model was presented by 
(Chen and Ting, 2004), who tried to relate the LDOS fea- 
tures probed by STM to the local magnetic state around 
the impurity (Fig. 39). They argued that the formation 
of a spontaneous S'^=l/2 moment was always associated 
with the substantial splitting of the small-?/ LDOS reso- 
nance, and concluded that since the resonances observed 
around Zn were unsplit (Pan et al., 2000b) -see Section 
VI. F. 2 below- they were necessarily associated with the 
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FIG. 40 Zero magnetic field phase diagram for a single im- 
purity of strength e — Vimp in the weak-coupling approach, 
with single band characterized by t' /t = —.2, n = 0.85, as 
a function of U. Values of S indicated are total values of S 
induced by an impurity in the paramagnetic state of a d-wave 
superconductor. After (Harter et al, 2007). 



Sz=0 state. The "local phase transition" from the 5^=0 
state to the S'^ = l/2 state, with concomitant spUtting of 
the LDOS resonance, is similar to the phenomenology of 
the singlet-doublet transition for a magnetic impurity in 
an s-wave superconductor (Sakai et ai, 1993; Sakurai, 
1970; Salkola et al, 1997). (Chen and Ting, 2004) also 
found in their numerical studies a Sz=0 phase which dis- 
played local SDW-type order, but (Harter et al., 2007) 
concluded that this phase was an artifact of the particu- 
lar finite-size system studied. The simpler phase diagram 
presented by the latter group is shown in Figure. 40, 
and indicates that the tendency to form a spontaneous 
extended moment increases either with correlations (U) 
or the bare impurity potential (Vimp)- It should be noted 
that mean field treatments of this type can at best de- 
scribe the tendency to form static moments. If the mean 
field model indicates an S=Q ground state, for example, 
this docs not preclude the existence of a fluctuating mo- 
ment at low energies. 

(Tsuchiura et al, 2001) performed similar calculations 
in a t— J model in zero field decoupled to allow rf-wave su- 
perconductivity, with a single site cut off from the rest of 
the system. The model was treated within the Gutzwiller 
approximation, amounting to assigning local bandwidth 
renormalizations to each bond. They found a sponta- 
neous moment phase with alternating magnetization near 
the impurity site, as well as a zero moment phase similar 
to the weak coupling calculations. In addition, they cal- 
culated the LDOS in each case. In the magnetized state, 
where a net magnetization corresponding to about 58% 
of the full moment 5=1/2 was found, they predicted 
a small 0(4 mcV) splitting of the resonance. Since this 
is not observed by STM (Section VI. F. 2), they also con- 
cluded that the impurity state at low T in the supercon- 
ductor does not carry a static moment. 

Similar conclusions were also reached by (Zhang et al, 



2002), who studied an Anderson Hamiltonian in the su- 
perconducting state with modified hybridizations con- 
necting the impurity with the four nearest neighbor sites 
of a Zn. A related model had been put forward ear- 
lier by (Polkovnikov et al, 2001), who began however 
from the Kondo limit, and argued that d-wave quasipar- 
ticles effectively screen the moment. These models and 
the controversial interpretation of the STM results have 
been reviewed in (Balatsky et al, 2006). 

Tsuchiura's work was criticized by (Wang and Lee, 
2002), who argued that the less than full moment found 
was indicative of a failure of the mean field procedure 
adopted. These authors also studied a vacancy in the 
superconducting t — J model, but used a slave bo- 
son mean field decomposition and spatially unrestricted 
Bogoliubov-de Gennes calculation of all local slave boson 
amplitudes. In addition, an additional mean field spin 
coupling aJ was added, so it is difficult to judge the ori- 
gin of differences between the two similar calculations. In 
any case, the (Wang and Lee, 2002) result showed either 
a singlet S = or doublet 5=1/2 ground state depend- 
ing on the additional spin coupling; no singlet phase with 
alternating magnetization was found. The low energy 
LDOS resonance was found in both cases, again slightly 
split in the doublet state. Wang and Lee also studied 
the behavior as a function of doping, and found that the 
static moment (doublet) state is much easier to stabilize 
in the underdoped phases. More recent numerical work 
on impurities in the superconducting t — J model by vari- 
ational MC arrived at a similar qualitative picture of the 
local spin correlations (Liang and Lee, 2002). 

Before closing this summary of the single impurity in 
zero applied field, we remind the reader that the 3^ = 
phase of Fig. 40, similar to the singlet phase of (Wang 
and Lee, 2002), while apparently less interesting than 
the "spontaneous moment" phase is in fact representative 
of most of the cuprate phase diagram including the su- 
perconducting phase (the exception is the phase labelled 
"spin glass" in Fig. 2). That is, as emphasized in Section 
IV already, there is no experimental evidence at optimal 
doping for static magnetic moments around nonmagnetic 
impurities in zero field. In highly underdoped samples at 
low temperatures, a glassy response is found by fiSR (see 
Section VI. E). In most cases of interest, however, we ex- 
pect that the situation is like that depicted in Fig. 39 (a), 
with zero magnetization until the field is applied. Then, 
as seen in Fig. 37, the presence of electronic correlations 
strongly enhances the natural Friedel-like response of the 
(paired) Fermi liquid such that the sign of the magnetism 
actually alternates from site to site. The integrated mag- 
netization in a field of a few Tesla may be less than the 
full spin- 1/2 moment, however. The reader may well ask 
what the theoretical justification is for using a Fermi liq- 
uid formalism to calculate the magnetic response in the 
normal state, where experimental evidence points to the 
absence of Landau quasiparticles. On the other hand, 
it is presumably more justified in the superconducting 
state due to the collapse of the scattering rate at the 
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FIG. 41 Magnetization histogram for a 34 x 34 system with 
U=1.75, and 1,4,8 and 12 impurities at T = 0.013 and 
giJ.BB/2 = 0.004. After (Harter et al, 2007). 



phasizcs the competing order aspect of the phenomenon, 
ascribing inhomogeneity at low doping to competing su- 
perconductivity and antiferromagnetism, modulated by 
disorder (Andersen et al, 2007; Atkinson, 2005, 2006; 
Harter et al, 2007; Kontani and Ohno, 2006; Mayr et al., 
2006) Other authors work at strong coupling and empha- 
size the effects of projection leading to a narrowing of the 
local band near the impurity (Garg et al., 2006; Wang 
and Lec, 2002). Although these works differ depending 
on details of parameter choice, whether or not long-range 
interactions are treated, etc., they agree that in the pres- 
ence of sufficiently strong interactions and sufficient dis- 
order strength, "patches" of (generally incommensurate) 
magnetic order can be nucleated by impurities. 



transition (Basov and Timusk, 2005) and the recovery of 
sharp quasiparticle spectral features seen, e.g. in ARPES 
(Damascelli et al., 2003). 



2. Many-impurity interference effects. 

(Chen and Ting, 2004) briefly investigated 2-impurity 
interference effects in the system with correlations 
treated in weak coupling, finding that the magnetization 
on sites near an impurity in the 5*2 = 1/2 phase could 
be enhanced or suppressed in the presence of a second, 
identical impurity, depending on its orientation with re- 
spect to the first. (Harter et al., 2007) explored this ques- 
tion further in the context of NMR. They looked first at 
the probability distribution of magnetizations for a single 
(field-induced) impurity, showing that sites ranging from 
those close to the impurity up to 12-15 lattice constants 
away contributed to pronounced satellite lines at certain 
magnetization values, but that these lines were immedi- 
ately smeared upon addition of as little as a few tenths of 
a per cent impurities (Fig. 41). The reason has to do with 
the fact that those satellites which, when broadened, ul- 
timately determine the width of the NMR line arise from 
sites located in the tails of the nodal quasiparticle bound 
state wavefunctions. As discussed in Sec. VI. C, these 
tail states interfere strongly with other such states when 
impurities are oriented close to 45° with respect to each 
other. This applies equally well to the spin as well as the 
spatial parts of these wavefunctions. This proposal will 
be discussed further in comparison to experimental data 
below. 

A number of authors have discussed the effects of finite 
disorder on d-wave superconductivity in the presence of 
electronic correlations. These issues are crucial to the de- 
bate about the origin of the pseudogap, the properties of 
the "spin glass phase" which coexists with superconduc- 
tivity in several of the cuprates, and the question of the 
origin of nanoscale gap inhomogeneity in these materials. 
One school begins from the weak coupling end, and em- 



E. Experiments probing disorder-averaged bulk properties 

The SCTMA and its extensions describe the low- 
energy resonant enhancement of the density of states by 
strong impurities in dirty rf-wave superconductor. We 
now explore the extent to which experiments can con- 
firm this general picture, and particularly those examples 
where the SCTMA, which neglects correlations, is insufS- 
cient. There are many fundamental experiments on bulk 
cuprate samples which probe the low-energy DOS in the 
presence of strong disorder, including specific heat, ther- 
mal conductivity, magnetic penetration depth and mi- 
crowave conductivity measurements. Many of these have 
been reviewed by (Hussey, 2002), who made the case for 
a qualitative overall agreement of the simple dirty d-wave 
ideas with experiment, but pointed out some discrepan- 
cies as well. 

As discussed in Section IV, early Tc suppression exper- 
iments were inconclusive regarding the symmetry of the 
superconducting order parameter in the cuprates. More 
decisive were penetration depth and NMR experiments 
which probe low energy excitations and revealed the ex- 
istence of the impurity band characteristic of gap nodes 
and strong potential scattcrcrs. Fig. 42 shows Lon- 
don penetration depth measurements reported on good 
quality Zn-substituted single crystals of YBCO by (Bonn 
et al., 1994), together with a fit to the theory of the pen- 
etration depth in a d-wave superconductor with unitary 
scatterers (Hirschfeld and Goldenfeld, 1993; Hirschfeld 
et al., 1994). The increase of the zero-T penetration 
depth and the T ^ crossover observed reflect directly 
the residual density of states created by strong impuri- 
ties. 

In principle, the asymptotic electronic linear-T term 
in specific heat experiments probes the residual DOS 
as well. From early days in cuprate physics, "anoma- 
lous" linear-T terms C(T) ~ T for T < in the 
superconducting state were reported for polycrystalline 
and later single crystal samples, and varying field de- 
pendences were reported as well (Hussey, 2002). These 
were shown to be qualitatively consistent with different 
levels of disorder by (Kiibert and Hirschfeld, 1997). Sys- 
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FIG. 42 Fits to penetration depth measurements of Hardy et 
al (Bonn et al., 1994) on crystals of near-optimally doped 
YBCO substituted with Zn concentrations 0, 0.15%, and 
0.31%, using dirty d-wave theory (Hirschfeld et al., 1994) with 
scattering parameters F as indicated. 



tcmatic studies with Zn doping in single crystals were 
performed by several groups (Chang et al., 2000; Loram 
et al, 1990b; Momono and Ido, 1996; Sisson et al., 2000), 
verifying the expected dependence of the residual density 
of states on impurity concentration. 

(Sisson et al., 2000) were the first group to do a sys- 
tematic comparison of the magnetic field dependence of 
Zn-doped samples with the Schottky anomaly expected 
on the basis of the assumption of local moments of or- 
der 0.3/xb, the typical order of magnitude found for these 
moments above Tc in optimally doped samples in early 
direct magnetization experiments (Mendels et al., 1994b; 
Zagoulaev et al., 1995) on polycrystalline samples. No 
such effect was observed, leading to the conclusion that 
any moments present must be Kondo screened by d-wave 
quasiparticles (Cassanello and Fradkin, 1997, 1996; In- 
gersent, 1996). 

We discuss now the effect of impurities on the super- 
fluid condensate from /xSR experiments. In high trans- 
verse magnetic fields, /zSR experiments allow one to mea- 
sure the siiperfluid density rig through the muon spin 
relaxation rate a due to the magnetic field penetra- 
tion in type II superconductors. At low temperatures, 
the ratio ns{T 0)/m*(, is observed to decrease with 
increasing impurity concentration, as shown in LSCO, 
YBCO, and YCaBCO for Ni, Zn (Bucci et al, 1994) 
(Bernhard et al., 1996) (Nachumi et al., 1996) and Li 
impurities (Bobroff, 2005). In (Bernhard et al., 1996), 
the impurity concentration dependence was observed to 
follow the expected behavior arising from pair-breaking 
due to Zn defects considered here in the unitary limit 
in a d-wave superconductor (Bernhard et al, 1996). In 
contrast, Uemura's group interpreted this decrease in a 



Bose condensate picture where Tc is proportional to rig 
(Uemura, 2004). The reduction of Tc and is then ex- 
plained in a "Swiss cheese" model where each Zn creates 
a non-superconducting region around it with an area of 
TT^^ with ^ being the superconducting coherence length 
(Nachumi et al., 1996). Even though this interpretation 
seems drastic, it fits the data well with no adjustable pa- 
rameters. However, local probes such as STM and NMR 
reveal a much more complex situation with a progres- 
sive spatial evolution of the electronic properties such as 
density of states or local magnetism around each Zn, as 
presented in the following section. This contradicts the 
simple granular picture with normal state islands around 
each Zn and superconducting condensate beyond. The 
"Swiss cheese" model is also in contradiction with the 
successful disorder-averaged (SCTMA) theory which pre- 
dicts square root dependences on impurity concentration, 
not linear. 

Finally, we discuss neutron experiments. The earli- 
est measurements on Zn-substituted YBCO crystals were 
performed by Sidis et al. (Sidis et al., 2000, 1996) In the 
presence of Zn, they observed a pronounced increase of 
low-energy spectral weight x(Q,ct;), where Q = 7r,7r, as 
well as a suppression and shift of the 40 meV resonance 
mode (Fig. 43). The spin gap at ~ 30 meV in optimally 
doped YBa2Cu307+5was found to be still present in the 
Zn-doped system. (Sidis et al., 1996) speculated that the 
sites far from the Zn would contribute to the resonance 
mode and spin gap, while the regions close to Zn would be 
responsible for the low energy new excitations. Indeed, 
these excitations can be associatcid with the enhanced 
relaxation rate observed both in the normal and super- 
conducting states in regions close to the impurity, as is 
seen in Fig. 19. However the distinction between the 
two regions is not straightforward (because of the long 
range DOS effects in the superconducting state), unless 
one uses a local probe such as NMR. This picture is con- 
sistent with the existence of a finite DOS at the Fermi 
level and with NMR findings detailed below. 

(Lake et al., 2002) reported evidence of dopant induced 
ordered static magnetism in "pure" underdoped LSCO 
samples with nft=10%. No such effects were observed 
near optimal doping (Lake et al., 2001). The basic results 
are summarized in Figure 44. In zero field, weak static 
magnetism with an ordered moment of /U ~ 0.1/xb/Cu 
was observed at a wave vector close to but not equal to 
(tt, tt). The strength of the moment corresponds roughly 
to a moment of spin-1/2 per dopant, but if the moment 
were concentrated exclusively on the random dopant site, 
would not lead to an observable signal. It therefore sug- 
gests that some fraction of the random dopants nucleate 
clouds of staggered polarization which then interfere to 
give an ordered state at low T. Note there is no evidence 
for such static spontaneous magnetism in any nominally 
pure YBCO sample except possibly at extremely low dop- 
ing (/uSr, Sonier...). 

(Kimura et al., 2003) studied explicitly the effect of Zn 
substitution in a large optimally doped LSCO crystal. 
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FIG. 43 Imaginary part of dynamical susceptibility x"{^ ~ 
TT,-K,uj) on YBCZnOr extracted from neutron scattering at 
T=10.5 K. Circles and squares are independent measurements 
in different energy ranges, open and closed symbols corre- 
spond to different methods of determining x" on same sam- 
ple. The results obtained at 4.5 K in pure sample at the same 
O content given for comparison by crosses. 
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FIG. 44 The dependence on temperature of the ordered mag- 
netic moment squared deduced from neutron measurements 
on LSCO (Lake et at, 2002). The zero-field signal (blue cir- 
cles) increases gradually below rc(II = T), and the dashed 
line is a guide to the eye. Also plotted is the field-induced 
signal (equal to the signal measured in-field minus the zero- 
field signal) for fields of H = 5 T (yellow circles) , f T (green 
circles) and f4.5 T (red circles). 



At zero and small Zn concentrations, no static signal was 
seen. However, when the Zn concentration reached 1.7%, 
they too found an elastic peak in the neutron response 
x" near (tt, tt) which increased by about 15% with de- 
creasing temperature. This group argued that the static 
correlations arose from the long range antiferromagnetic 
coherence among the induced moments around different 
Zn ions. A magnetic signal was also reported in Zn- 
substituted optimally doped LSCO by (Nachumi et al, 
1996) at very low temperatures. Thus it appears that 
moment formation and staggered ordering is significantly 
more favored in LSCO than in YBCO, and the tendency 
thereto is enhanced as one underdopes the samples. It 
is interesting to ask what the main difference between 
the two materials is. One possibility is the different way 
in which the two are doped, leading to an "intrinsic" Sr 
disorder in LSCO at the 10% level (see Section IV). Re- 
cently, (Andersen et ai, 2007) showed that in an increas- 
ingly random environment the probability of the local 
electronic structure being favorable at certain locations 
for induced impurity moments increases, and that when 
the concentration of these "active" defects is sufficient, 
an ordered state appears. The "order by disorder" phe- 
nomenon observed by neutrons, the distinction between 
YBCO and LSCO and other intrinsically doped cuprates, 
and the enhancement of the spin glass behavior by the 
addition of Zn, emerges from this analysis. More theoret- 
ical work on the magnetic properties of heavily disordered 
materials with out-of-plane defects is needed, however. 

F. Experiments probing local properties 

Local measurements such as NMR allow one to distin- 
guish between regions close and far from the impurity. 
However, in the case of Cu and O Nuclear Magnetic Res- 
onance, Gd Electron Spin Resonance, or Yb Mossbauer 



experiments, it is not straightforward to separate the two 
contributions, but requires careful analysis. By contrast, 
STM or Li NMR measurements give a direct insight into 
the charge and spin behavior on the nearest-neighbor 
sites. It is only the comparison between the former and 
latter types of experiments which will permit the reso- 
lution of spatially-dependent effects of an impurity in a 
superconducting host. 

1. Far from the defect 

In the normal state, the existence of local paramagnetic 
moments around Zn, Li and other impurities substituted 
at Cu site was shown by NMR (Mahajan et al., 1994) and 
has been reviewed in Section IV. In the superconducting 
state, similar experiments are not as straightforward be- 
cause of the presence of the vortex lattice induced by the 
applied magnetic field of a few Tesla. This vortex lattice 
results in a distribution of the local fields at the nuclei 
sites, which mixes with the impurity effect itself. In the 
case of Cu NMR, such an effect is negligible since the 
Cu hyperfine coupling is very large, so that any impurity 
effect on the DOS or the magnetism dominates that of 
the vortex lattice. The Cu NMR shift measured in Zn 
and Ni substituted YBCO is reported in Fig. 45 (Ishida 
et al., 1993). 

As already demonstrated in the normal state, the NMR 
shift of Cu (or O or Y) probes regions far from the defects, 
while regions close to the defects contribute to the width 
and/or induce additional satellites to the main NMR line. 
Above Tc, in the presence of Zn or Ni, no additional shift 
of the main line could be observed, indicating that re- 
gions far from the defects are not affected by substitu- 
tion. In contrast, below Tc in the case of Zn only. Fig. 
45 shows an additional shift induced by the impurities. 
This shift has been interpreted as due to an induced DOS 
close to Fermi Level by Zn. Indeed, NMR static measure- 
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FIG. 45 T-dependence of ''^Cu Knight shift in YBCO7 with 
(left) Zn or (right) Ni substituted for Cu (Ishida et al, 1993). 



ments allow one to probe the DOS at the Fermi Level, as 
demonstrated in the case of metals where the NMR shift 
measures the Pauli susceptibility (Slichter, 1998). In the 
present case, this additional shift merely represents the 
residual DOS far from impurities, which confirms that 
only Zn induces a sizable DOS at the Fermi level. Similar 
results were obtained using Gd electron-spin resonance 
(Janossy et al., 1994). Dynamical measurements such as 
longitudinal nuclear relaxation Ti which can be measured 
on ^"^Cu nuclei or on the Yb nuclei by Mossbauer effect 
also reveal this DOS induced effect by Zn (Hodges et at, 
1995; Hotta, 1993; Ishida et al., 1993). (Ishida et al, 
1993) analyzed Ti as originating from two contributions, 
from regions either close or far from Zn. They interpreted 
the Ti variations with impurity concentration in a d-wave 
model where the impurities induce some residual density 
of states. However, as already argued in the normal state 
(Section IV.C.4) (Itoh et al., 2003), the identification of 
the two contributions is not clear. Furthermore, an ad- 
ditional contribution to the relaxation may be expected 
from antiferromagnetic spin fluctuations enhanced in the 
vicinity of the impurity even in zero applied field (Itoh 
et al., 2003, 2001) which has the same origin as that stud- 
ied on Li or Y in Section IV. 




FIG. 46 Experimental STM impurity spatial patterns at res- 
onance and associated spectra, (a) observed differential con- 
ductance around Zn impurity at -1.2 mV; (b) comparison of 
spectra directly above Zn impurity site (red) and at site far 
from impurity (blue); (c) observed STM intensity around Ni 
impurity site at 9 mV (inset: pattern at -9 mV); (d) com- 
parison of spectra directly above Ni impurity site (red) and 
away from Ni (blue); (e) observed STM intensity around na- 
tive defect, possibly Cu vacancy at -0.5 mV; (f) comparison 
of spectra directly above native defect (thick solid line) and 
away from defect (thin solid line). Note that in the images 
(a), (c), and (f) the Cu-O bond is at 45° with respect to the 
horizontal. 



2. Near the defect 

The most local picture we have of impurity states in 
cuprates is provided in principle by scanning tunnelling 
microscopy (STM). Many aspects of the STM results 
have been reviewed recently by (Balatsky et al., 2006). 
The first great success of the STM technique was the ob- 
servation of resonant states around magnetic impurities 
in conventional superconductors (Yazdani et al., 1997), 
followed by nonmagnetic defects at low temperatures in 
the cuprate Bi2Sr2CaCu208 (Hudson et al., 1999; Pan 
et al., 2000a; Yazdani et al., 1999). The latter discovery 
confirmed earlier theoretical proposals that such states 
should be reflected in the LDOS of d-wave superconduc- 
tors (Balatsky et al, 1995; Byers et al, 1993). The HTS 
compound Bi-2212 has been used for most STM studies 
because of its high-quality surfaces when cleaved between 



BiO layers. Scanning spectroscopy measurements were 
performed on Zn and Ni impurities (Hudson et al, 2001; 
Pan et al, 2000b), as well as on native defects (Hudson 
et al, 2003, 1999) which apparently occur naturally in as- 
grown crystals. These results are summarized in Fig. 46, 
where it is seen that the different defects have different 
resonant energies within the superconducting gap, and in 
the case of Ni a double-peak structure is observed. In ad- 
dition, the spatial patterns when the surface is scanned 
at the resonant energy appear to differ from impurity to 
impurity, and in the case of Ni show a striking rotation 
between positive and negative bias. 

''Li NMR experiments allow further insight into the 
defect's immediate neighborhood. Li substituted at the 
Cu site is identical to Zn in the sense that it induces sim- 
ilar Tc reduction and local induced magnetism above 
(Bobroff et al, 1999) (see Section IV. The Li NMR shift 
has been demonstrated to probe only the local spin den- 
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FIG. 47 T variation of the Li NMR shift '^K in YBCO (Bo- 
broff et ai, 2001). Top: underdoped 06.6 sample, with Curie- 
Weiss fit (solid line). Bottom: optimally doped 06.97 sam- 
ples, Kondo model susceptibility fit to data above (solid line) 
and below Tc (dashed line). 



sity on its four neighboring Cu, both in the underdoped 
and optimally doped regime. It displays a Curie- Weiss 
like behavior which signals the appearance of the induced 
staggered magnetism around the nonmagnetic impurity. 
In the superconducting state of underdoped YBCO, the 
Curie behavior observed above indicating the presence 
of such local moments on nearest neighbor Cu persists 
below Tc with no change, as shown in Fig. 47 (Bobroff 
et ai, 2001). This is clear evidence that induced mag- 
netism on the nearest neighbor site is not affected in the 
underdoped case by superconductivity. By contrast, a 
dramatic increase of the Li shift is observed below Tc in 
an optimally doped sample (Fig. 47). This was originally 
interpreted as a reduced Kondo screening of the moment 
in the superconducting state since it can be fitted with 
the same Curie- Weiss behavior as above Tc but with a re- 
duced Curie- Weiss temperature (Bobroff et ai, 2001). As 
will be shown below, the LDOS peak detected by STM 
experiments can contribute as well to the induced local 
susceptibility near the defect. 

The LDOS pattern predicted for a strong impurity in 
a d-wave superconductor with noninteracting quasiparti- 
cles was shown in Fig. 36. The crucial qualitative dif- 
ference between the effect predicted and observed is the 
existence of an intensity maximum on the central site of 
the impurity pattern in the case of Zn, which is impossi- 
ble in the naive theory for potentials Vimp larger than the 
bandwidth, since electrons are effectively excluded from 
this site. Theoretical efforts to understand this paradox 
have fallen into three categories. The first is specific to 
the STM technique and relies on the fact that the im- 
purity states are localized in the Cu02 plane, two layers 
below the BiO surface; the intervening layers are then ar- 
gued to provide a primary tunnelling path which samples 
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FIG. 48 ^^O NMR lines in YBCO7 at T=15 K for three dif- 
ferent in-plane concentrations of Zn 2;p=1.5, 3 and 6% shown 
relative to the pure case (in gray). From Ref. (Ouazi et ai, 
2006). 



not the Cu below the tip, but preferentially the 4 nearest 
neighbors (Martin et al., 2002; Zhu et ai, 2000). Some in- 
direct support for this point of view has been provided by 
density functional theory (Wang et ai, 2005). A second 
class of approaches (Polkovnikov et ai, 2001) attributes 
the redistribution of spatial spectral weight to nonlocal 
Kondo screening, an effect neglected in the calculations 
described in Section VI. D. It was argued in (Balatsky 
et ai, 2006), where the controversy is reviewed in much 
more detail, that models of this kind require unphysi- 
cally small nonmagnetic compared to magnetic interac- 
tions to explain STM spectra. Finally, one can obtain 
LDOS patterns similar to experiment for both Zn and Ni 
by assuming a distribution of site potentials to tune the 
weights of on-site and nearest-neighbor LDOS (Flatte, 
2000; Tang and Flatte, 2002, 2004). While the bare im- 
purity potential is much shorter range, it may be that 
the phenomenological parameters used in these models 
are approximations to dynamically generated ones in a 
more complete theory, see Section V. There is currently 
no consensus as to which of these explanations, if any, 
is appropriate to describe both NMR and STM impurity 
results. Resolution of this "paradox" is crucial to the in- 
terpretation of STM images and the exploitation of their 
huge potential to provide local information on cuprates 
and other strongly correlated materials. 

In order to resolve the impurity effect beyond the near- 
est neighbor Cu sites, one has to use ®^Y, ^'^Cu or ^^O 
NMR as demonstrated in the normal state (see Section 
IV) (Ouazi et ai, 2006). In the superconducting case, 
only a ^^O NMR study has been performed, leading to 
typical spectra displayed in Fig. 48. They show that a 
nonmagnetic impurity such as Zn or Li induces an asym- 
metric broadening of the line together with a shift to- 
ward higher frequencies (Ouazi et ai, 2006). The latter 
effect is similar to that measured with Cu NMR and is 
related to the residual DOS close to Ep induced by de- 
fects sitting far from the measured nuclei . On the other 
hand, the Zn-induced broadening is due to the induced 
magnetism and LDOS closer to the defect. For a more 
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quantitative analysis, the vortex lattice effects have to be 
separated from the impurity effects by carrying out the 
experiment at different applied fields. Indeed, the field 
distribution due to the vortex lattice leads to a roughly 
field-independent broadening of the line, while both im- 
purity induced LDOS and paramagnetic effects lead to 
a broadening linear in field. The fact that Zn or Li 
broaden the line on both high and low frequency sides 
is evidence for the alternating behavior of the induced 
cloud, similarly to the normal state. However, the broad- 
ening is not symmetric anymore as in the normal state, 
but asymmetric towards high frequencies. This can be 
understood as due to the additional LDOS contributions 
which leads to a positive shift. This is confirmed by the 
fact that the Ni induced broadening does not show such 
an asymmetry, as expected since it does not induce any 
LDOS close to Fermi Level in contrast with Zn (see (d) 
of Fig. 46). In summary, Cu, O and Li NMR experiments 
imply that a nonmagnetic impurity induces simultane- 
ously a staggered magnetization similar to that observed 
in the normal state with typical extension of a few cell 
units, together with a spatially decaying LDOS contribu- 
tion which is still sizable far from the impurity. 

The experiments presented above in both the super- 
conducting and normal states clearly require an expla- 
nation in terms of a combination of LDOS enhancement 
of the spin susceptibility and an induced magnetic mo- 
ment around the impurity sites. While it is clear that 
due to the resonant nature of the LDOS effects in the 
superconducting state, these effects will be more impor- 
tant below Tc, attempts to explain the superconducting 
state data without invoking induced moments must en- 
counter difficulties. For example, (Williams et al., 2000) 
proposed that quasiparticle resonant states, correspond- 
ing to those imaged by STM experiments around Zn 
atoms in BSCCO-2212 (Pan et ai, 2000b), might be 
entirely responsible for the enhanced magnetic response 
near Zn seen in NMR. (Chang et al, 2004) then argued 
that for a single nonmagnetic impurity, the temperature 
dependence of the observed spin-lattice relaxation time 
and Knight shift could indeed be qualitatively under- 
stood in terms of LDOS enhancement due to impurity 
bound states alone. First of all, significant T-dependent 
enhancements of local susceptibilities near nonmagnetic 
impurities occur in the normal state of optimally doped 
cuprates as well. An "LDOS-only" approach cannot ac- 
count for this since impurities do not produce LDOS res- 
onances in the normal (metallic) state. Secondly, the 
NMR experiments on these materials clearly show that 
the magnetization near a Zn alternates in sign in both the 
superconducting and normal states, as discussed above. 
At low T, where the homogeneous polarization of the d- 
wave superconductor is negligible, the Priedel oscillation 
in the LDOS induced by an impurity must perforce be 
entirely positive, i.e. cannot give rise to the observed 
broadening on the negative side of the line. Finally, the 
"Knight shift" calculated in Ref. (Chang et ai, 2004) 
is defined to be a local susceptibility enhancement very 



near the impurity; in fact the measured Knight shift in 
experiments where the nucleus is distinct from the impu- 
rity itself is the shift of the total NMR line, determined 
by sites far from the impurities. 

These two effects can be understood in a unified pic- 
ture in terms of the enhanced local susceptibility due to 
the effect of correlations on the resonant impurity state 
(Harter et ai, 2007). In this study, from the magnetiza- 
tion histogram computed and displayed in Fig. 41, oxy- 
gen NMR lines could be computed for various impurity 
concentrations including the effect of vortices. A quanti- 
tative agreement is found with the experimental results 
shown in Fig. 48, tending to confirm the validity of this 
picture. 



G. Summary 

We have reviewed the simple theory of nonmagnetic 
impurities in d-wave superconductors, discussed some 
preliminary attempts which have been made to include 
electronic correlations beyond BCS pairing in this frame- 
work, and compared the results with experiments. The 
most important difference between the superconducting 
state and the normal state of optimally to overdoped 
cuprates is the existence of low energy states of a quasi- 
particle virtually bound to nonmagnetic potentials as, 
e.g., Zn impurities. These bound states should exist even 
in the d-wave superconductor without additional corre- 
lations, and have weak paramagnetic character to the 
extent that they enhance the LDOS locally around the 
impurity sites. Such LDOS enhancements have indeed 
been observed by STM experiments near impurities (al- 
though the spatial distribution of LDOS weight is not 
completely understood). NMR experiments find, how- 
ever, that the induced spin polarizations near these sites 
are much larger than expected according to any calcula- 
tion relying on BCS-type models alone, and we are forced 
to the conclusion that the usual "dirty d-wave theory" 
must be supplemented by the interactions which lead 
to paramagnetic moment formation around the impurity. 
As discussed in Section IV for the normal state as well, 
electronic correlations of antiferromagnetic character are 
essential to account for the observations. 

For optimally doped cuprates, the moments are rela- 
tively weak and well screened, but as the superconduct- 
ing state is entered NMR measurements of the local sus- 
ceptibility show that the opening of the superconducting 
gap reduces this screening, as expected. This behavior is 
reminiscent of Kondo screening, but has also been repro- 
duced by theories which neglect processes which lead to 
the usual Kondo effect. It appears furthermore that a de- 
scription beyond the single-impurity picture, accounting 
for interference of individual impurity states, is necessary 
to understand the details of experiments with impurities 
at per cent level or higher. Experiments such as specific 
heat, NMR Knight shift, and penetration depth which 
probe the average density of states primarily far from 
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the impurity sites are sensitive to the details of this in- 
terference and to the formation of an "impurity band" in 
the superconducting state. These appear to be in rough 
agreement with each other and with the predictions of 
the SCTMA. 

On the other hand, at low doping it is to be ex- 
pected that discrepancies will increase in underdoped 
samples as the magnetic character of the disorder is en- 
hanced. In this case, strong moment formation has al- 
ready taken place by the time the superconducting tran- 
sition is reached, and the local susceptibility is found to 
be Curie-like down to the lowest temperatures, by NMR 
measurements. Fewer experimental data are available to 
support this premise, however. In disordered cuprates 
such as LSCO, a range of experiments inchiding NMR 
and /iSR have observed the freezing of spins in an exten- 
sion of the "spin glass" phase into the superconducting 
states. These appear to be due to disorder-induced mag- 
netic centers which at sufficient concentration or correla- 
tion strength (undcrdoping) interfere coherently to create 
long-range magnetic order visible to neutrons. The mag- 
netic disorder should scatter quasiparticles efficiently as 
in the normal state, and lead to a modification of "dirty 
d-wave" predictions of superconducting transport at low 
T, as already observed for thermal conductivity measure- 
ments (Andersen and Hirschfeld, 2006). 



VII. CONCLUSIONS 

The study of strongly correlated electron systems is 
quite challenging and requires the most sophisticated ex- 
perimental and theoretical techniques. The primary ap- 
proach emphasized throughout this manuscript has been 
the investigation of the response of such systems to the 
perturbation introduced by point defect impurities. As in 
the classic case of Friedel oscillations in normal metals, 
the study of defects is valuable because they have the po- 
tential to probe aspects of the underlying pure system, 
e.g. magnetic and charge correlation lengths, sometimes 
difficult to do in situations when the fundamental na- 
ture of the states being probed is not always understood. 
Thus one of our primary motivations to study defects in 
such systems has been to try to understand what impu- 
rities tell us about the pure materials. This has conse- 
quently led us to show how extrinsic properties can be 
"subtracted" from an analysis of experiments on "pure" 
samples, in order to reveal the generic features. 

Impurity states in strongly interacting materials are 
most unusual, and manifest a rich variety of phenom- 
ena depending on dimensionality, strength of correla- 
tions, carrier doping, temperature and magnetic field. 
The quantum impurity problem of a spinless potential 
in a correlated host has already attracted a great deal of 
theoretical attention, and promises to pose as much of 
a challenge as the Kondo problem did in its day. The 
largest effects of such a defect arc naturally found in 
its immediate neighborhood, within a correlation length 



of order a few lattice spacings, such that local probes 
such as NMR, NQR, ^SR or STM arc particularly valu- 
able tools to investigate the fingerprints of correlations at 
small length scales. On the other hand, impurities also 
influence magnetic and transport properties of bulk ma- 
terials, and therefore to a lesser extent we also considered 
information provided by techniques which measure global 
responses of these systems (neutron scattering, photoe- 
mission, resistivity, Hall and optical conductivities, and 
Raman scattering). 

In a strongly correlated electron system, new length 
scales appear which reflect the interactions present, and 
which can be probed by impurities. For example, the 
magnetic response of an impurity in an interacting ID 
system is not only much larger than the usual Friedel re- 
sponse to an impurity in a free Fermi gas, but also mani- 
fests a spatial decay which reflects directly the magnetic 
correlation length ^ of the pure system. In an attempt 
to gain insight into the physics of correlated 2D systems, 
we have considered first the low energy states of ID sys- 
tems, where many models are exactly soluble (Section 

III) . There, a small change in a single hopping amplitude, 
or magnetic coupling or local potential is enough to pro- 
duce total backscattering in the spin or charge channels 
at T = (the two situations can be discussed separately 
because of spin-charge separation). Applying a uniform 
magnetic field if to a spin chain or ladder containing 
a non magnetic site yields a local alternating magnetic 
polarization with an amplitude proportional to H and a 
spatial extent generically of the order of the correlation 
length ^ of the pure system. 

Turning to the 2D (cuprate) case, experiments per- 
formed above Tc show that spinless defects substituting 
for Cu atoms in the planes produce a large alternating 
magnetic polarization on the sites surroimding the impu- 
rity. The size of the moment is orders of magnitude larger 
than that observed in conventional 3D metals with isova- 
lent spinless impurities, where the effect is attributable 
to simple Friedel-type oscillations in the LDOS (Section 

IV) . With the exception of the so-called spin glass phase, 
NMR experiments indicate that the moments are strictly 
paramagnetic, i.e. not frozen at zero field. The maximum 
amplitude of the staggered polarization is found to be of 
the form CH/ {T + <d), and the extention ^j^p of this po- 
larization away from the defect shows a similar T depen- 
dence. At optimal doping inYBCO, ^imp ~ 3 unit cells. 
This value is quite consistent with that deduced from INS 
measurements of Zn doped YBCO, and it compares well 
with the correlation length ^ of the pure system, as one 
would expect by analogy with the ID case. In 2D, rig- 
orous proofs of this notion are not available, but all ap- 
proximate methods, including strong-coupling RVB-type 
or nearly antiferromagnetic Fermi liquid scenarios are, 
in different regimes of the phase diagram, quite success- 
ful in producing magnetic and transport responses close 
to those found experimentally (Section V). This under- 
scores the importance of correlations in these materials, 
and the ability of impurity probes to reveal correlations 



58 



which can be hard to detect in the homogeneous state. 

The existence of a characteristic temperature scale Q 
in the magnetic response evokes Kondo physics, and we 
have argued that many of the experimental observations 
from NMR and transport properties display a consistent 
phenomenology with the picture of a screened moment 
formed around a nonmagnetic impurity. For increas- 
ing doping, beyond the optimally doped regime 6 gets 
larger and larger and exceeds 200K. As a result, both 
the amplitude of the alternating polarization and ^imp 
decrease. Similarly the large — log T upturns of the re- 
sistivity which initially scale with defect content in the 
underdoped regime progressively evolve upon overdop- 
ing towards the weak localization behavior expected for 
a Fermi liquid metal. 

Nevertheless, we have shown that a number of consid- 
erations invite caution against taking the Kondo analogy 
too literally. In particular, the evolution towards the 
overdoped regime, as the physics gets closer to the Fermi 
liquid limit, is yet to be understood. In a regular metal, 
a non magnetic impurity produces potential, rather than 
magnetic, scattering. Hence one does not expect mag- 
netic signatures, when interactions become weak enough, 
and for spinless impurities. NMR and /uSR suggest that 
the nature of the staggered magnetization evolves from 
a resonance to a bound state as hole doping is reduced, 
bringing the system close to the AF boundary. This be- 
havior is precisely that found in theoretical modelings of 
spinless impurities in the framework of t — J or Hub- 
bard models, which ignore processes leading to Kondo 
screening. Theoretical analysis based on slave bosons 
shows that the response of a strongly correlated system 
to a nonmagnetic defect is qualitatively similar in 2D 
and ID. Related theoretical treatments of transport ne- 
glecting Kondo screening also predict such upturns in 
the presence of strong Coulomb repulsions and suggest 
the need to be close to the AF instability. 

Ultimately, these questions regarding impurity states 
may not be resolved until consensus is reached on the 
peculiar normal state of the "pure" cuprates itself. We 
believe, nevertheless, that the evidence of the impurity 
experiments is providing important clues as to the na- 
ture of this state. For instance, while descriptions of 
the impurity-induced moment phenomenon from both 
the weak and strong coupling sides have proven some- 
what successful, the strength of the Coulomb repulsion 
U or exchange J has to be tuned to a value large enough 
to bring the system close to a magnetic instability. It is 
tempting to speculate that this is not an accident, but 
indicative of a renormalization group flow. In ID, where 
U simply controls the magnitude of the effect or rate of 
approach to the fixed point, the effect is clear, but in 2D 
the connection remains elusive. 

Impurity studies have been also shown to offer insight 
into fundamental aspects of the cuprate phase diagram 
and offer points of reference against which theories may 
be compared. For example, in the pseudogap phase, mag- 
netic properties of Cu sites located far from the impurity 



sites are quantitatively similar to those of the pure state. 
In particular, the opening of the pseudogap observed by 
NMR above Tc is not affected by spinless impurities. Fur- 
ther confirmation of this observation comes from resistiv- 
ity measurements of underdoped samples. 

In the superconducting state of the cuprates, a great 
deal was learned early on about the d-wave nature of the 
cuprate state from impurity substitutions (Section VI). 
This is because impurities in such states break Cooper 
pairs and mix order parameters with different signs on 
different parts of the Fermi surface, acting as a sim- 
ple phase-sensitive probe of order parameter symmetry. 
Thus nonmagnetic impurities depress the critical tem- 
perature of the cuprates and reduce the area of the su- 
perconducting dome. In the optimal to overdoped cases, 
several scenarios have been advocated to explain the ob- 
served linear decrease of Tc, down to zero, with the defect 
concentration, but a consensual theoretical description is 
still lacking. 

Single impurities in d-wave superconductors give rise 
to Andreev bound states, which overlap when sufficiently 
close, giving rise to a d-wave "impurity band". In the 
cuprates, NMR experiments show a paramagnetic re- 
sponse to the impurity qualitatively similar to the normal 
state, albeit with a smaller 6 at a given doping. However, 
due to impurity resonances, the E = response (LDOS 
effect) gives a significant contribution to the magnetiza- 
tion. At present the NMR experiments on Zn and other 
in-plane Cu subtituents arc not obviously consistent with 
STM spectroscopy on Zn impurities in BSCCO. On the 
one hand, the absence of a split LDOS peak is consistent 
with the paramagnetic (singlet) nature of the impurity 
state; on the other hand, the LDOS spatial pattern is at 
odds with what one would anticipate from NMR, namely 
a maximum LDOS on the near neighbors (and no LDOS 
on the impurity site) . 

In closing, we wish to stress again that our conclu- 
sions regarding impurity-induced magnetism, carefully 
developed in controlled situations in quasi-lD systems, 
as well as cuprate materials with impurity substitutions 
and electron irradiation, also have implications for the 
nominally pure cuprates. The same kind of analysis can 
be potentially applied to a wide class of other disordered 
strongly correlated systems, such as heavy fermions and 
2D organic conductors, provided materials problems are 
mastered to allow introduction of controlled defects. 

Generically, the structures of the materials we consider 
consist of active layers (chains, ladders, planes) sepa- 
rated from each other by buffer regions which may con- 
tain charge reservoirs. The thickness of the buffer region 
and the nature of the molecular bonding determine the 
strength of the couplings between the active layers. In 
quasi-lD materials it is difficult to change the carrier con- 
tent of the active layer, so that hydrostatic or chemical 
pressure are used to change the bandwidth and hence the 
ratio of kinetic to Coulomb energies. In quasi-2D cuprate 
compounds, it is usually possible to tune the carrier con- 
centration of the active layer through heterovalent charge 
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substitutions in the buffer layers. Thus, in most such ma- 
terials, the very process allowing the change of the ratio 
of kinetic to correlation energies brings in uncontrolled 
disorder into the system and affects magnetic and trans- 
port properties of strongly correlated systems, according 
to the ideas discussed here. 

We have shown that the magnetic effects due to a single 
nonmagnetic impurity are enhanced by increasing elec- 
tronic correlations and decreasing temperature, and that 
important new phenomena arise due to the interference 
of these states in dirty systems. These effects may ex- 
plain the obviously disordered "spin glass" phase of the 
underdoped cuprates, but could also under ly some of the 
observed behavior in the superconducting and pseudogap 
phase which have been previously attributed to the com- 
petition of homogeneous superconducting and antiferro- 
magnetic states, or to some more exotic cause. In the su- 
perconducting state, well-defined quasiparticles allow the 
most controlled theoretical treatments. Here, the effects 
of strong correlations and impurity-induced magnetism 
are small for optimal to overdoped systems, accounting 
for the success of the "dirty d-wave" approach, but they 
become more prominent in the imdcrdopcd regime. INS 
and ijSR experiments on YBCO suggest that the alter- 
nating polarization fluctuates at very low frequencies, but 
that disordered moments freeze in LSCO, BSCCO and 
other intrinsically disordered systems, creating glassy be- 
havior and, under some circumstances, ordered antifer- 
romagnetism coexisting with superconductivity. 

We believe that one should be careful not to ignore dis- 
order and inhomogeneity entirely but not to go too far 
either in assigning the origin of high-temperature super- 
conductivity itself to such effects. One can learn a great 
deal from studying the "dichotomy" between those ma- 
terials which are clean because they can be doped and 
rendered superconducting in stoichiometric form - pri- 
marily Y-123 and Y-248 and the majority of materials 
which are intrinsically disordered. These out of plane de- 
fects require much further study; they may act as simple 
scatterers and pairbrcakers, or they may influence the 
pair interaction in more subtle ways. Unquestionably, 
the phase diagrams of the nominally pure materials are 
strongly influenced by the type of doping and defects, 
and one important possible route to higher T^'s will be 
to understand and control the role of disorder. 
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